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A  COMPARISON  OF  DETERMINISTIC  LOT  SIZING  TECHNIQUES 
USING  FOCUS  FORECASTS  OF  STOCHASTIC  DEMAND  DATA 

CHAPTER  I 

INTRODUCTION 

.  \  *  r'  r.-r  . 

/ 

A  basic  concern  of  any  organization  which  manages 
production  or  inventory  is  the  question  '*How  much?^,  i.e., 
how  much  to  produce  or  how  much  (inventory)  to  order?  It  is 
a  very  easy  question  to  ask  but  not  quite  as  easy  to 
answer*  (Saunders,  1987) 

^The  difficulty  stems  from  the  nature  of  ^consumer" 
demand.  Specifically,  future  demand  is  seldom  known  with 
any  degree  of  certainty  { Tsadey 19£5a^ .  Anticipated  demand 
is  determined  as  best  as  possible  using  any  one  of  a 
multitude  of  forecasting  techniques  and  only  then  ^pluggedf"; 
into  a  production  lot  size  heuristic.  Unfortunately,  if  one 
subscribes  to  the  theory  that  forecasts  are  usually  wrong, 
then  the  old  adage,  "garbage  in,  garbage  out",  would  tend 
to  suggest  there  can  never  be  an  optimal  solution. 

^Most  research  in  this  area  has  therefore  concentrated 
on  developing  and/or  modifying  production  lot  size 
heuristics  in  the  hopes  of  providing  the  next  best  thing, 
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^east  wroni/**) 


i.e.,  the  --least  wrong*')  answer .  The  result  has  been  quite 
an  array  of  techniques  varying  in  both  size  (complexity) 
and  scope  (Ritchie  and.  Isado^-JLS&S^ — 

The  problem  left  to  industry  is  one  of  choice.  Which 
heuristic  is  best?  Several  studies  have  been  performed  in 
an  effort  to  answer  this  question  as  well*  Some  of  these  ' 
works  include  Benton  and  Whybark  (1982);  Callarman  and 
Hamrin  (1979,  1984);  De  Bodt  (1983);  De  Bodt  and  Wassenhove 
(1983a,  b);  Tsado  (1985a);  Ritchie  and  Tsado  (1986); 

Wemmerlov  and  Whybark  (1984).  Most  of  these  works, 
however,  deal  only  with  simulated  demand  data.  And  only 
Tsado  (1985a)  uses  empirical  demand  data  to  "validate"  the 
results  obtained  from  simulated  and  published  data. 
Unfortunately,  he  generates  the  forecasts  for  demand 
"artificially",  i.e.,  given  an  analysis  of  the  demand 

pattern  over  the  entire  demand  history. 

The  importance  of  validating  theoretical  results 
(either  analytical  or  simulated)  should  not  be 
underestimated.  For  example,  Flores  and  Whybark  (1986)  in 
their  study  of  forecasting  techniques  have  shown  that 
significant  differences  can  occur  between  the  results  found 
from  synthetic,  i.e.,  simulated,  data  and  those  obtained 
from  empirical  data.  They  state,  "...the  message  to 
researchers  rings  clear:  be  careful  in  drawing  "real-world" 


conclusions  from  laboratory  data."  Amar  and  Gupta  (1986) 
state  very  much  the  same  thing  regarding  their  study  of 
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simulated  and  empirical  demand  on  production  scheduling 
algorithms:  "Final  claims  about  the  superiority  of  [a] 
proposed  methodology...  can  be  settled  only  after 
sufficient  experience  with  real  life  situations." 

Further,  all  of  these  studies  combined  shortage  costs 
(if  included)  and  inventory  and  setup  costs  in  the  total 
cost  calculation.  As  shortages  may  imply  different  "costs" 
to  different  organizations,  it  would  be  interesting  to 
analyze  these  two  types  of  costs  separately.  The  need  for 
further  validation  of  previous  studies  of  lot  sizing 
techniques  is  therefore  justified. 


CHAPTER  II 


I-C.  O  ; 


LITERATURE  REVIEW 


A/  XT 

This  _chapter4'provides  a  review  of  the  literature  on 
basic  lot  sizing  techniques  and  their  app'  ication  to 
stochastic  demand.  Material  covering  the  lot  size 
algorithms  and  forecast  models  used  in  this  study  is 


presented  in  Chapters  III  and  IV,  respectively.  ^ 

/UcvAkiY  •.^•VkoA  (VVO'Uv'-,,  .. .  +I»l  .  ,  , 


,  !  ■  j 


Deterministic  Demand  Models  • 

In  order  to  determine  what  the  optimal  answer  is  to  the 
question,  "How  much?",  one  must  first  examine  the  type  of 
demand  to  be  modeled.  The  most  tractable  demand  model  is, 
of  course,  constant  or  level  demand.  Therefore,  if  the 
relevant  fixed  costs  (generally  order  or  setup  costs)  and 
variable  costs  (generally  inventory  holding  costs)  are 
known,  and: 

1.  demand  is  constant  and  deterministic, 

2.  the  order  quantity  is  assumed  a  continuous 
variable, 

3.  there  are  no  quantity  price  breaks, 

4.  costs  are  relatively  stable. 
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5.  replenishment/production  lead  time  is  zero,  and 

6.  no  shortages  or  back  orders  are  allowed, 

the  optimal  (most  economic)  production  or  order  quantity  is 
easily  derived  by  minimizing  the  total  relevant  costs  (TRC) 
per  unit  time,  i.e., 

TRC(Q)  =  variable  cost  +  fixed  cost  =  Q  h/2  +  A  D/Q 

where  Q  is  the  order  or  production  quantity,  h  is  the 
inventory  holding  cost  expressed  as  cost  per  unit  per 
period,  or  $/unit/period,  A  is  defined  as  the  fixed  set-up 
or  order  cost,  and  D  is  the  demand  rate  of  the  item  in 
units  per  unit  time  (from  Silver  and  Peterson,  1985). 
Specifically,  this  economic  order  quantity,  or  EOQ,  is 
given  as: 


EOQ  =  (2  A  D  /  h)* 

Once  we  relax  the  assumption  of  level  demand  and  allow 
time  variance,  however,  the  EOQ  is  no  longer  guaranteed  to 
provide  an  optimal  solution. 

This  assumption  of  constant  demand  is  one  of  the  first 
problems  we  encounter  with  the  basic  EOQ  model.  Few 
manufacturers,  suppliers,  or  retailers  can  expect  to  have 
requirements  for  exactly  N  units  of  a  product  every  period 
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over  the  entire  product's  life  cycle.  In  fact,  one  would 
expect  demand  to  (hopefully)  increase  from  zero  when  a 
product  is  introduced,  stabilize  once  initial  demands  are 
satisfied,  and  then  (unfortunately)  decrease  as  the  product 
becomes  obsolete.  And,  generally,  this  is  the  case. 

Hofer  (1977)  defines  the  fundamental  stages  of 
product/market  evolution  similarly.  However,  all  stages 
basically  fall  into  three  categories:  linearly  increasing 
demand,  level  demand,  and  linearly  decreasing  demand 
(Chalmet,  De  Bodt,  and  Van  Wassenhove,  1985).  These  time 
varying  levels  of  demand  render  the  once  optimal  EOQ  model 
to  the  level  of  a  mere  heuristic.  (A  heuristic  is  an 
algorithm  which  gives  near  optimal  problem  solutions.) 

Although  the  Wagner  and  Whitten  (1958)  dynamic 
programming  approach  to  the  time  varying  demand  model  is 
guaranteed  to  provide  an  optimal  solution,  many  authors 
feel  this  method  is  too  complicated  for  general  industrial 
use  (McLaren,  1977;  Silver,  1981;  Wagner,  1980).  In 
addition,  the  Wagner-Whitten  algorithm  may  provide 
sub-optimal  results  when  used  in  the  context  of  a  rolling 
demand  horizon  as  normally  used  in  industry  (De  Bodt  and 
Wassenhove,  1983a).  This  "sub-optimality"  results  from 
violation  of  the  assumption  that  demand  after  the  last 
period  in  the  horizon  is  zero.  It  is  interesting  to  note, 
however,  that  objections  to  use  of  the  Wagner-Whitten 
technique  have  steadily  declined  in  the  past  several  years 
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primarily  due  to  the  general  increase  in  the  power  of 
microcomputers  and  the  advent  of  such  efficient  high  level 
languages  as  C  and  PASCAL  (Saydam  and  McKnew,  1987  ;  Evans, 
1985) . 

As  a  result,  most  work  in  this  area  has  involved  the 
development,  test  and  evaluation  of  a  multitude  of 
(generally)  simple  heuristic  policies  in  an  attempt  to 
provide  solutions  as  "near-optimal"  as  possible.  Some  of 
these  heuristics  include  EOQ,  Silver-Meal,  part-period 
balancing  (Eisenhut),  and  least  total  cost  to  name  a  few. 

All  heuristics  can  be  divided  into  three  general 
classifications,  specifically: 

1.  EOQ  rules  (which  trade  off  order  costs  and 
holding  costs  per  unit  time,  e.g.,  discrete  EOQ), 

2.  Marginal  cost  rules  (which  equate  marginal 
order  and  holding  costs  per  period,  e.g.. 

Silver  and  Meal,  1973),  and 

3.  Target  rules  (which  set  holding  costs  equal  to  a 
target,  e.g.,  the  part-period  balance  algorithm 
which  increases  the  production  lot  size  until  the 
holding  cost  reaches  a  target  equal  to  the 
ordering  cost) 

(From  Wemmerlov  and  Whybark,  1984). 

Several  studies  have  shown  that  various  heuristics 
will  perform  differently  under  different  types  of 
conditions.  In  particular,  Ritchie  and  Tsado  (1986)  have 
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shown  the  best,  most  robust  lot-size  heuristics  to  be 
marginal  cost  (Groff,  1979),  the  simplified  part-period 
(balance),  and  Silver-Meal  algorithms.  Further,  they  show 
that  switching  from  one  rule  to  another  is  not  worthwhile, 
meaning  it  is  generally  better  to  use  a  good  rule  to  begin 
with. 

Therefore,  given  deterministic,  time  varying  demand, 
the  question  of  how  much  to  produce  or  order  seems  to  be 
relatively  easy  to  answer.  Unfortunately,  we  encounter  a 
second  problem  with  our  assumptions  on  demand. 
One  rarely  knows  with  any  certainty  what  future  demands 
will  be  (Tsado,  1985a). 

Stochastic  Demand  Models 

Solution  of  the  EOQ  or  lot  size  problem  for  stochastic 
demand  is  very  difficult.  As  a  result,  most  research  in 
this  area  has  focused  on  the  relative  performance  of 
algorithms  designed  for  deterministic  demand  as  applied  to 
stochastic  demand  over  a  rolling  horizon  (Benton  and 
Whybark,  1982;  Callarman, 1979 ;  Callarman  and  Hamrin,  1979, 
1984;  De  Bodt  and  Wassenhove,  1983;  Tsado,  1985a;  Wemmerlov 
and  Whybark,  1984). 

By  rolling  horizon,  we  mean  that  "lot  sizing  takes 
place  over  a  fixed  number  of  periods,  the  forecast  horizon, 
and  that  only  the  first  period's  decision  is  implemented. 
Next  period,  a  new  fixed  horizon  problem  is  made,  etc. 
(Baker  and  Peterson,  1979)." 


(From  Wemmerlov  and  Whybark, 
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1984.)  The  forecast  horizon,  in  turn,  is  based  on  mean 
time  between  orders,  or  TBO,  given  by: 

TBO  =  [ (2  A)/(D  h) ]2 

where  A  is  the  setup  or  holding  cost,  D  is  the  average 
forecasted  demand,  and  h  is  inventory  holding  cost.  TBO, 
therefore,  is  dependent  on  the  inventory  ratio,  i.e.,  A/h. 

There  is  no  one  consensus  on  what  values  of  A/h  to 
use.  In  a  problem  presented  by  Berry  (1972)  and  later  used 
by  many  others,  the  value  was  152.5.  De  Matteis  (1968), 
on  the  other  hand,  used  a  factor  of  100.  (From  Heemsbergen, 
1987.)  Tsado  (1985)  used  a  wide  range  of  values, 
specifically: 

7.5  10.0  30.0  50.0  70.0 
90.0  110.0  130.0  150.0  170.0 

Ritchie  and  Tsado  (1986)  used  a  value  of  400!  The  reasons 
why  the  literature  is  so  inconsistent  are  not  quite  clear, 
however  the  reasons  for  using  such  a  broad  range  are. 
DeBodt  and  Van  Wassenhove  (1983)  have  shown  that  various 
ranges  of  TBOs  will  lead  to  different  costs  even  when 
forecast  error  is  small.  Specifically,  smaller  values  of 
TBO  lead  to  larger  percentage  cost  increases.  Studies  by 
Blackburn  and  Millen  (1980)  suggest  that  a  forecast  horizon 
of  3  TBO  is  sufficient  to  minimize  cost  increases  due  to  a 
small  forecast  horizon,  but  only  for  heuristic  procedures. 
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Due  to  the  sensitive  nature  of  the  Wagner-Whitten 
algorithm,  Lundin  and  Morton  (1975)  suggest  a  forecast 
horizon  of  5  TBO  to  ensure  a  cost  performance  that  is 
within  1%  of  the  "optimal"  for  an  infinite  horizon 
(Wemmerlov  and  Whybark,  1984). 

Callarman  and  Hamrin  (1979),  in  one  of  the  earliest 
studies  of  stochastic  demand,  compared  the  relative 
performance  of  six  well-known  heuristics  (such  as  the  EOQ, 
part-period,  and  Silver-Meal  algorithms)  under  conditions 
of  uncorrelated  forecast  errors  and  fixed  lead  times  while 
using  the  coefficient  of  demand  variation  (s/m)  and  time 
between  orders  (TBO)  as  experimental  factors.  Their  basic 
conclusion  was  that  no  single  lot  sizing  rule  was  "best" 
under  all  conditions.  They  did  rank  the  heuristics, 
however,  with  Wagner-Whitten  coming  out  on  top,  followed 
closely  by  the  EOQ,  and  ending  up  with  Silver-Meal  as  one 
of  the  poorer  performers.  They  also  noted  that  total  costs 
tended  to  increase  with  forecast  error  resulting  in  smaller 
differences  in  the  relative  performance  of  the  heuristics. 
Callarman  (1979)  reaffirmed  these  results  for  an  inventory 
model  which,  unlike  the  previous  study,  explicitly  included 
stockouts.  (From  Tsado,  1985a,  and  Wemmerlov  and  Whybark, 
1984. ) 

Benton  and  Whybark  (1982)  confirmed  the  results  of 
Callarman  (1979)  and  Callarman  and  Hamrin  (1979). 
Specifically,  their  study  of  three  lot  sizing  techniques 
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(using  uncorrelated  forecast  errors  and  varying  such  system 
parameters  as  level  of  uncertainty),  showed  a  negative 
correlation  between  relative  heuristic  performance  (in 
terms  of  cost)  and  forecast  error,  i.e.,  as  forecast  error 
increased,  the  differences  in  cost  performance  of  the  three 
heuristics  decreased.  (From  Wemmerlov  and  Whybark,  1984.) 

De  Bodt  and  Wassenhove,  in  two  separate  studies, 
reported  findings  similar  to  those  previously  mentioned. 
Their  first  study  (1981)  examined  the  relative  performance, 
both  analytically  and  via  simulation,  of  the 
Wagner-Whitten,  Silver-Meal,  and  least  unit  cost 
algorithms.  Using  simulated  constant  demands  injected 
with  white  noise  and  forecasted  using  exponential 


smoothing. 

they 

showed  that 

cost 

differences  were 

negligible 

even 

when 

the  amount 

of 

forecast  error  was 

small.  Although 

the 

assumptions 

on 

demand  were  rather 

restrictive,  the  simulation  results  bore  out  the  analytical 
results  regarding  expected  cost  increases  due  to  forecast 
error . 

Their  second  study  (De  Bodt  and  Van  Wassenhove,  1983a) 
used  actual  demand  data  but  was  also  a  simulation  effort  in 
that  the  forecast  error  was  generated  artificially.  Their 
conclusions  were  essentially  the  same  (i.e.,  Silver-Meal, 
part-period,  least  unit  cost,  and  EOQ,  adjusted  to  cover 
integral  periods  of  demand,  performed  equally  well), 
however  they  did  state  a  preference  for  the  basic  EOQ  model 
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when  used  in  a  multi-stage  environment.  These  results  are 
interesting  in  that  some  of  the  operating  conditions  were 
different  from  those  used  in  earlier  research. 
Specifically,  the  authors  assumed  zero  lead  time  and  that 
the  forecast  error  for  the  next  period's  demand  was  zero 
(implying  zero  probability  of  a  stockout).  These 
assumptions  were  made  in  order  to  make  the  analytical  study 
tractable. 

The  following  year,  Wemmerlov  and  Whybark  (1984) 
presented  a  comprehensive  study  of  14  single-stage  lot 
sizing  techniques  using  demand  uncertainty  in  the  form  of 
forecast  errors  introduced  via  simulation.  The  operating 
conditions  are  similar  to  those  used  by  Benton  and  Whybark 
(1982),  Callarman  (1979),  and  Callarman  and  Hamrin  (1979). 
However,  they  do  incorporate  non-zero  lead  times.  Their 
most  important  results  are  quoted  as  follows: 

"1.  Relative  cost  performance  is  strongly  affected  by 
the  introduction  of  forecast  errors.  The  magnitude 
of  these  errors,  however,  is  not  significant. 

2.  The  Wagner-Whitten  procedure  loses  its  position  as 
the  least  cost  rule  (as  in  the  deterministic 
demand  model). 

3.  Only  two  rules,  [Wagner-Whitten]  and  WMR3 
[suggested  by  Wemmerlov  (1981)],  remain  on  the 
list  over  the  six  best  rules  overall  (from  the 
list  of  best  performers  in  the  deterministic 
demand  model). 

4.  The  relative  advantage  of  [Wagner-Whitten] 
compared  to  the  other  rules  decreases. 
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5.  The  performance  of  the  EOQ  rule  improves 
dramatically.  This  is,  no  doubt,  due  to  the 
non-discrete  character  of  this  rule,  leading  to 
the  ordering  of  a  larger  quantity  than  what  is 
needed  over  an  integer  number  of  periods.  In 
effect,  then,  the  EOQ  rule  carries  with  it  its  own 
safety  stock. 

6.  A  wider  choice  of  lot-sizing  rules  is  available 
when  compared  to  the  'no  uncertainty'  case.  Not 
only  are  there  no  differences,  from  a  statistical 
standpoint,  among  the  six  best  rules,  but  the  cost 
penalties  for  several  of  the  other  heuristics  are 
quite  small.  This  can  be  contrasted  to  the  case 
with  no  demand  uncertainty,  for  which  (Wagner- 
Whitten)...  emerge[s]  as  being  significantly 
better  than  the  other  rules." 


They  further  point  out  that  their  simulation  results 
seem  to  justify  current  industry  practice.  Specifically, 
Wagner-Whitten  is  not  applied  in  industry  (primarily  due  to 
complexity  and  "system  nervousness"  as  shown  in  their 
study).  The  EOQ  and  Eisenhut  algorithms,  on  the  other 
hand,  are  widely  used.  Previous  studies  involving 
deterministic  demand  would  lead  one  to  believe  this  is  bad 
practice.  However,  "if  it  is  acknowledged  that  the 
'forecasts  are  always  wrong’,  then  current  industry 
practice  seems  to  be  justified."  In  other  words,  the 
question  of  which  lot  sizing  technique  is  the  "best” 
becomes  moot;  a  simple  technique  will  probably  suffice. 

Tsado  (1985a)  concurs  with  the  results  of  Wemmerlov  and 
Whybark  (1984).  However,  he  recognizes  a  serious 
limitation  to  their  work  and  to  the  work  of  those  that 
preceded  them.  All  of  these  studies  involved  the  use  of 
simulated  demand  data  and/or  simulated  forecast  errors. 
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"While  simulation  is  an  important  tool  for  analysing 
problems  that  requires  (sic)  complex  mathematical 
solutions,  it  could  lead  to  different  results  by 
different  users,  if  there  are  differences  in  the  way 
the  data  was  simulated,  or  in  the  demand 
characteristics  of  the  data.  Moreover,  simulations 
cannot  always  explain  all  the  peculiarities  of  real 
life  experience." 


Additional  problems  with  the  previous  works  are: 

1.  Each  study  uses  different  versions  of  the 
part-period  balance  (see  Heemsbergen,  1987). 

2.  There  were  contradictions  in  some  of  the  initial 
assumptions  (discussed  previously). 

3.  All  forecast  errors  were  assumed  to  be  unbiased. 
However,  in  actual  practice  forecasts  may  be 
biased. 

4.  Forecasting  techniques  (exponential  smoothing, 
regression,  etc.)  were  not  used. 

Tsado's  (1985a)  study  therefore  attempted  to  examine 
the  possible  interactions  between  demand  pattern  and  lot 
size  performance,  lot-size  technique  and  forecasting 
technique  (or  forecast  parameters),  and  uncorrelated 
forecast  errors  and  lot-size  performance.  His  results 
follow: 

"1.  Forecast  errors  have  tremendous  influence  on 
the  performance  of  the  heuristic  policies 
even  when  these  forecast  errors  are  small. 

2.  With  the  exception  of  the  incremental  cost 
approach,  the  cost  differences  between  a 
number  of  heuristic  policies  is  small.  This 
contrast(s)  with  the  case  of  deterministic 
time  varying  demand  function  for  which  there 
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were  significant  differences  between  the 
performance  of  the  heuristic  policies. 

3.  The  magnitude  of  the  trend  in  demand  and  the 
type  of  forecasting  technique  used  seem  to 
have  insignificant  influence  on  performance...." 

Tsado  (1985a)  hoped  to  validate  the  previous  work  on 
lot-sizing  techniques  (as  well  as  compare  his  own  heuristic 
developed  specifically  for  stochastic  demand),  and,  on  the 
surface,  it  appears  that  he  did.  In  the  case  of  simulated 
data,  he  used  (linear)  exponential  smoothing.  He  broadened 
his  scope  on  the  published  data  by  including  the  use  of 
Winter's  seasonal  forecasting  model.  He  restricted  his  use 
to  (linear)  exponential  smoothing  once  again  for  the  actual 
demand  data.  His  reasoning  was  that  "these  forecasting 
methods  were  (appropriate)  because  they  provided  reasonable 
forecasts."  Unfortunately,  Tsado  limited  his  application 
of  these  forecasting  techniques  by  using  the  entire  demand 
history  available  to  him  to  fit  his  forecast. 

Industry,  on  the  other  hand,  does  not  have  this 
ability.  In  other  words,  forecasts  are  based  on  a  limited 
demand  history  (if  at  all)  and  are  then  updated 
continuously  over  the  "rolling  horizon",  i.e.,  from  period 
to  period.  Tsado' s  methods  therefore  do  not  seem 
reasonable.  What  methods,  then,  are  reasonable? 

Makridakis,  Andersen,  Carbone,  Fildes,  Hibon, 
Lewandowski,  Newton,  Parzen,  and  Winkler  (1982)  established 
the  following: 
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1.  Knowledge  of  the  underlying  demand  pattern  of  a 
time  series  does  help  in  choosing  a  model. 

2.  Simple  models  seem  to  work  well,  especially  when 
the  basic  series  is  changing  or  in  the  absence 
of  prior  knowledge  as  to  the  underlying 
structure  of  the  demand  pattern. 

3.  Under  the  conditions  where  simple  models  work  well, 
the  average  of  the  forecasts  from  several  simple 
models  was  superior  to  the  forecast  from  a  single 
model . 

Flores  and  Whybark  {1986),  however,  proposed  a 
different  method.  A  practitioner  developed  approach,  this 
method,  called  focus  forecasting,  involves  the  selection  of 
the  one  forecasting  model  which  would  have  performed  the 
best  in  the  recent  past  to  make  the  next  forecast.  As  a 
result  of  continuous  updating  of  all  forecasting  models, 
the  choice  of  forecasting  method  may  vary  from  time  to 
time. 

They  compared  both  techniques  (average  vs.  focus 
forecasting)  using  both  synthetic  (simulated)  and  empirical 
(actual)  demand  data.  The  method  of  averaging  performed 
best  on  the  synthetic  data.  More  importantly,  there  was  no 
significant  difference  in  the  relative  performance  of  focus 
forecasting  and  forecast  averaging  when  used  on  empirical 
data.  The  authors  believe  this  is  due  to  the  higher  mean 
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average  deviations  (MADs)  characteristic  of  actual  demand. 
In  other  words,  "empirical  time  series  are  far  more 
difficult  to  forecast  than  the  synthetic". 

It  is  important  to  note  that  none  of  the  seven 
forecasting  techniques  used  for  both  focus  forecasting  and 
forecast  averaging  used  any  form  of  regression,  exponential 
smoothing,  Winter's  method,  or  ARIMA  modeling  (although  a 
simple  moving  average  of  3  and  6  months  was  used).  They 
did,  however,  compare  the  focus  and  averaging  techniques 
with  exponential  smoothing  (as  a  common  basis  of 
comparison)  and  observed  that  exponential  smoothing 
generally  outperformed  both  although  the  significance  was 
not  as  great  for  the  empirical  demand  data.  It  would 
therefore  be  interesting  to  apply  focus  forecasting  to  the 
more  sophisticated  exponential  smoothing  models. 

Additionally,  all  of  the  aforementioned  studies 
combined  shortage  costs  with  inventory  holding  and  setup 
costs.  (Some  authors  incorporated  an  arbitrary  service 
level  using  a  predetermined  amount  of  safety  stock.) 
Wemmerldv  and  Whybark  (1984)  point  out  two  approaches.  One 
is  to  set  stockout  costs  as  a  separate  factor.  However, 
they  state  that  the  results  might  not  be  meaningful  when 
compared  to  other  studies.  The  other,  used  by  all  of  the 
studies  cited,  sets  service  levels  via  an  appropriate 
amount  of  safety  stock  and  quantifies  only  the  inventory 
holding  and  setup  costs. 
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A  basic  assumption  for  this  method  to  be  valid  is  that 
the  "cost"  of  a  stockout  is  the  same  for  everyone.  Another 
is  that  each  heuristic  employed  performs  the  same  for 
relative  shortage  costs  as  they  do  for  relative  inventory 
costs.  These  assumptions  may  not  be  valid,  therefore  it 
may  be  appropriate  to  look  at  these  two  types  of  costs 
separately. 

Bookbinder  and  H'ng  (1986)  employed  both  forecasting 
methodology  and  stockout  costs  (separate  from  the  standard 
inventory  costs)  in  their  study  of  rolling  horizon 
production  planning.  Their  main  emphasis,  however,  was  on 
the  procedure  for  probabilistic  production  planning  and  not 
on  the  relative  performance  of  the  lot  sizing  rule  employed 
within  the  procedure. 

And  finally,  the  baseline  used  in  the  previous  studies 
on  lot  size  algorithms  to  determine  relative  cost  increases 
for  each  heuristic  is  questionable.  Most  of  these  works 
(including  Wemmerlov  and  Whybark,  1984)  used  the 
WagnerWhitten  heuristic  (i.e.,  given  stochastic  demand)  as 
the  baseline.  Tsado  (1985a)  employed  the  EOQ  "heuristic". 
But  this  is  like  trying  to  measure  distance  with  a  rubber 
ruler ! 

It  is  suggested  that,  in  order  to  minimize  the  "error" 
inherent  in  such  an  approach,  the  optimal  Wagner-Whitten 
solution  given  a-priori  knowledge  of  the  demand  "history" 
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should  be  used.  The  increase  in  cost  due  to  the  heuristic 
using  forecast  demand  can  then  be  thought  of  as  the 
expected  value  of  perfect  information  (EVPI).  EVPI  may  be 
thought  of  as  the  maximum  amount  of  money  one  would  be 
willing  to  pay  for  perfect  knowledge  of  the  future.  This 
approach  should  not  only  minimize  the  error  in  the  design 
model,  but  should  also  be  intuitively  more  appealing.  (See 
Raiffa,  1968.) 

Research  Goals 

There  is  a  need  for  a  study  which  forecasts  empirical 
demand  in  the  same  manner  i**  hich  the  lot  sizing 
algorithms  are  implemented,  i.e.,  over  a  rolling  horizon. 
For  the  purposes  of  this  study,  a  system  of  focus 
forecasting  is  used. 

Further,  shortage  costs  need  to  be  analyzed  separately 
from  inventory  and  setup  costs  since  (1)  shortages  have 
a  "variable"  cost,  and  (2)  the  various  algorithms  may 
perform  differently  when  shortages  are  treated  as  a 
separate  entity. 

Finally,  validation  of  the  stochastic  heuristic 
presented  by  Tsado  (1985a,  b)  is  required. 


CHAPTER  III 


LOT  SIZE  HEURISTICS 

The  lot  size  heuristics  described  are  those  developed 
for  deterministic  or  discrete  demand.  As  stated,  there  are 
three  basic  categories  of  lot  size  algorithms:  EOQ-based, 
marginal  cost-based,  and  target-based.  The  algorithms  used 
in  this  study  for  each  category  are  the  standard  EOQ, 
Silver-Meal,  and  part-period  balance  methods,  respectively. 

Also  presented  are  Wagner  and  Whitten’s  (1958)  dynamic 
programming  method  which  is  used  as  both  an  optimal 
baseline  for  cost  comparison  and  as  a  separate  lot  size 
heuristic  (when  solved  for  forecast  demand)  and  Tsado's 
(1985a,  b)  stochastic  lot  size  heuristic.  Although  various 
refinements  exist  for  all  the  heuristics  listed,  the 
simpler  versions  were  used  in  the  study. 

Eisenhut 

The  part-period  balance  algorithm,  hereinafter  referred 
to  as  Eisenhut' s  lot  size  heuristic,  determines  the  number 
of  periods  to  order  or  produce  by  selecting  that  period  for 
which  holding  cost  most  closely  approximates  the  setup 
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or  ordering  cost. 

Using  the  example  provided  by  Silver  and  Peterson 
(1985),  let  the  setup  cost  =  $54,  holding  cost  =  $0.40  per 
unit  per  period,  and  demand  be  given  by  Di  =  {10,  62,  12, 

130,  154,  129}  for  the  first  6  periods.  The  algorithm 

yields 

T=l:  Holding  =  0 

T=2 :  Holding  =  D2h  =  $24.80  <  $52.00 

T=3 :  Holding  =  $24.80  +  2D3h  =  $34.40  <  $52.00 

T=4 :  Holding  =  $34.40  +  3D4h  =  $190.40  >  $54.00 

Therefore  since 

| 34 . 40  -  52 . 00  I  =  17.60  <  | 52 . 00  -  190. 40  |  =  138.40 

we  produce  for  3  integral  periods. 

EOQ 

The  economic  order  quantity,  derived  earlier,  is 
non-optimal  for  the  case  of  non-constant,  or  time-varying, 
demand.  To  be  used  as  a  heuristic  in  the  case  of 

stochastic  demand,  the  average  of  the  forecasted  demand  is 
used  in  the  model. 

Using  the  same  example  where  is  a  6-period 

forecast,  average  demand  is  approximately  83  units. 


Therefore 
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Q*  =  (2  x  52.00  x  83  /  0.40)*  =  147  units 

where  Q*  is  the  "optimal"  order  quantity  defined  by  the 
standard  EOQ  formula.  Notice  that  the  simple  form  of  EOQ 
provides  a  non-integer  time  supply  which,  in  effect,  acts 
as  an  automatic  safety  stock  (Tsado,  1985a). 

Silver-Meal 

Silver  and  Meal  (1973)  proposed  a  heuristic  for 
time-varying,  deterministic  demand  which  uses  the  concept 
of  marginal  cost,  i.e.,  it  attempts  to  minimize  total 
relevant  costs  per  unit  time  (a  quantity  which  we  will 
refer  to  as  TRCUT ) .  Expressed  as  a  function  of  time, 

TRCUT(T)  =  [A  +  T  ( ( t-1 )  Dt  h) ]  /  T 

where  A  is  the  setup  cost  and  t  ((t-1)  Dt  h)  is  the 
total  carrying  cost  to  the  end  of  period  T.  Selection  of 
the  "optimal"  number  of  periods  to  include  in  the 
replenishment  occurs  when  TRCUT(T+1)  >  TRCUT ( T ) .  Using  our 
previous  example: 

T=1 :  TRCUT ( 1 )  =  54.00/1  =  $54.00 

T=2 :  TRCUT ( 2 )  =  (54.00  +  1  x  62  x  0.04)  /  2  =  39.40 

T=3 :  TRCUT ( 3 )  =  (78.80  +  2  x  12  x  0.04)  /  3  =  29.47 

T=4 :  TRCUT ( 3 )  =  (88.40  +  3  x  130  x  .04)  /  4  =  61.10 
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and  we  select  a  replenishment  quantity  which  will  cover  3 
periods. 

Tsado 

Tsado's  (1985a)  stochastic  heuristic  is  primarily  a 
modification  of  the  EOQ  which  incorporates  the  idea  of 
minimizing  total  relevant  costs  for  a  given  replenishment 
cycle  while  keeping  track  of  previous  costs.  As  this 
method  is  generally  unknown,  more  will  be  said  regarding 
its  derivation. 

The  assumptions  used  in  the  derivation  are  (1)  no 
shortages  are  allowed,  (2)  demand  for  the  next  period  is 
known  with  certainty,  (3)  all  other  periods  are  forecast, 
(4)  a  replenishment  occurs  in  period  t  if  demand  cannot 
be  satisfied  for  period  t+1,  and  (5)  demand  is  assumed  to 
be  steady  and  continuous.  The  first  two  assumptions  are 
basically  equivalent  and  neither  is  used  in  this  research, 
i.e.,,  we  allow  shortages. 

Tsado  (1985a)  first  derives  an  equation  for  the 
expected  increase  in  stockholding  costs,  Stc,  given  that 
(1)  lead  time  is  zero,  (2)  replenishment  occurs 
instantaneously,  and  (3)  stock  at  the  end  of  the 
replenishment  interval  is  zero.  Specifically, 

Stc  =  h  D  L2  /  2 
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where  h  is  the  inventory  holding  cost  in  $/unit/period,  D 
is  the  rate  of  demand  (continuous),  and  L  is  the  length 
of  the  replenishment  interval.  Note  that,  although  the 
formula  is  derived  for  the  continuous  model,  the  heuristic 
is  applied  discretely,  i.e.,  to  periodic  demand. 

He  then  shows,  given  L  =  T  -  t  (since  we  wish  to 
satisfy  demand  up  to  the  horizon,  T),  that 

Stc  =  L2  D  h  /  2  =  (T  -  t)2  D  h/2 

where  Stc  and  L  are  as  previously  defined,  D  is  the 
current  forecast  of  demand  or  its  average,  h  is  as 
previously  defined,  T  is  the  last  period  in  the  forecast 
horizon,  and  t  is  the  period  of  the  present  setup. 

This  implies  that  total  relevant  costs  at  time  T  may 
be  written  as 

TRCUT(T)  =  [  Zt  +  (T  -  t)2  D  h/2  +  S  ]  /  T 

=  (Zt  +  S)/T  +  [  T  -  2  t  +  t2/T  ]  D  h/2 

where  Z^  is  the  total  inventory  cost  (holding  and  setup) 
up  to  time  t  and  S  is  the  fixed  cost  of  the  setup. 

Taking  the  derivative  with  respect  to  T, 

=  -  (Zt  +  S)/T2  +  (1  -  t2/T2)  D  h/2 


dTRCUT ( T ) /dT 
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which  set  to  zero  yields 

T  =  [  t2  +  2(Zt  +  S)/Dh  ] 2 

Since  L  =  T  -  t  and  the  replenishment  quantity,  Q, 
is  equal  to  the  average  forecast  of  demand,  D,  times  the 
replenishment  interval,  L,  Tsado's  lot  size  formula 
becomes 


Q  =  DL  =  D  (T-t)  =  DT  -  Dt  =  - ( Dt )  +  DT 
=  (-Dt)  +  [  (Dt)2  +  (2D  ( Zt  +  S))  /  h)2 

When  t=0,  this  equation  reduces  to  the  simple  EOQ 
formula,  therefore  the  first  setup  for  our  example  will  be 
identical  to  that  obtained  previously. 

Waqner-Whitten 

Wagner  and  Whitten's  (1958)  algorithm  is  a  dynamic 
program  which  provides  an  optimal  solution  to  the  discrete, 
time-varying  lot  size  problem.  When  used  as  a  heuristic 
for  the  stochastic  demand  model,  the  algorithm  computes  the 
"optimal''  solution  over  the  forecast  horizon  using  the 
forecasted  demand.  Using  our  example,  the  Wagner-Whitten 
procedure  and  solution  are  shown  in  Table  1. 

The  column  signifies  the  period  where  a  setup  occurs 
and  the  row  gives  the  current  period  we  are  looking  at. 
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1  N 

1 

2 

3 

4 

5 

6  1 

1 

52.00* 

2 

52.00 

24.80, 

76.80 

52.00 

52.00 

104.00 

3 

76.80 

9.60. 

88.40 

104.00 

4.80 

108.80 

76.80 

52.00 

128.80 

4 

88.40 

113.60 

128.80 

88.40 

156.00 

244.40 

104.00 

217.60 

52.00 

180.80 

52.00, 

140.40 

5 

140.40 

61.60 

202.00 

140.40 
52.00, 

192.40 

6 

192.40 

51.60, 

244.00 

192.4 

52.00 

244.40 

Table  1.  Wagner-Whitten  Procedure 


The  figure  in  row  1  and  column  1,  referred  to  as  (1,1), 
gives  the  total  cost  in  period  1,  i.e.,  the  cost  of  the 
setup. 

Cost  (2,1)  gives  us  the  total  cost  if  we  produce  enough 
in  period  1  to  cover  periods  1  and  2,  whereas  (2,2)  is  the 
total  cost  if  we  produce  in  both  periods.  The  first  number 
at  each  point  in  the  matrix  gives  the  total  inventory  cost 
(setup  and  holding)  from  the  previous  period.  The  second 
number  provides  the  inventory  holding  cost  of  the  current 
period  (if  we  had  produced  in  the  previous  setup  for  that 
period's  demand)  or  the  setup  cost  as  required. 
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The  asterisk  shows  the  "optimal"  policy  for  the  current 
period,  and,  except  for  the  policy  where  we  "test"  for  a 
setup  (which  is  given  on  the  diagonal),  that  cost  is 
carried  over  to  the  next  period  we  wish  to  look  at  (every 
(i,j)  "period"  where  j  <  i). 

The  Wagner-Whitten  theorem  states  that,  once  the 
optimal  policy  occurs  in  column  j,  calculations  may  be 
discontinued  for  any  column  i  where  i  <  j. 

In  general,  a  cost  C(i,j),  means  that  we  set  up  in 
period  j  to  produce  demands  for  periods  j,  j+1,...,  i, 
and  that  the  demands  for  periods  1,  2,...,  j-1  are 
produced  by  an  optimal  policy. 

Unlike  the  previous  heuristics,  the  Wagner-Whitten 
algorithm  must  be  used  over  the  entire  forecast 
(time-rolling)  horizon,  but,  like  the  others,  only  the 
first  period's  decision  is  implemented.  Cost  calculations 
are  carried  out  in  the  same  manner  as  in  the  other 
heuristics . 

The  "optimal"  Wagner-Whitten  solution  is  read  backwards 
through  the  matrix  using  the  "starred"  costs  as  a  guide. 
For  our  example,  we  show  setups  in  periods  1,  4,  and  5  and 
would  produce  in  the  first  period  for  periods  1,  2,  and  3. 

Note  that  optimality  of  the  Wagner-Whitten  procedure  is 
no  longer  guaranteed  since  the  assumptions  of  deterministic 
demand  and  zero  demand  after  period  T,  where  T  is  the 
planning  horizon,  are  violated. 


CHAPTER  IV 


THE  FORECAST  MODEL 

Time  series  analysis  involves  developing  forecasts  of  a 
variable  entirely  from  its  past  history.  These  techniques 
generally  model  the  variable  in  such  a  way  that  past 
patterns  in  the  data  series  are  used  to  help  modify  the 
mean  and  thereby  predict  future  values.  Although  past 
performance  is  no  guarantee  of  future  performance,  time 
series  methods  are  generally  successful  in  statistically 
stable  conditions,  for  short-term  forecasts  where  there  is 
insufficient  time  for  substantial  change  barring 
catastrophes  (as  in  our  current  study),  as  a  base  forecast 
for  judgemental  models,  and  for  screening  data  in  order  to 
better  understand  the  variable  being  forecasted  (Barron  and 
Targett,  1985). 

Introduction 

Jenkins  (1979)  describes  five  classes  of  time  series 
models.  These  are: 

1.  univariate  models  in  which  a  single  variable  is 
forecast  from  its  own  past  history, 
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2.  transfer  function  models  which  add  inputs  from 
other  variables, 

3.  intervention  models  which  represent  unusual  events 
such  as  strikes,  etc., 

4.  multivariate  stochastic  models  which  represent 
several  series  with  mutual  interaction,  and 

5.  multivariate  transfer  function  models  which 
relate  several  output  variables  to  several  input 
variables  in  which  a  relationship  exists. 

Univariate  models,  although  of  an  elemental  nature,  are 
important  from  a  forecasting  viewpoint  for  several  reasons. 
First,  they  may  be  the  only  model  which  is  the  only 
practical  approach  based  on  the  magnitude  of  the  problem. 
Second,  ''it  may  be  impossible  to  find,  or  there  may  not 
exist,  variables  related  to  the  one  being  forecast.  Third, 
when  multivariate  models  exist,  the  univariate  model  may  be 
used  as  a  baseline  to  measure  the  other's  performance.  And 
finally,  the  presence  of  large  residuals  (the  difference 
between  actual  values  and  the  "stationary  mean")  may 
correspond  to  strikes,  faulty  data,  etc.  and  therefore  act 
as  a  tool  to  screen  data.  In  spite  of  these  points. 
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however,  it  must  be  recognized  that  univariate  models  are 
generally  valid  for  short-term  forecasts  only  (Ibid.)*  As 
all  of  the  studies  mentioned  in  Chapter  II  utilized 
univariate  forecasting  methodology,  this  discussion  will 
be  restricted  to  procedures  in  this  area. 

Univariate  models  are  classified  by  Barron  and  Targett 
(1985)  according  to  the  type  of  series  to  which  they  can 
be  applied.  These  are: 

1.  stationary  (random  variation  about  a  mean  or  a 
series  which  may  be  modeled  as  a  stochastic 
"random  walk"), 

2.  trending  (a  consistent  movement  either  upwards  or 
downwards  in  the  series), 

3.  seasonal/cyclical  (a  series  which  exhibits  a 
pattern  over  a  number  of  time  periods  where 
seasonal  implies  a  period  of  a  year  or  less  and 
cyclical  refers  to  a  pattern  greater  than  one 
year),  and 

4.  seasonal  and/or  cyclical  with  a  trend  (a  complex 
of  seasonal  and/or  cyclical  patterns  and  trends). 
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The  last  three  classifications  may  be  grouped  under  the 
heading  of  non-stationary  time  series.  Johnson  and 
Montgomery  (1976)  state  that  the  basic  goal  in  any 
univariate  time  series  method  is  to  reduce  the  residuals, 
or  error,  to  a  normally  distributed  random  variable  with 
mean  zero  and  constant  variance  (also  known  as  "white 
noise").  In  other  words  we  seek  a  stationary  model  from  a 
non-stationary  time  series.  (See  Hoel,  Port,  and  Stone, 
1972. ) 

There  are  two  general  types  of  time  series  forecast 
methods,  those  involving  smoothing  techniques  and  those 
involving  autoregressive  parameters,  generally  referred  to 
as  ARMA  (p,q)  models.  Johnson  and  Montgomery  (1976) 
suggest  that  ARMA  models  should  be  considered  only  when 
there  exists  a  sufficient  amount  of  demand  history  for 
analysis,  typically  around  36  periods  or  more.  Since  large 
amounts  of  demand  history  from  the  same  environment  may  not 
always  be  available  and  previous  studies  have  not  utilized 
ARMA  models,  ARMA  models  were  not  considered  for  use  in  the 
current  study. 

Further,  since  lot  sizing  is  performed  on  a  rolling 
horizon,  forecasting  should  be  performed  on  the  same  basis. 
Therefore,  we  will  restrict  ourselves  to  exponential 
smoothing  models  when  applied  to  the  concept  of  focus 
forecasting. 
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Exponential  Smoothing 

Smoothing  techniques  (or  models)  replace  the  original 
time  series  by  a  "smoothed”  one,  i.e.,  one  produced  from 
statistical  or  weighted  averages  of  values  from  the 
original  series  in  an  attempt  to  reduce  or  discount  the 
random  fluctuations  or  variance.  Generally,  the  last 
smoothed  value(s)  provide(s)  the  forecast  for  all  future 
time  periods  in  the  (rolling)  forecast  horizon. 

Simple  Exponential  Smoothing  Model 

The  simplest  case  is,  of  course,  when  the  time  series 
is  already  stationary,  i.e.,  it  may  be  represented  by 

xt  =  m  +  et 

where  m  (or  mu)  is  the  statistical  mean  of  the  time 
series  and  et  is  the  error  or  difference  between  the  mean 
and  the  actual  value  of  the  data  point.  Two  techniques 
which  deal  with  such  stationary  models  are  moving  averages 
(not  discussed)  and  simple  exponential  smoothing. 

Exponential  smoothing  assumes  that  recent  data  is  more 
important  than  old  data;  a  concept  which  is  rather 
intuitively  appealing.  Then,  based  on  the  relative  value 
attached  to  the  significance  of  the  residuals,  it  computes 
a  smoothed  "average"  of  the  data.  Specifically,  the  model 


states 
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St  ~  ( 1  "  3 )  St— 1  +  3  Xt 

where  S*-  is  the  new  smoothed  value  at  time  t,  St-i  is 
the  old  smoothed  value  at  time  t-1,  xt  is  the  most 
recent  actual  value,  and  a  (or  alpha)  is  a  weight  chosen 
by  the  forecaster  such  that  0  <  a  <1.  Obviously,  the 
larger  the  value  of  alpha  the  more  weight  will  be  attached 
to  the  most  recent  data  point. 

To  see  this,  one  merely  needs  to  expand  the  equation 
for  all  "N"  which  yields 

St  =  a  xt  +  a  (1-a)  xt-i  +  a  (1-a)2  xt-2  +  . .. 

where  the  weights  given  to  the  data  points  from  the  most 
recent  to  the  most  distant  are  a,  a  (1-a),  a  (1-a)2,  and 
so  forth  (Barron  and  Targett,  1985).  Since  both  a  and 
(1-a)  are  less  than  one,  the  weights  are  decreasing 
monotonically  with  time. 

Silver  and  Peterson  (1985)  rewrite  the  exponential 
smoothing  model  to  obtain 

St  =  St-i  +  a  (x(t)  -  St-!)  =  St-i  +  a  et 

where  all  variables  are  as  previously  defined.  This 
implies  the  new  forecast  value  is  equal  to  the  old  forecast 
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value  minus  a  fraction  of  the  most  recent  error.  In  other 
words,  the  exponential  smoothing  model  assumes  that  a 
portion  of  the  last  forecast  error,  namely  (1-a)  is  due  to 
some  random  fluctuation  and  the  other  portion,  namely 
alpha,  is  due  to  some  real  shift  in  the  value  of  the 
estimate.  In  practice,  the  value  of  alpha  usually  ranges 
between  0.1  and  0.4  (Ibid.). 

Holt's  Exponential  Smoothing  Model 
Now  consider  time  series  which  are  initially 
non-stationary  but  which  can  be  made  stationary  by 
differencing.  By  differencing  we  mean 


DEL  St  =  St  - 

Sf-d 

where 

DEL 

is  the  differencing 

operator 

and  d  is  the 

period 

Of 

differencing.  For  a 

strictly 

trending  time 

series. 

a 

difference  of  d=l  will 

yield  a 

stationary  time 

series . 

To 

see 

this,  one  must  first 

examine 

the  statistical 

significance  of  trending  and  seasonal  data.  When  a  time 
series  trends,  the  values  between  successive  data  points 
are  highly  correlated.  (Since  the  time  series  is 
correlated  with  itself,  a  more  appropriate  term  is 
"autocorrelated" . )  The  same  is  true  for  seasonal  time 
series  where  the  autocorrelation  occurs  at  lag  d,  i.e., 
for  time  series  values  S^  and  St-^.  The  differencing 
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operator  therefore  yields  white  noise,  i.e.,  the  residuals 
are  normally  distributed  with  zero  mean  and  constant 
variance  (Jenkins,  1979).  As  mentioned  earlier,  the 
objective  of  all  time  series  analyses  is  to  fit  a  model 
such  that  the  residuals  yield  white  noise  (Johnson  and 
Montgomery,  1976). 

Due  to  the  nature  of  how  the  forecasting  methodology 
was  implemented,  exponential  smoothing  models  which  account 
for  seasonality  were  ignored.  Trend,  however,  is  accounted 
for  through  the  use  of  Holt's  exponential  smoothing  model. 
(Linear  regression  may  also  be  used  but  was  not  employed  in 
this  study.) 

A  strictly  trending  (linear)  time  series  will  take  the 

form 


xt  =  m  +  A  t  +  et 

where  lit  defines  a  linear  trend  (as  a  function  of  t) 
with  slope  A.  Other  trends  are  possible.  However,  our 
discussion  is  limited  to  linear  trends.  Successfully 
differencing  a  time  series  more  than  once  for  trend  is  a 
good  indication  the  trend  is  non-linear  (Ibid.). 

Let  the  trend  at  time  t  be  given  by  Tt  =  St  -  St-i» 
Since  St  is  a  random  variable,  Tt  is  also  a  random 
variable.  Therefore,  using  the  same  logic  as  simple 
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exponential  smoothing,  we  can  smooth  the  trend  by  the 
following: 


Tt  =  (1  "  9)  Tt-1  +  9  <xt  “  xt-l)) 

i.e.,  the  smoothed  trend  is  equal  to  a  portion  of  the 
previous  smoothed  trend  plus  a  portion  of  the  most  recently 
observed  trend.  The  selection  of  g  (or  gamma)  is  made  in 
the  same  manner  as  alpha. 

Using  this  estimate,  we  can  modify  St-i  in  the  simple 
exponential  smoothing  model  to  obtain 

St  =  (1  -  a)  (St-i  +  Tt-i )  +  a  x(- 


or  more  generally 


Ft+i  =  St  +  i  Tt 

where  Ft+i  is  the  forecast  for  the  t+i'th  period 
(Barron  and  Targett,  1985) 

Focus  Forecasting 

Flores  and  Whybark  (1986)  studied  two  forecasting 
systems,  "one  recommended  by  practitioners  for  use  in 


inventory  management,  and  the  other  the  result  of  an 
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international  forecasting  competition  among  academics." 
These  are  the  methods  of  focus  forecasting  and  forecast 
averaging,  respectively. 

Although  this  topic  was  touched  upon  briefly  in  the 
literature  review,  we  would  like  to  say  a  little  more 
regarding  the  aforementioned  study  and  our  proposed 
extension. 

The  forecast  procedures  used  in  the  comparison  were 
very  simplistic,  e.g.,  "the  forecast  for  the  next  month  is 
the  actual  demand  for  the  same  month  last  year....  [or] 
...is  one-sixth  of  the  total  actual  demand  for  the  last  six 
months  (a  two-quarter  moving  average)."  Another,  slightly 
convoluted  approach  was  "if  the  demand  in  the  last  six 
months  is  more  than  2.4  times  the  demand  for  the  six  months 
preceding  that,  the  forecast  for  the  next  month  is 
one-third  of  the  demand  for  the  same  three  month  period 
last  year  (i.e.,  we  are  starting  into  the  downside  of  a 
seasonal  swing) . " 

The  focus  and  averaging  techniques  were  then  compared 
to  each  other  and,  most  importantly,  to  exponential 
smoothing,  i.e.,  exponential  smoothing  provided  the 
"baseline"  for  comparison. 

Although  averaging  performed  better  than  focus 
forecasting  on  the  simulated  data  (there  was  no  statistical 
difference  for  the  empirical  data),  neither  procedure 
performed  better  than  exponential  smoothing.  In  fact. 
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exponential  smoothing  was  significantly  better  than  either 
of  the  other  two  procedures. 

Exponential  smoothing  would  then  seem  to  be  the  obvious 
choice.  The  next  question,  however,  is  the  selection  of 
alpha  and  gamma,  i.e.,  the  forecast  parameters.  Past 
studies  have  "fit”  the  parameters  over  the  entire  demand 
history  of  each  empirical  data  set  (when  empirical  data 
were  used).  Industry,  of  course,  doesn't  have  this  type  of 
clairvoyance;  they  would  have  to  take  an  educated  guess 
given  a  limited  demand  history  and  monitor  the  forecast 
model  to  make  appropriate  changes  when  necessary.  But, 
since  this  study  was  "automated",  we  did  not  have  this 
"luxury"  either. 

It  therefore  makes  sense  to  either  (1)  average  the 
exponential  smoothing  forecasts  from  varying  parameter 
levels  or  (2)  use  the  focus  forecasting  approach.  We 
selected  the  focus  forecasting  approach  as  it  seems  to  be 
the  most  appealing  (intuitively).  The  idea  is  to  keep 
track  of  the  mean  absolute  deviation  and  bias  of  a  set  of 
exponentially  smoothed  forecasts  and  select  the  one  best 
forecast  for  the  next  planning  horizon. 

Silver  and  Peterson  (1985),  however,  argue  that 
changing  the  smoothing  constants  (what  they  refer  to  as 
"adaptive"  smoothing),  while  having  considerable  intuitive 
appeal,  is  "not  necessarily  better  than  regular. 
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non-adaptive  smoothing."  (See  Ekern,  1981;  Flowers,  1980; 
and  Gardner  and  Dannenbring  (1980))  Specifically,  they 
feel  the  resulting  forecasts  would  be  excessively 
"nervous" . 

Fortunately,  the  lot  sizing  problem  only  requires  use 
of  an  extended  forecast  about  every  "TBO"  periods.  For  our 
purposes,  the  focus,  or  adaptive,  approach  should  be  quite 
reasonable.  In  fact,  comparisons  of  the  mean  absolute 
deviation  (MAD)  for  the  adaptive  procedure  to  the  MADs  of 
each  individual,  static  procedure  (tested  during  program 
development)  were  quite  favorable  and  tend  to  support  this 
position. 

Separate  research  regarding  the  relative  merit  of  focus 
or  "adaptive"  and  averaged  exponential  smoothing  techniques 
(as  used  in  automatic  forecasting)  is  probably  warranted. 


CHAPTER  V 


THE  EXPERIMENT 

This  chapter  constitutes  the  bulk  of  this  research  and 
is  divided  as  follows:  Sample  Data,  Assumptions, 
Performance  Criteria,  Computer  Model,  Experimental  Design, 
Results,  and  Analysis.  Concluding  remarks  are  contained  in 
Chapter  VI. 

Sample  Data 

Data  was  obtained  from  two  separate  industrial  sources. 
The  first  group  originally  contained  500  data  sets 
consisting  of  52  weekly  periods,  however  only  207  of  these 
proved  suitable  for  our  purposes.  Specifically,  all  data 
sets  which  contained  an  alphanumeric  or  zeroes  were 
discarded.  The  second  group  was  very  limited  at  5  data 
sets,  however  each  data  set  consisted  of  78  (monthly) 
periods. 

From  the  first  group  of  207,  36  were  selected  randomly 
for  the  study.  They  were  then  classified  according  to  the 
coefficient  of  variation  and  data  type. 

The  variability  of  each  data  set  was  determined  to  be 
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either  low  (0  <  s/m  <  0.5),  medium  (0.5  <  s/m  <  1.0),  or 
high  (s/m  >  1.0),  where  s  (or  sigma)  is  the  standard 
deviation.  Selection  of  the  "cut-offs"  were  arbitrary. 

Data  type  consisted  of  two  classifications:  linear  and 
non-linear.  The  reason  for  this  was  two-fold.  First,  a 
relatively  small  number  of  data  sets  were  selected  for  the 
study.  Second,  the  forecast  model  used  simple  exponential 
smoothing  as  well  as  Holt's  exponential  smoothing  model  for 
a  linear  trend,  i.e.,  it  wasn't  "designed"  to  handle  a 
non-linear  demand  pattern.  It  was  therefore  necessary  to 
account  for  the  possibility  the  forecast  model  would 
perform  worse  for  the  non-linear  case. 

An  ARIMA  "identify"  was  performed  on  each  data  set 
using  the  Statistical  Analysis  System  (SAS)  ((c)  1985  by 
SAS  Institute  Inc.)  in  order  to  determine  which  demand 
patterns  could  not  be  considered  level,  i.e.,  as  white 
noise.  A  second  "identify"  using  a  differencing  of  1 
determined  which  demand  patterns  could  be  considered 
non-linear . 

A  complete  listing  of  the  36  demand  patterns  in  the 
first  group  is  given  in  Appendix  A,  however  we  will  discuss 
a  few  selected  patterns  here. 

Figure  1  shows  a  plot  of  the  first  data  set.  Although 
the  demand  series  is  generally  linear  with  a  relatively 
small  variance,  large  outliers  occurring  at  periods  14  and 
16  "inflate"  the  coefficient  of  variation  to  1.49.  Outliers 
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such  as  these  posed  somewhat  of  a  problem  for  the  automatic 
forecast  model;  a  procedure  discounting  such  outliers  was 
developed  and  will  be  discussed  at  a  later  point  in  this 
chapter. 

The  demand  series  depicted  in  Figure  2  also  constitutes 
white  noise,  however  "spikes"  occur  at  periods  2,  16,  32, 
33,  and  46.  If  one  considers  periods  32  and  33  to  be 
"split",  then  the  spikes  occur  about  1  every  15  cycles. 
Like  the  previous  demand  series,  these  spikes  inflate  the 
coefficient  of  variation  to  about  1.68. 

Figure  3  on  the  other  hand  shows  no  significant  spiking 
when  compared  to  the  general  variability  of  the  data  set. 
This  demand  pattern  qualified  as  white  noise  and  showed  a 
moderate  coefficient  of  variation  of  about  0.49. 

A  demand  series  showing  a  slight  downward  trend  (after 
an  initial  upswing)  and  moderate  variance  (coefficient  of 
variation  of  approximately  0.33)  is  shown  in  Figure  4. 
Figure  5  shows  a  demand  series  with  a  definite  drop  in 
demand  in  period  8  followed  by  an  upward  trend.  Coefficient 
of  variation  for  this  series  is  slightly  higher  as  a  result 
(about  0.53).  Both  series  are  considered  linear  (non¬ 
constant)  . 

Non-linear  demand  sets  are  given  in  Figures  6  and  7. 
Figure  6  shows  a  rough  "concave"  pattern  which  is  somewhat 
obscured  in  spite  of  the  relatively  low  coefficient  of 
variation  (0.35).  The  non-linear  pattern  of  Figure  7,  on 
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Figure  5 
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Non-Linear  Demand 

Example  1 
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Figure  6 


49 


Non-Linear  Demand 

Example  2 
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Figure  7 
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the  other  hand,  is  clearly  convex  and  could  very  well  be 
seasonal  (although  we  can't  be  sure  due  to  the  limited 
history).  The  coefficient  of  variation  for  this  set  is  a 
slightly  higher  0.46. 

Table  2  shows  the  exact  break-out  of  the  36  demand  sets 
in  this  group. 


|  Structure 

Coef  of  Var  j 

Linear 

|  Non-Linear 

Low  j 

12 

1  8 

Medium  | 

6 

1  5 

High  | 

5 

1  o 

Table  2.  Data  Classification  —  Group  1 

The  data  in  the  second  group  is  given  in  Appendix  B. 
Due  to  the  limited  number  of  data  sets  in  this  group,  they 
were  not  classified  by  data  type  or  degree  of  variation.  It 
should  be  noted,  however,  that  significantly  more  variation 
in  structure  can  exist  for  these  longer  demand  sets  as  they 
span  a  period  of  over  5  years. 


Assumptions 

Assumptions  used  to  develop  the  single-stage, 
production  lot  size  problem  are  similar  to  those  employed 
by  other  research  and  are  as  follows: 
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1.  Demand  is  probabilistic  and  is  forecast  using  a 
limited  amount  of  prior  history. 

2.  A  fixed  cost  is  incurred  for  each  setup. 

3.  The  inventory  holding  cost  is  a  function  of  the 
amount  of  inventory  on  hand  at  the  end  of  a  given 
period. 

4.  Production  lead  time  is  zero  (i.e.,  we  have  enough 
inventory  at  the  end  of  a  production  period  to 
meet  that  period's  demand). 

5.  All  demands  are  met  at  the  end  of  each  period. 

6.  There  is  no  safety  stock  except  that  inherent  in 
a  particular  lot  size  heuristic. 

7.  Back  orders  are  allowed. 

8.  There  is  no  monetary  penalty  for  shortages  in  the 
cost  calculations  (i.e.,  shortages  are  handled  as 
a  separate  criteria). 

9.  Demand  for  the  next  period  is  not  known  with 
certainty. 
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10.  An  updated  forecast  is  available  for  any  period. 

Performance  Criteria 

There  are  three  types  of  criteria  (dependent 
variables)  of  interest  in  this  study.  These  are  cost, 
number  of  stockouts,  and  the  amount  short  per  stockout. 

Relative  Cost 

Previous  studies  have  used  the  Wagner-Whitten  procedure 
when  used  as  a  heuristic  as  the  baseline  for  cost 
comparisons.  Unfortunately,  the  Wagner-Whitten  procedure 
is  suboptimal  in  the  case  of  a  rolling  horizon  and 
probabilistic  demand. 

Arguments  for  the  use  of  the  Wagner-Whitten  heuristic 
as  the  baseline  are: 

1.  Wagner-Whitten  is  the  baseline  used  for  the 
deterministic  case. 

2.  It's  not  known  before  hand  which  rule  will 
outperform  the  others. 

3.  Use  of  the  Wagner-Whitten  "heuristic"  will  make 
the  study  more  easily  comparable  to  previous 


works. 
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We  feel  these  reasons  do  not  justify  the  use  of  one 
heuristic  as  a  basis  of  comparison.  It's  true  that  we  do 
not  know  what  the  optimal  cost  of  a  probabilistic  lot  size 
problem  will  be  until  the  demands  have  already  been 
satisfied,  i.e.,  we  don't  know  what  our  future  demands  will 
be.  However,  by  comparing  the  cost  obtained  through  the 
use  of  a  heuristic  (when  demand  is  considered  stochastic) 
with  the  optimal  cost  obtained  by  Wagner-Whitten  over  the 
entire  demand  "history"  (when  considered  deterministic),  we 
obtain  a  true,  fixed  reference  for  comparison. 

The  key  is  the  interpretation  of  the  cost  comparison. 
Specifically,  this  difference  in  cost  may  be  thought  of  as 
the  maximum  amount  of  money  we  would  be  willing  to  pay  for 
perfect  knowledge  of  our  future  demand  (referred  to  as  the 
expected  value  of  perfect  information  or  EVPI).  (See 
Raiffa,  1968.) 


Number  of  Stockouts 

Wemmerlov  and  Whybark  (1984),  Tsado  (1985a),  and  others 
arbitrarily  set  service  levels  in  order  to  handle  the 
question  of  stockouts.  By  service  level,  we  mean  that 
there  exists  enough  safety  stock  to  assure  demands  are  met 
at  least  percent  of  the  time.  Generally,  levels 
between  90  and  99.999  percent  have  been  chosen.  As  a 
result,  the  stockout  question  is  largely  ignored. 

Since  we  assume  that  stockouts  have  a  "variable"  cost, 
i.e.,  the  cost  of  a  stockout  to  one  organization  may  be 
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quite  less  than  that  perceived  by  another,  setting  an 
arbitrary  service  level  may  not  be  appropriate.  Further, 
by  pre-determining  a  service  level,  the  effects  of  a  lot 
size  algorithm  on  inventory  (holding  and  setup)  costs  and 
stockout  costs  may  be  confounded. 

In  a  manner  similar  to  that  employed  by  Bookbinder  and 
H'ng  (1986),  we  chose  to  "count”  the  number  of  times  a  lot 
size  heuristic  produced  a  stockout.  Obviously,  this  number 
will  vary  according  to  the  TBO  level,  therefore  we  chose  to 
compute  the  stockout  "cost"  as  the  number  of  times  a 
stockout  occurred  expressed  as  a  percentage  of  the  number 
of  replenishments  made. 

For  example,  given  a  52  period  demand  "history"  with  a 
TBO  level  of  2,  then  5  stockouts  out  of  26  replenishments 
(approximately)  will  yield  a  stockout  "cost”  of  0.1923, 
i.e.,  about  19.23%  of  the  replenishments  made  experienced  a 
stockout.  For  a  TBO  of  6,  5  stockouts  would  imply  a  "cost" 
of  57.69%. 


Percent  Short  per  Stockout 

Another  factor  in  the  stockout  question  is  the  amount 
of  shortage  when  a  stockout  occurs.  The  average  number 
short  per  stockout  is  therefore  an  important  "cost" 
consideration,  however,  an  average  shortage  of  N  units 
doesn't  tell  us  much. 

There  are  two  ways  of  handling  this  problem.  One  is  to 
express  the  shortage  as  a  percentage  of  average  demand. 
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another  is  to  express  it  as  a  percentage  of  the  demand  for 
the  period  in  which  we  were  short.  We  chose  the  later. 

Justification  for  our  selection  is  as  follows.  Consider 
an  average  demand  of  500  units.  If  we  were  to  have 
forecast  a  demand  of  550  units  where  actual  demand  was  600 
units,  then  our  percentage  short  is  only  8.3%  of  the  actual 
demand.  If  we  had  used  average  demand,  we  would  have  shown 
a  shortage  of  10%.  Now  assume  an  average  demand  of  50 
units.  Similarly,  assume  a  forecasted  demand  of  100  units 
and  an  actual  demand  of  150  units.  But  now  we  show  a 
shortage  of  33.3%  of  actual  demand  and  a  misleading  100%  of 
average  demand. 

In  both  cases  the  forecast  was  50  units  greater  than 
average  demand,  and  actual  demand  was  50  units  greater  than 
forecasted  demand.  Obviously,  shortage  "cost"  expressed  as 
a  percentage  of  actual  demand  is  a  more  accurate  estimate 
of  the  true  "cost"  associated  with  a  shortage. 

Computer  Model 

Although  not  a  simulation  study,  a  computer  model  was 
used  to  generate  forecasts,  compute  production  policies  via 
the  various  lot  size  heuristics  (including  the  optimal 
Wagner-Whitten  cost),  and  to  compute  the  costs  associated 
with  each  policy.  This  section  discusses  the  issues  of 
program  development  and  validation. 
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Program  Development 

Both  the  forecasting  procedure  and  lot  size  procedure 
were  automated  via  a  program  written  in  MICROSOFT 
QuickBASIC  (R)  and  run  on  an  IBM  XT  (R)  compatible 
microcomputer.  The  program  listing  is  given  in  Appendix  C. 

For  purposes  of  clarity,  this  section  is  further 
subdivided  into  2  groups.  The  first  discusses  the  forecast 
algorithm;  the  second  addresses  the  lot  size  procedure. 

The  Forecast  Procedure.  The  complete  forecast  is 
generated  over  the  entire  demand  history  of  each  data  set 
(on  a  rolling  horizon  basis)  prior  to  implementation  of  the 
lot  size  procedure.  Estimates  of  level  demand  and  trend 
(when  Holt's  exponential  smoothing  model  is  used)  are 
stored  in  memory.  Although  the  forecast  for  each  period  is 
used  in  the  lot  size  procedure,  extended  forecasts  are  only 
developed  when  required  by  the  particular  lot  size 
heuristic  employed. 

To  provide  a  compact  computer  algorithm,  the  simple 
exponential  smoothing  procedure  was  incorporated  into 
Holt's  procedure  by  setting  the  trend  parameter,  g,  to 
zero. 

Focusing  is  carried  out  by  keeping  track  of  each 
individual  or  single  forecast's  mean  absolute  deviation  and 
smoothed  error  tracking  signal  (bias).  (Estimates  of  the 
MAD  are  also  exponentially  smoothed.)  The  forecast  with 
the  best  current  MAD  is  selected  for  the  focused  model  if 
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the  bias  is  within  acceptable  limits,  specifically  between 
-0.8  and  0.3. 

Silver  and  Peterson  (1985)  argue  that  a  negatively 
biased  forecast,  i.e.,  where  forecast  exceeds  demand,  is 
preferable  to  a  positively  biased  forecast,  i.e.,  where 
demand  exceeds  the  forecast,  since  being  a  few  items 
overstock  is  preferable  to  consistently  being  short 
(causing  too  many  premature  setups). 

Wemmerlov  and  Whybark  (1984)  specifically  avoid  the  use 
of  biased  forecasts  by  adjusting  the  average  actual  demand 
per  period  to  equal  the  average  forecast  demand  per  period. 
While  easily  done  for  simulated  demand  data,  it's  generally 
not  appropriate  for  empirical  demand  forecast  on  a  rolling 
horizon  basis. 

Research  by  Lee,  Adam,  and  Ebert  (1987)  show  that  "bias 
is  the  only  measure  that  satisfactorily  reflects  inventory 
carrying  cost...  [and]  only  bias  displays  any  reasonable 
association  with  the  shortage  cost  and  shortage  units...." 

Since  carrying  cost  is  caused  by  over  forecasting  (what 
they  refer  to  a  positive  bias)  and  shortage  costs  are 
caused  by  under  forecasting  (referred  to  as  negative  bias), 
the  use  of  an  unbiased  forecast  (as  used  by  Wemmerlov  and 
Whybark  (1984))  might  seem  reasonable.  The  research  by 
Lee,  et  al.  (1987),  however,  shows  that  "the  structures  of 
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these  two  component  costs  may  not  be  symmetrical  about  the 
zero  bias  level."  Unfortunately,  they  do  not  provide 
guidelines  as  to  what  the  nominal  bias  levels  may  be. 

The  specific  bias  levels  used  in  the  forecast  model 
were  determined  in  conjunction  with  an  outlier  discounting 
criterion.  An  example  data  set  which  exhibited  a  steep 
downward  trend  due  to  large  upward  spikes  (outliers)  was 
used.  The  steep  downward  trend  was  "leveled"  somewhat  by 
discounting  the  outliers  (more  on  this  in  a  moment)  and 
then  varying  the  bias  criteria  in  an  effort  to  eliminate  a 
large  series  of  zero  forecasts  caused  by  the  initial 
"trend".  (The  data  set  used  is  shown  in  Figure  1.) 

Outliers  were  discounted  by  keeping  track  of  the 
average  demand  and  standard  deviation  of  the  series  at  each 
point  in  the  forecast  "cycle".  If  an  outlier  exceeded  4 
standard  deviations,  the  actual  demand  was  reduced  to  the 
mean  plus  4  standard  deviations  for  forecast  purposes.  This 
provided  a  stabilizing  influence  on  the  forecast  which 
otherwise  would  have  to  have  been  provided  by  human 
intervention.  On  the  downside,  the  forecast  model  would 
lag  slightly  behind  a  true  shift  in  the  mean  of  the  demand 
series.  (This  type  of  lag  is  a  standard  "penalty"  for 
exponentially  smoothed  forecast  procedures.) 

The  Lot  Size  Procedure.  Other  than  the  Wagner-Whitten 
algorithm,  the  other  heuristics  are  simple  to  use  and  will 
not  be  discussed  here  (please  refer  to  Appendix  C  for  more 
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information).  Our  discussion  will  be  limited  to  that  part 
of  the  procedure  which  determines  our  production  policy. 

Research  on  lot  size  procedures  has  been  performed  by 
Silver  (1978),  Askin  (1981),  Bookbinder  and  Tan  (1983),  and 
Bookbinder  and  H'ng  (1986).  Our  procedure,  while  developed 
prior  to  our  knowledge  of  the  previous  works,  is  similar  to 
that  suggested  by  Bookbinder  and  Tan. 

Our  procedure  is  as  follows: 

1.  Use  a  focus  forecast  from  simple  exponential 
smoothing  and  exponential  smoothing  with  trend 
models  for  demands  over  the  rolling  horizon. 

2.  Treat  the  forecast  demands  as  deterministic  and 
employ  a  specific  lot  size  heuristic. 

3.  If  on-hand  inventory  is  positive,  the  amount 
produced  will  be  the  amount  obtained  from  the 
lot  size  heuristic  minus  the  on-hand  inventory. 

4.  If  on-hand  inventory  is  negative,  i.e.,  a  stockout 
has  occurred,  the  amount  produced  will  be  the 
amount  obtained  from  the  lot  size  heuristic  plus 
the  amount  backordered. 

5.  Each  period,  the  on-hand  inventory  is  compared  to 
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the  forecast  for  the  next  period.  If  our  forecast 
exceeds  our  inventory  position,  we  schedule  a 
setup  for  the  next  period,  otherwise  we  continue. 

6.  When  the  next  period's  demand  is  realized,  we 
either  meet  demand  or  we're  short.  If  a  shortage 
occurs,  a  setup  is  scheduled  for  the  next  period, 
otherwise  we  look  at  next  period's  forecast 
(Step  5) . 

7.  We  develop  an  extended  forecast  only  when  a  setup 
is  scheduled. 

8.  Continue  this  procedure  until  we  exhaust  all 
available  demand  data. 

9.  Discount  the  inventory  holding  cost  for  all  on- 
hand  inventory  used  to  satisfy  demand  beyond  the 
last  period  in  the  data  set. 

Figure  8  provides  a  flowchart  depicting  the  logic  of 
the  lot  size  procedure  employed. 

Program  Validation 

Verification  of  the  computer  model  was  obtained  through 
hand  calculations  and  analysis  of  the  results.  (Discussed 
in  a  separate  section.) 


Inv  -  Inv  -  0(t)  8top  If  Demand 
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Figure  8  —  Production  Procedure  (Flow  Chart) 
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The  optimal  Wagner-Whitten  procedure  and  all  lot  size 
heuristics  were  validated  by  hand  using  a  data  set  from  the 
first  group.  The  forecast  model  used  for  hand  verification 
of  the  heuristics  was  simple  exponential  smoothing  with  an 
alpha  parameter  of  0.2.  The  code  used  for  the 
Wagner-Whitten  procedures  was  an  equivalent  branch  and 
bound  algorithm  published  by  Jacobs  and  Khumawala  (1987). 
Verification  of  the  code  was  accomplished  by  comparing  the 
results  with  solutions  obtained  using  the  Wagner-Whitten 
algorithm  by  hand. 

The  forecast  procedure  was  also  validated  by  hand, 
however,  the  overall  focus  forecasting  policy  was  not. 
Instead,  we  validated  the  model  during  program  development 
by  comparing  the  focus  MAD  with  each  individual  MAD  for 
several  data  sets.  The  focus  forecast  compared  very 
favorably,  i.e.,  while  not  the  best,  it  was  significantly 
better  than  most. 

Verification  of  the  general  lot  size  procedure  was  also 
obtained  by  hand.  Specifically,  the  heuristics  were  run 
using  the  example  forecast  and  the  results  for  each 
transaction  printed  out  for  verification.  The  policy  was 
then  computed  by  hand  and  compared  with  the  printout. 

EXPERIMENTAL  DESIGN 

The  experimental  designs  for  each  group  of  data  was 
different  due  to  the  limited  number  of  data  sets  available 
in  the  second  group. 
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The  first  group  uses  an  unbalanced  5  factor  design  with 
5  performance  criteria  (3  of  which  are  the  costs  outlined 
previously;  the  other  2  are  measurements  involving  the  mean 
absolute  deviation  of  the  forecast  series).  The  design  is 
unbalanced  since  3  of  the  5  factors,  data  set,  data  type, 
and  degree  of  variation,  are  attributes  associated  with  the 
data  set.  (By  data  set,  we  mean  the  specific  data  set  of 
which  there  are  36.  Data  type  and  degree  of  variation  are 
as  defined  earlier  in  our  discussion  of  the  sample  data  and 
are  nested  within  data  set.)  The  other  two  factors  are,  of 
course,  the  lot  size  algorithm  and  TBO. 

The  TBO  factor  was  set  at  5  levels:  2,  4,  6,  8,  and  10. 
To  do  this,  we  set  the  TBO  level  a-priori  and  determined 
the  appropriate  A/h  ratio  based  upon  the  mean  or  average 
demand  of  each  data  set.  Our  procedure  is  therefore 
similar  to  the  studies  performed  by  Berry  (1972),  Callarman 
and  Hamrin  (1979),  Wemmerlov  and  Whybark  (1984),  and 
others . 

All  interactions  are  considered  except  the  5-way 
interaction  (as  it's  equivalent  to  the  error  term).  The 
second  group  was  handled  slightly  differently  in  that  only 
the  primary  factors,  lot  size  algorithm  and  TBO,  are  used 
in  the  ANOVA,  i.e.,  we  employ  a  simple  3-factor  balanced 


design  with  interaction. 
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RESULTS 

The  results  for  both  data  groups  are  very  similar, 
however  the  3  factor  design  for  the  second  group  wasn't 
able  to  discriminate  as  well  as  the  5  factor  unbalanced 
design  of  the  first.  The  ANOVA  results  are  presented  in 
Table  3. 


Variable 

Significance 

Group  1  | 

Group  2 

Cost 

0.001 

0.001 

Short 

0.001 

0.001 

%  Short 

0.001 

0.006 

Table  3.  Basic  ANOVA 


As  you  can  see,  both  ANOVAs  are  significant.  Tables  4 
and  5  provide  a  "breakdown"  of  the  significance  for  each 
factor  combination  for  Groups  1  and  2,  respectively.  The 
asterisk  denotes  significance  at  the  0.01  level. 

The  3-  and  4-way  interactions  are  generally  significant 
for  cost  although  not  for  shortages  or  amounts  short  in  the 
Group  1  ANOVA.  Results  are  similar  for  the  2-way 
interactions  in  Group  2.  The  results  for  all  common 
factors  and  their  interactions  are  generally  the  same  for 
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both  groups,  e.g.,  the  LOT(size)  factor  is  significant 
whereas  the  TYPE  factor  is  not.  (Note:  This  is  not 
necessarily  true  for  all  dependent  variables  or  "costs".) 


Factor 


Significance 


Cost 

Short 

%Short 

LOT 

0.0001* 

0.0001* 

0.0001* 

TBO 

0.0001* 

0.0001* 

0.0006* 

TYPE 

0.9926 

0.3008 

0.0001* 

VAR 

0.0001* 

0.0001* 

0.0001* 

SET ( VARxTYPE ) 

0.0001* 

0.0001* 

0.0001* 

LOTxTBO 

0.0001* 

0.0001* 

0.0001* 

LOTxTYPE 

0.5518 

0.4180 

0.6018 

LOTxVAR 

0.0030* 

0.0489 

0.1056 

LOTxSET ( VARxTYPE ) 

0.0001* 

0.0001* 

0.0009* 

TBOxTYPE 

0.4834 

0.3676 

0.3004 

TBOxVAR 

0.0001* 

0.0001* 

0.0037* 

TBOxSET ( VARxTYPE ) 

0.0001* 

0.0036* 

0.0001* 

VARxTYPE 

0.0126 

0.0010* 

1.0000 

LOTxTBOxTYPE 

0.0693 

0.8679 

0.5761 

LOTxTBOxVAR 

0.0001* 

0.2513 

0.5739 

LOTxTYPExVAR 

0.0001* 

0.2142 

0.7851 

TBOxTYPExVAR 

0.0001* 

0.0734 

0.0005* 

LOTxTBOxTYPExVAR 

0.0001* 

0.6336 

0.4499 

Table  4.  Detailed  ANOVA  (Group  1) 
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Factor 


Significance 


Cost 

Short 

%Short 

LOT 

0.0001* 

0.0001* 

0.0047* 

TBO 

0.0001* 

0.0001* 

0.0180 

SET 

0.0124 

0.3534 

0.0023* 

LOTxTBO 

0.0001* 

0.0004* 

0.4349 

LOTxSET 

0.0001* 

0.1247 

0.0146 

TBOxSET 

0.0011* 

0.2736 

0.0067* 

Table  5 

Detailed 

ANOVA  (Group  2) 

Tables  6  and  7  give  the  results  of  the  Tukey  multiple 
range  tests  for  each  single  factor  of  interest. 

Means  with  the  same  letters  are  not  significantly 
different  and  are  listed  from  high  to  low.  Note  that  the 


Factor 

j  Cost 

|  Short 

j  %Short 

Lotsize 

5  13  2  4 

A  B  CCCCC 

4  2  3  1  5 
AAA  C 

BBB  D 

4  2  3  1  5 
AAA  CCC  D 
BBB 

TBO 

X  8  2  6  4 

A  BBB  CCC 

X  8  6  4  2 

A  BBBBB  C 

X  4  2  6  8 
AAA 

BBBBBBB 

Type 

Lin  Non 

AAAAAAA 

Lin  Non 

AAAAAAA 

Lin  Non 

A  B 

Var 

3  2  1 

ABC 

3  2  1 

ABC 

3  2  1 

ABC 

Table  6.  Single  Factor  Tukey  Results  (Group  1) 
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Factor 

|  Cost 

|  Short 

|  %Short 

Lotsize 

5  14  3  2 

2  3  4  1  5 

4  2  3  1  5 

A  BBBBBBB 

AAAAAAA 

BBB 

AAAAAAAAA 

TBO 

X  8  2  6  4 

X  8  6  4  2 

X  8  6  2  4 

A  BBBBBBB 

AAAAAAA 

BBB 

AAAAA  CCC 
BBB 

Table  7.  Single  Factor  Tukey  Results  (Group  2) 


lotsize  heuristics  are  classified  as  before,  i.e.,  1  = 
Eisenhut,  2  =  EOQ,  3  =  Silver-Meal,  4  =  Tsado's  method,  and 
5  =  Wagner-Whitten  (non-optimal) .  The  X  denotes  a  TBO  of 
10,  "Lin"  is  short  for  linear,  and  non-linear  is 
abbreviated  as  "Non". 

There  are  8  2-factor  interactions  of  interest  which  are 
significant  in  the  Group  1  ANOVA:  LOTxTBO  and  TBOxVAR  for 
cost,  number  of  shortages,  and  amounts  short,  LOTxVAR  for 
cost,  and  VARxTYPE  for  shortages.  Figures  9  through  16 
depict  these  interactions.  Figures  17  and  18  provide  the 
results  of  2  significant  2-factor  interactions,  LOTxTBO  for 
both  cost  and  number  of  shortages  from  the  Group  2  ANOVA. 

All  interactions  involving  the  (data)SET  factor  are 
omitted  as  differences  in  cost  due  to  the  demand  series  is 
expected. 

The  reader  is  referred  to  Appendices  D  and  E  for 
additional  information. 
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Figure  9 
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2-Factor  Interaction  (Group  1) 
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Figure  10 
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Figure  11 
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Figure  17 
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Figure  18 
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ANALYSIS 

It's  obvious  from  Tables  6  and  7  that  different  lotsize 
algorithms  perform  differently  for  standard  inventory 
(holding  and  setup)  costs  than  for  shortage  related 
"costs" . 

In  this  study,  we've  shown  that  the  EOQ,  Silver-Meal, 
and  Tsado  algorithms  perform  significantly  better  for 
inventory  costs  than  do  the  Eisenhut  or  Wagner-Whitten 
heuristics.  Further,  the  inventory  costs  for  the  LOTxTBO 
interaction  depicted  in  Figure  9  are  relatively  stable  for 
these  “near-optimal"  performers. 

Tsado' s  algorithm  performed  best  for  cost  (and  worse 
for  shortages)  in  the  Group  1  ANOVA,  although  the 
differences  were  not  significant  between  the  3  best.  The 
Group  2  ANOVA  placed  Tsado 's  heuristic  in  third,  but  again, 
the  differences  among  the  best  performers  were  not 
significant . 

Conversely,  the  Eisenhut  and  Wagner-Whitten  heuristics 
perform  significantly  fewer  shortages  and  less  items  short 
per  stockout  (expressed  as  a  percentage  of  actual  demand 
for  the  stockout  period).  Figures  10  and  11  show  the 
relatively  wide  range  of  performance  for  the  various 
heuristics . 

Eisenhut  tends  to  be  the  most  stable  for  TBO  with 
shortages  occurring  between  12  and  18  percent  of  the 
replenishments  made.  Wagner-Whitten,  while  slightly  less 
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stable,  is  the  overall  best  performer.  It's  interesting  to 
note  that  performance  of  the  W-W  heuristic  as  a  "function 
of"  TBO  is  essentially  inverted.  Figure  11  shows  a  much 
greater  degree  of  interaction  between  all  five  heuristics 
for  amounts  short. 

Note  that  a  TBO  of  10  provides  the  W-W  solution  with 
virtually  no  stockouts  and  less  than  10  percent  of  actual 
demand  short  when  a  stockout  does  occur.  The  price, 
however,  is  an  average  cost  over  3  times  as  great  as  the 
true  optimal  solution.  In  all  cases,  the  Wagner-Whitten 
algoritnm  performed  worse  for  cost  and  best  for  shortages, 
although  significance  was  not  shown  for  amounts  short  in 
the  Group  2  ANOVA  (see  Table  7). 

It  therefore  appears  the  Wagner-Whitten  and  Eisenhut 
algorithms  maintain  a  significant  amount  of  inherent  safety 
stock  whereas  the  others  tend  to  "run  lean".  Additional 
inventory  wou.1 1  drive  up  the  inventory  holding  costs  while 
reducing  the  number  of  stockouts  due  to  being  a  few  items 
short . 

Table  6  shows  the  type  of  data  set,  i.e.,  whether  it's 
linear  or  non-linear,  does  not  affect  either  cost  or  the 
number  of  shortages  significantly.  Non-linearity  does  seem 
to  affect  the  amount  short  per  stockout.  It  appears  that 
non-linear  demand  has  a  smaller  percentage  short  per 
stockout  due  to  "over-forecasting"  by  our  linear  forecast 
model.  The  degree  of  variance,  however,  is  significant 
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for  all  dependent  variables  with  lower  overall  costs 
associated  with  lower  variability  (as  would  be  expected) . 

Although  the  significance  level  was  set  at  0.01,  Figure 
13  (TBOxVAR  -  Cost)  shows  interaction  between  TBO  levels  of 
2  and  4.  Although  lower  variance  generally  implies  lower 
overall  costs  (i.e.,  both  inventory  and  shortage 
related),  the  reverse  is  true  for  the  TBO  factor  at  a  level 
of  2.  We  would  like  to  point  out  that  actual  cost 
performance  at  this  level  may  not  be  that  "significant" 
(41.3,  41.7,  and  39.0  for  low,  medium,  and  high  variances, 
respectively) . 

Cost  performance  as  a  function  of  TBO  tends  to  validate 
the  "rule  of  thumb"  advocated  by  Wemmerlov  and  Whybark 
(1984)  which  states  the  length  of  the  rolling  horizon 
should  be  approximately  3  times  the  average  time  between 
orders.  A  TBO  of  4  was  best  using  our  rolling  horizon  of 
12  periods.  Performance  for  number  of  stockouts  and 
amounts  short  per  stockout,  however,  do  not  support  this 
assertion. 

We  should  note  at  this  point  that  Wemmerlov  and  Whybark 
(1984)  also  advocate  a  rolling  horizon  of  at  least  5  times 
the  average  time  between  orders  for  the  Wagner-Whitten 
algorithm  when  used  as  a  heuristic.  The  argument  in  this 
case  was  the  fact  the  W-W  heuristic  utilized  the  entire 
length  of  the  rolling  horizon  in  order  to  make  its  initial 
production  decision.  From  Appendix  D,  we  see  that  the  mean 
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cost  performance  for  the  W-W  heuristic  for  a  TBO  of  2  is 
41.015926  and  for  a  TBO  of  4  is  32.546886.  This  tends  to 
discredit  their  assertion. 

There  may  be  2  reasons  for  this.  First,  previous  cost 

criteria  "consisted"  of  standard  inventory  costs  and  an 

% 

"artificial"  service  level.  The  mean  "cost"  performance 
for  stockouts  is  better  for  the  Wagner-Whitten  heuristic 
with  a  TBO  of  2  than  a  TBO  of  4  (9.1463710  and  14.3588956, 
respectively),  therefore  the  "total  costs"  due  to  inventory 
and  stockouts  may  "average"  out.  (Note  there  was  no  real 
difference  in  the  amounts  short  for  either  case.)  Second, 
forecasts  generally  tend  to  "worsen"  as  they  extend  further 
into  the  future.  This  would  tend  to  imply  that  using 
increasingly  "bad"  data  in  order  to  make  the  initial 
production  decision  results  in  a  "worse"  decision. 

The  cost  results  for  the  TBO  factor  (where  a  TBO  which 
is  one-third  of  our  rolling  horizon  is  "optimal")  tend  to 
validate  the  production  procedure  and  the  computer  code  in 
general. 

Other  results  which  tend  to  validate  this  research  are 
as  follows: 

1.  The  difference  between  all  heuristic  costs  and 
the  optimal  cost  are  strictly  positive. 

2.  The  average  difference  between  the  mean  absolute 
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deviation  of  the  lotsize  forecasts  and  the 
focused  forecast  for  each  data  set  is  strictly 
positive. 

3.  The  amount  of  lotsize  forecast  error  decreases 
with  TBO. 

4.  Lotsize  forecast  error  is  greater  for  non-linear 
than  for  linear  data  sets. 

5.  Lotsize  forecast  error  is  generally  greater  for 
data  sets  with  a  higher  degree  of  variance. 

And  finally,  we  would  like  to  note  that  the  Group  2 
design  did  not  discriminate  as  well  as  the  Group  1  design. 
This  is  probably  due  to  2  factors.  First,  the  number  of 
datasets,  and  therefore  observations,  is  smaller  in  Group  2 
than  Group  1.  Second,  the  design  for  Group  1  is  more 
"complex",  i.e.,  it  accounts  for  more  error  than  the 
simpler  Group  2  model.  It's  interesting  to  note,  however, 
that  Figures  17  and  18  show  essentially  the  same  thing  as 
Figures  9  and  10,  i.e.,  both  data  groups  (Groups  I  and  II) 

show  ai roximately  the  same  average  performance,  both 
absolute  (numerically)  and  relative  (to  each  other),  for 
all  lotsize  algorithms  at  all  factor  levels. 


CHAPTER  VI 


CONCLUDING  REMARKS 

The  purpose  of  this  research  was  three-fold.  First, 
there  existed  a  need  to  perform  a  study  of  lotsize 
heuristic  performance  which  forecasts  empirical  demand  data 
in  the  same  manner  in  which  the  lotsize  heuristics  are 
implemented,  i.e.,  over  a  rolling  horizon.  Second,  we 
wished  to  analyze  shortage  costs  separately  from 
traditional  inventory  costs.  And  finally,  we  wished  to 
validate  the  heuristic  presented  by  Tsado  (1985). 

The  lack  of  significance  for  the  LOTxVAR  interaction 
validates  the  assertion  of  previous  studies  that,  although 
higher  variance  increases  costs,  it  doesn't  alter  the 
relative  performance  of  the  lot  size  heuristics. 

We  have  shown,  however,  that  results  from  previous 
studies  were  confounded  due  to  the  way  shortage  costs  were 
handled,  i.e.,  by  setting  an  arbitrary  service  level. 
Whereas  previous  studies  showed  no  significant  difference 
in  lot  size  heuristic  performance,  we  have  shown  that  some 
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lot  size  heuristics  perform  better  for  traditional 
inventory  costs  and  others  perform  better  with  regard  to 
the  average  number  of  stockouts  expressed  as  a  percentage 
of  the  total  number  of  replenishments  made  and  average 
amounts  short  as  a  percentage  of  actual  demand  for  the 
stockout  period. 

Whether  this  result  has  any  application  to  the  way 
lotsize  heuristics  are  chosen  and  implemented  or  not  is  for 
industry  to  decide.  (Specifically,  a  company  will  probably 
set  a  service  level  based  upon  its  own  perception 
of  the  "cost"  associated  with  a  shortage.)  However,  we 
believe  this  result  is  important  in  that  it  shows  that 
researchers  should  carefully  choose  their  basic 
assumptions.  In  this  case,  previous  researchers  would  have 
obtained  results  similar  to  this  research  if  they  had  set 
service  level  as  a  factor  in  their  experimental  designs. 
But  by  setting  a  single  service  level,  they  "confounded" 
the  relative  performance  of  the  various  lot  size  heuristics 
examined. 

And  finally,  we've  shown  that  Tsado's  algorithm 
performed  very  well  for  traditional  inventory  holding  and 
setup  costs  and  have  therefore  validated  his  1985  study. 

This  study  was  by  no  means  all-encompassing.  The 
following  are  suggestions  for  further  research: 
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1.  A  similar  study  is  required  using  a  balanced 
experimental  design  similar  to  that  used  by 
Group  1  and  using  a  more  complete  array  of 
lotsize  techniques. 

2.  Separate  research  regarding  the  relative  merit 
of  focus  or  "adaptive"  and  averaged  exponential 
smoothing  techniques  (as  used  in  automatic 
forecasting  over  a  rolling  horizon)  is  warranted. 
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GROUP  2  DATA 


PERIOD/SET 

1 

2 

3 

4 

5 

1 

1053 

441 

548 

783 

176 

2 

1254 

558 

795 

511 

92 

3 

1566 

402 

993 

2116 

44 

4 

1659 

234 

1103 

2998 

84 

5 

1143 

363 

927 

954 

124 

6 

1374 

306 

831 

1299 

254 

7 

1353 

189 

845 

1245 

63 

8 

1296 

114 

665 

1274 

141 

9 

1920 

387 

888 

1151 

147 

10 

1248 

444 

539 

1316 

86 

11 

1104 

327 

363 

804 

175 

12 

909 

318 

370 

488 

34 

13 

945 

239 

425 

832 

131 

14 

870 

376 

375 

2076 

47 

15 

1247 

399 

451 

835 

154 

16 

990 

398 

348 

1124 

68 

17 

1422 

281 

461 

1436 

68 

18 

1359 

282 

459 

1267 

137 

19 

1113 

227 

437 

2020 

150 

20 

1479 

340 

465 

832 

73 

21 

1455 

380 

401 

909 

113 

22 

1029 

671 

426 

718 

65 

23 

998 

562 

102 

646 

60 

24 

966 

522 

0 

656 

61 

25 

627 

277 

1136 

633 

24 

26 

723 

341 

821 

744 

64 

107 


108 


PERIOD/SET 

1 

2 

3 

4 

5 

27  ' 

885 

512 

848 

632 

36 

28 

1323 

431 

1231 

676 

27 

29 

1137 

361 

928 

685 

106 

30 

1389 

386 

867 

819 

79 

31 

1455 

3880 

656 

402 

72 

32 

1110 

555 

1084 

519 

38 

33 

1470 

1583 

1302 

694 

61 

34 

867 

2667 

1263 

531 

46 

35 

1095 

2626 

997 

486 

68 

36 

861 

2584 

1208 

517 

26 

37 

465 

1907 

1353 

247 

20 

38 

771 

1540 

1038 

404 

86 

39 

756 

1163 

1207 

330 

29 

40 

1038 

908 

1223 

245 

22 

41 

1356 

863 

1002 

213 

43 

42 

1521 

751 

1173 

209 

34 

43 

1194 

683 

802 

103 

90 

44 

1158 

2231 

1548 

190 

15 

45 

1569 

2578 

2309 

309 

33 

46 

957 

922 

882 

133 

28 

47 

798 

1085 

1094 

99 

15 

48 

686 

1059 

1305 

168 

15 

49 

633 

1162 

1918 

266 

214 

50 

639 

973 

1216 

180 

120 

51 

894 

749 

2409 

274 

108 

52 

837 

971 

1907 

331 

106 

109 


PERIOD/SET 

1 

2 

3 

4 

5 

53 

918 

903 

1370 

174 

136 

54 

987 

974 

1136 

120 

96 

55 

1035 

979 

936 

78 

31 

56 

840 

858 

1284 

98 

72 

57 

1107 

1105 

2340 

96 

156 

58 

801 

1419 

1444 

121 

106 

59 

714 

1119 

1530 

146 

102 

60 

696 

963 

1376 

171 

148 

61 

570 

919 

1319 

109 

96 

62 

357 

937 

1480 

125 

100 

63 

525 

1588 

1291 

201 

108 

64 

696 

876 

1407 

127 

120 

65 

870 

656 

1108 

110 

113 

66 

570 

508 

978 

139 

139 

67 

657 

581 

1149 

55 

151 

68 

750 

867 

776 

132 

92 

69 

669 

1198 

1566 

84 

96 

70 

474 

1368 

1216 

180 

102 

71 

402 

1104 

1389 

231 

148 

72 

414 

1190 

1719 

135 

194 

73 

315 

894 

882 

127 

145 

74 

402 

883 

1266 

177 

153 

75 

504 

887 

1499 

116 

147 

76 

429 

906 

1141 

175 

144 

77 

372 

833 

1153 

113 

157 

78 

315 

697 

996 

158 

197 

APPENDIX  C 


PROGRAM  LISTING 


DECLARE  SOB  eoq  ( ) 

DECLARE  SUB  tsado  { ) 

DECLARE  SUB  wwheuristic  ( ) 

DECLARE  SUB  eoqtimesupply  ( ) 

DECLARE  SUB  eisenhut  ( ) 

DECLARE  SUB  silver meal  ( ) 

DECLARE  SUB  wagnerwhitten  ( ) 

DECLARE  SUB  forecaster  { ) 

COMMON  SHARED  /Adata/  cost,  count,  datatype,  holding, 

horizon,  inventory,  k9 

COMMON  SHARED  /Bdata/  lastn,  mad,  n,  production,  setup, 

sigma,  xbar 

DEFSNG  A-Z 

x  =  10  'Dummy  variable  to  allow  dynamic  dimensioning  of 
demand(n)  and  forecast (n) 

DIM  SHARED  D(x),  demand(x),  forecast(x,  x),  lotdata(x) 

DIM  SHARED  R(x),  a(x,  x),  0(x) 

'Iterate  through  all  'DATA  sets 

FOR  dataset  =  1  TO  5 

READ  xbar,  sigma,  datatype,  n 

horizon  =12  '12  weeks  (3  months)  for  weekly  DATA  or  12 

months  for  yearly  data 
ERASE  demand,  forecast,  D,  lotdata 

REDIM  SHARED  D(horizon),  demand(n),  forecast(n  +  1,  2), 
lotdata(n  +  horizon) 

FOR  i  =  1  TO  n:  READ  demand(i):  NEXT  i 

IF  n  =  52  THEN  lastn  =  6  ELSE  lastn  =  12 

'Compute  variation  as  Low  (1),  Medium  (2),  or  High  (3) 

IF  sigma  /  xbar  <=  .5  THEN 
variation  =  1 

ELSE 

IF  sigma  /  xbar  >  1!  THEN 
variation  =  3 

ELSE 


variation  =  2 
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END  IF 
END  IF 

'Generate  optimal  forecast 
CALL  forecaster 

'Begin  iterations  through  all  values  of  TBO 
FOR  TBO  =  2  TO  10  STEP  2 

setup  =  INT( . 5  +  .5  *  xbar  *  TBO  ~  2):  holding  =  1 

IF  n  =  52  THEN  lastn  =  6  ELSE  lastn  =  12 

start  =  lastn 

'Get  optimal  solution 

CALL  wagnerwhitten 

'Begin  iterations  through  all  lotsize  algorithms 
FOR  lotsize  =  1  TO  5 

IF  n  =  52  THEN  lastn  =  6  ELSE  lastn  =  12 

'Start  basic  algorithm 

lastn  =  lastn  +  1:  short  =  0:  count  =  0 

58  IF  lastn  >  n  THEN  2101 

FOR  i  =  1  TO  horizon 

D ( i )  =  forecast ( lastn,  1)  +  (i  -  1)  *  forecast (lastn,  2) 

IF  D(i)  <  0  THEN  D(i)  =  0 
lotdata( lastn  +  i  -  1)  =  D(i) 

NEXT  i 

D( 1 )  =  D( 1 )  +  short 

'Select  appropriate  lotsize  technique 

ON  lotsize  GOTO  665,  666,  667,  668,  669 

665  CALL  eisenhut 
GOTO  701 

666  CALL  eoq 
GOTO  701 

667  CALL  silvermeal 
GOTO  701 

668  CALL  tsado 
GOTO  701 

669  CALL  wwheuristic 

'Compute  inventory  level  and  cost  after  production 

701  inventory  =  inventory  +  production 

cost  =  cost  +  setup:  count  =  count  +  1 

'Compute  inventory  and  holding  costs  after  demand  for 

'current  period  is  satisfied 

61  IF  lastn  >  n  GOTO  2101 

inventory  =  inventory  -  demand( lastn)  -  short 
cost  =  cost  +  inventory  *  holding 
IF  lastn  >=  n  GOTO  2101 
lastn  =  lastn  +  1 

'Determine  if  forecast  demand  exceeds  inventory 

IF  inventory  >  lotdata(lastn)  THEN  'Inventory  exceeds 

forecast 

IF  demand ( lastn)  >  inventory  THEN  'Demand  too  large 

—  shortage 

short  =  demand (lastn)  -  inventory 
inventory  =  0 

percentshort  =  percentshort  +  short  /  demand( lastn) 
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shortcounter  =  shortcounter  +  1 
lastn  =  lastn  +  1 

GOTO  58  'Setup  a  new  production  run 
ELSE  'Demand  was  less  than  our  inventory  —  deliver 
current  demand 
short  =  0 

GOTO  61  'Deliver  next  period's  demand 
END  IF 

ELSE  'Next  period's  forecast  for  demand  exceeds  current 
inventory 

short  =  0 

GOTO  58  'Setup  a  new  production  run 
END  IF 

'Calculate  the  true  mean  absolute  deviation  for  all  periods 
'included  in  the  lotsize  problem 
2101  FOR  1  =  start  +  1  TO  n 

lotsizeforecasterror  =  lotsizeforecasterror  + 

ABS(lotdata(l)  -  demand{l)) 

lotsizetrackingsignal  =  lotsizetrackingsignal  +  demand(l)  - 
lotdata( 1) 

NEXT  1 

lotsizeforecasterror  =  lotsizeforecasterror  /  (n  -  start) 
lotsizetrackingsignal  =  lotsizetrackingsignal  / 

(lotsizeforecasterror  *  (n  -  start)) 
'Subtract  holding  costs  for  periods  beyond  N  from  the  total 
'cost  —  this  will  reduce  the  variability  between  TBOs 
' ( treatments) . 

cost  =  cost  -  inventory  *  holding 

'Compute  PERCENTSHORT  if  SHORTCOUNTER  is  nonzero  (to 
'prevent  division  by  0). 

IF  shortcounter  >  0  THEN  percentshort  =  percentshort  / 
shortcounter  *  100 

'Change  SHORTCOUNTER  into  the  fraction  of  times  short  to 
'number  of  replenishments  made  (COUNT  =  (approx.)  TBO  — 
'this  will  reduce  the  variability  due  to  TBO  (e.g.  a  model 
'with  a  TBO  of  2  can  havequite  a  few  more  shortages  in  a 
'given  time  period  than  a  model  with  a  TBO  of  10). 
shortcounter  =  100  *  shortcounter  /  count 
'Output  results  to  'DATA  file  in  current  directory. 

OPEN  "THESIS2 .OUT"  FOR  APPEND  AS  #1 

WRITE  #1,  dataset,  datatype,  variation,  sigma  /  xbar,  TBO, 
lotsize,  100  *  (cost  -  k9)  /  k9,  shortcounter, 
percentshort,  100  *  (lotsizeforecasterror  -  mad) 

/  mad,  100  *  lotsizeforecasterror  /  xbar,  100  * 
lotsizetrackingsignal 

CLOSE  #1 

'Note:  The  'DATA  is  given  in  percentages  to  provide 
'non-exponential  format  in  the  output  (in  order  to  allow 
'SAS  to  read  the  DATA  directly  (after  the  commas  are 
' removed) ) . 

'Zero  appropriate  variables. 

inventory  =  0:  production  =  0:  lotsizeforecasterror  =  0: 
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lotsizetrackingsignal  =  0 

shortcounter  =  0:  percentshort  =  0:  cost  =  0 

ERASE  lotdata 

REDIM  lotdata(n  +  horizon) 

NEXT  lotsize 
NEXT  TBO 
NEXT  dataset 
END 

'DATA  statements 
SUB  forecaster 

'This  subroutine  computes  a  focused  forecast  using  Holt's 
'exponential  smoothing  model  and  simple  exponential 
'smoothing  (special  case  of  Holt's  model  where  gamma  =  0). 
DIM  holt(4,  5,  6),  bestholt(6),  stat(n) 

'Copy  the  demand  'DATA  to  a  new  matrix  for  the  algorithm 

'which  smoothes  out  outliers 

FOR  i  =  1  TO  n:  stat(i)  =  demand(i):  NEXT  i 

'Store  the  various  levels  of  Alpha  and  Gamma  in  the 

'forecast  matrix  called  HOLT 


FOR  j  =  1  TO  4:  FOR 
NEXT  k:  NEXT  j 

k  =  1  TO  5: 

holt( j. 

k. 

1)  =  j  * 

.1: 

FOR  k  =  1  TO  5:  FOR 
.Is  NEXT  j:  NEXT  k 

j  =  1  TO  4: 

holt ( j , 

k , 

2)  =  (1  - 

k)  * 

w  =  .05  'Smoothing  constant  for  the  exponential  smoothing 
form  of  MAD  averagedemand  =  (demand(l)  + 
demand(2))  /  2  'Initialize  averge  demand 
'Compute  the  first  forecast  for  all  Alpha  and  Gamma  levels 
'as  the  average  of  the  first  two  periods  of  demand. 

FOR  j  =  1  TO  4 
FOR  k  =  1  TO  5 

holt(j,  k,  3)  =  averagedemand 

holt(j,  k,  5)  =  ABS(holt(j,  k,  3)  -  demand(3))  'Initialize 
individual  MADs 
NEXT  k 
NEXT  j 

'Begin  forecast  procedure 

FOR  h  =  3  TO  n  'Iterate  through  N  periods 
'Compute  current  average  and  std  deviation  of  the  demand 
'sequence  for  H-l  periods 
averagedemand  =  0:  sum  =  0 

FOR  ij  =  1  TO  h  -  1:  sum  =  sum  +  stat(ij):  NEXT  ij 
averagedemand  =  sum  /  (h  -  1) 
sumsquare  =  0 

FOR  ik  =  1  TO  h  -  1:  sumsquare  =  sumsquare  +  (stat(ik)  - 

averagedemand)  ~  2: 

NEXT  ik 

deviation  =  (sumsquare  /  (h  -  2))  ~  .5 

bestholt(5)  =  9999999999#  'Set  high  initial  value  of  MAD 
FOR  j  =  1  TO  4  'Iterate  through  all  4  Alpha  levels 

FOR  k  =  1  TO  5  'Iterate  through  all  5  Gamma  levels 

residual  =  demand(h)  -  holt(j,  k,  3)  'Store  difference 


between  last  estimate 
and  actual  demand 

'Discount  outliers  (virtually  all  points  should  fall  within 
' 3  or  4  std  deviations  if  the  points  are  distributed 
’normally) . 

IF  datatype  =  1  THEN  'Only  valid  for  constant  and 
slightly  trending  demand 
IF  deviation  >  0  THEN 

IF  ABS( residual )  /  deviation  >  4  THEN 

residuaL  =  SGN{ residual)  *  4  *  deviation 
stat(h)  =  residual  +  averagedemand 

ELSE 
END  IF 

ELSE 
END  IF 

ELSE 
END  IF 

temp  =  holt(j,  k,  3)  'Store  last  estimate  (the  one  for 
current  period's 
demand ) 

'Compute  new  estimate. 

holt ( j ,  k,  3)  =  (1  -  holt ( j ,  k,  1))  *  (holt(j,  k,  3)  + 

holt(jr  k,  4))  +  holt(j,  k ,  1)  *  stat(h) 

IF  holt ( j ,  k,  3)  <  0  THEN  holt ( j ,  k,  3)  =  0 
'Compute  new  trend. 

holt(j,  k,  4)  =  (1  -  holt(j,  k,  2))  *  holt(j,  k,  4)  + 

holt(j,  k,  2)  *  (holt(j,  k,  3)  -  temp) 
'Compute  current  MAD  and  current  smoothed  forecast  error, 
holt ( j ,  k,  5)  =  (1  -  w)  *  (holt ( j ,  k,  5))  +  w  * 

ABS( residual) 

holt(j,  k,  6)  =  (1  -  w)  *  (holt(j,  kf  6))  +  w  *  residual 
'Determine  the  best  forecast  to  date, 
trackingsignal  =  holt(j,  k,  6)  /  holt(j,  k,  5) 

'As  defined  by  Silver  and  Peterson  (1985) 

'Note:  A  negatively  biased  forecast  (i.e.,,  where  forecast 
'exceeds  demand)  is  preferable  to  a  positively  biased 
'forecast  (i.e.,,  where  demand  exceeds  forecast)  since 
being 

'a  few  items  overstock  ("safety  stock")  is  preferable  to 
'consistently  being  a  few  items  short  (causing  too  many 
'premature  setups).  [See  EOQ  example.  Ibid.] 

IF  trackingsignal  <  .3  AND  trackingsignal  >  -.9  THEN 
IF  holt ( j ,  k,  5)  <  bestholt( 5)  THEN 
bestholt(l)  =  holt ( j ,  k,  1) 

'Alpha  of  best  forecast  for  the  period 
bestholt(2)  =  holt(j,  k,  2) 

'Gamma  of  best  forecast  for  the  period 
bestholt(3)  =  holt(j,  k,  3) 

'Best  smoothed  estimate 
bestholt(4)  =  holt(j,  k,  4) 

'Best  smoothed  trend 
bestholt(5)  =  holt(j,  k,  5) 
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'Best  MAD  for  the  period 
bestholt{6)  =  holt(j,  k,  6) 

'Best  smoothed  forecast  error  for  the  period 

ELSE 
END  IF 

ELSE 
END  IF 

NEXT  k  'Next  Gamma 
NEXT  j  'Next  Alpha 

'Now  store  the  best  forecast  for  this  period  (which  is  the 
'forecast  for  'demand  in  period  H+l,  i.e.,  the  next 
period) . 

forecast(h  +  1,  1)  =  INT(.5  +  bestholt(3)) 
forecast (h  +  1,  2)  =  INT(.5  +  bestholt(4)) 

NEXT  h  'Next  period 

'Compute  the  MAD  for  the  entire  forecast  over  N  -  Lastn 

'periods. 

mad  =  0 

FOR  i  =  lastn  +  1  TO  n:  mad  =  mad  +  ABS( forecast ( i,  1)  - 

demand ( i ) ) :  NEXT  i 

mad  =  mad  /  (n  -  lastn) 

'Compute  the  tracking  signal  for  the  entire  forecast  over  N 
'-  Lastn  periods, 
track  =  0 

FOR  i  =  lastn  +  1  TO  n:  track  =  track  +  demand(i)  - 
forecast(i,  1):  NEXT  i 

trackingsignal  =  track  /  (mad  *  (n  -  lastn)) 

END  SUB 

SUB  silvermeal 

DIM  trc( 24) ,  trcut ( 24 ) 

trc(l)  =  0:  trcut(l)  =  0:  production  =  0 
trc(l)  =  setup:  trcut(l)  =  trc(l) 

FOR  kk  =  2  TO  horizon 

trc(kk)  =  trc(kk  -  1)  +  (kk  -  1)  *  D(kk)  *  holding 

trcut(kk)  =  trc(kk)  /  kk 

IF  trcut ( kk)  >  trcut(kk  -  1)  THEN  55 

NEXT  kk 

55  FOR  1  =  1  TO  kk  -  1 
production  =  production  +  D(l) 

NEXT  1 

production  =  production  -  inventory 
END  SUB 

SUB  eisenhut  'part-period  balancing 
DIM  trc( 24) ,  trcut ( 24 ) 
trc(l)  =  0:  production  =  0 
FOR  kk  =  2  TO  horizon 

trc(kk)  =  trc(kk  -  1)  +  (kk  -  1)  *  D(kk)  *  holding 
'determine  first  period  where  accumulated  holding  costs 
'exceeds  setup 
IF  trc(kk)  >  setup  THEN 
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'determine  which  integer  period  is  closer  to  actual  setup 
'costs 

IF  ABS(trc(kk)  -  setup)  >  ABS(trc(kk  -  1)  -  setup) 
THEN 

lastperiod  =  kk  -  1 

ELSE 

lastperiod  =  kk 
END  IF 
GOTO  565 

ELSE  'continue  to  next  period,  i.e.,  next  kk 
END  IF 
NEXT  kk 

565  FOR  1  =  1  TO  lastperiod 
production  =  production  +  D(l) 

NEXT  1 

production  =  production  -  inventory 
END  SUB 

SUB  eoq 

production  =  0:  avg  =  0 
FOR  i  =  1  TO  horizon 
avg  =  avg  +  D(i) 

NEXT  i 

avg  =  avg  /  horizon 

production  =  (2  *  avg  *  setup  /  holding)  "  .5  +  short 
production  =  INT(.5  +  production)  -  inventory 
END  SUB 

SUB  tsado 

production  =  0:  avg  =  0 
FOR  i  =  1  TO  horizon 
avg  =  avg  +  D(i) 

NEXT  i 

avg  =  avg  +  short:  avg  =  avg  /  horizon 
time  =  lastn  -  7 

production  =  -avg  *  time  +  {(avg  *  time)  A  2  +  (2  *  avg  * 
(cost  +  setup)  / 
holding))  *  .5 

production  =  INT(.5  +  production)  -  inventory 
END  SUB 

SUB  wagnerwhitten 
5  ERASE  R,  a,  O 

10  REDIM  R(n) ,  a(5000,  5),  0{n) 

60  m  =  n  -  lastn 

110  FOR  i  =  1  TO  m:  R(i)  =  demand (i  +6):  NEXT  i 

120  S  =  setup 

130  C  =  holding 

150  a(l,  1)  =  1 

160  a(l,  2)  =  1 

170  a(l,  3)  =  S 

180  a(l,  4)  =  0 
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190  a(l,  5)  =  1 

200  N1  =  0 

210  N2  =  1 

220  FOR  i  =  1  TO  m 

230  IF  R ( i )  =0  THEN  530 

240  k9  =  999999999# 

250  K8  =  0 

260  FOR  k  =  1  TO  N2 

270  IF  a(Nl  +  k,  3)  >  k9  THEN  320 

280  K8  =  N1  +  k 

290  k9  =  a(K8,  3) 

300  K7  =  k9  +  S 
320  NEXT  k 
330  N3  =  1 
340  j  =  N1  +  N2 
350  a(j  +  N3 ,  1)  =  i 

360  a(j  +  N3 ,  2)  =  1 

370  a(j  +  N3 ,  3)  =  K7 

380  a(j  +  N3 ,  4)  =  K8 

390  a(j  +  N3 ,  5)  =  i 

400  FOR  k  =  1  TO  N2 

410  Cl  =  (i  -  a(Nl  +  k,  5))  *  C  *  R(i) 

420  IF  Cl  >  S  THEN  500 

430  IF  a(Nl  +  k,  4)  +  Cl  >  K7  THEN  540 

440  N3  =  N3  +  1 

450  a(j  +  N3 ,  1)  =  i 

460  a(j  +  N3,  2)  =  0 

470  a(j  +  N3 ,  3)  =  a(Nl  +  k,  3)  +  Cl 

480  a(j  +  N3 ,  4)  =  N1  +  k 

490  a(j  +  N3 ,  5)  =  a(Nl  +  k,  5) 

500  NEXT  k 

510  N1  =  j 

520  N2  =  N3 

530  NEXT  i 

540  k9  =  999999999# 

550  K8  =  0 

560  FOR  k  =  1  TO  N2 

570  IF  a(Nl  +  k,  3)  >  k9  THEN  600 

580  K8  =  N1  +  k 

590  k9  =  a(K8,  3) 

600  NEXT  k 

610  'Solution  Cost  =  k9 
620  IF  a(K8,  2)  =  0  THEN  640 
630  0(a(K8,  1) )  =  1 
640  K8  =  a(K8,  4) 

650  IF  a{K8,  4)  =  0  THEN  670 

660  GOTO  620 

670  FOR  i  =  1  TO  m 

680  IF  i  =  1  THEN  720 

690  IF  O(i)  =  0  THEN  730 

720  Q  =  0 

730  Q  =  Q  +  R{ i ) 
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740  NEXT  i 
END  SUB 

SUB  wwheuristic 
1105  ERASE  R,  a,  0 

1110  REDIM  R(horizon),  a(1000,  5),  O(horizon) 
1115  production  =  0:  flag  =  0 
1160  m  =  horizon 

11110  FOR  i  =  1  TO  m:  R{i)  =  D(i):  NEXT  i 

11120  S  =  setup 

11130  C  =  holding 

11150  a(l,  1)  =  1 

11160  a(l,  2)  =  1 

11170  a(l,  3)  =  S 

11180  a(I,  4)  =  0 

11190  a(l,  5)  =  1 

11200  N1  =  0 

11210  N2  =  1 

11220  FOR  i  =  1  TO  m 

11230  IF  R(i)  =  0  THEN  11530 

11240  k999  =  999999999# 

11250  K8  =  0 

11260  FOR  k  =  1  TO  N2 

11270  IF  a(Nl  +  k,  3)  >  k999  THEN  11320 
11280  K8  =  N1  +  k 
11290  k999  =  a(K8,  3) 

11300  K7  =  k999  +  S 
11320  NEXT  k 
11330  N3  =  1 
11340  j  =  N1  +  N2 
11350  a(j  +  N3 ,  1)  =  i 

11360  a(j  +  N3 ,  2)  =  1 

11370  a(j  +  N3 ,  3)  =  K7 

11380  a(j  +  N3 ,  4)  =  K8 

11390  a(j  +  N3 ,  5)  =  i 

11400  FOR  k  =  1  TO  N2 

11410  Cl  =  (i  -  a(Nl  +  k,  5))  *  C  *  R(i) 

11420  IF  Cl  >  S  THEN  11500 

11430  IF  a(Nl  +  kf  4)  +  Cl  >  K7  THEN  11540 

11440  N3  =  N3  +  1 

11450  a ( j  +  N3 ,  1)  =  i 

11460  a( j  +  N3 ,  2)  =  0 

11470  a( j  +  N3 ,  3)  =  a(Nl  +  k,  3)  +  Cl 

11480  a( j  +  N3 ,  4)  =  N1  +  k 

11490  a( j  +  N3,  5)  =  a(Nl  +  k,  5) 

11500  NEXT  k 

11510  N1  =  j 

11520  N2  =  N3 

11530  NEXT  i 

11540  k999  =  999999999# 

11550  K8  =  0 

11560  FOR  k  =  1  TO  N2 
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11570  IF  a (N1  +  k,  3)  >  k999  THEN  11600 
11580  K8  =  N1  +  k 
11590  k999  =  a(K8,  3) 

11600  NEXT  k 

11620  IF  a(K8,  2)  =  0  THEN  11640 
11630  0{a(K8,  1) )  =  1 
11640  K8  =  a(K8,  4) 

11650  IF  a{K8,  4)  =  0  THEN  11670 
11660  GOTO  11620 
11670  FOR  i  =  1  TO  in 
11680  IF  i  =  1  THEN  11720 
11690  IF  O(i)  =  0  THEN  11730 

11705  IF  flag  =  1  THEN  11720 

11706  production  =  Q 

11707  flag  =  1 
11720  Q  =  0 
11730  Q  =  Q  +  R(i) 

11740  NEXT  i 

11800  production  =  production  -  inventory 
END  SUB 


APPENDIX  D 


ANOVA  RESULTS  —  GROUP  I  DATA 


The  following  text  provides  the  code  used  in  the  SAS 
routine.  Output  consists  of  all  subsequent  pages. 


DATA; 

INPUT  SET  TYPE  VAR  COFVAR  TBO  LOT  COST  SHORT 

PERSHORT  DELERR  PERERR  BIAS; 

DROP  COFVAR  BIAS; 

CARDS; 

PROC  ANOVA; 

CLASS  SET  TBO  LOT  VAR  TYPE; 

MODEL  COST  SHORT  PERSHORT  DELERR  PERERR  =  LOT  TBO 
TYPE  VAR  SET (TYPE  VAR)  LOT*TBO  LOT*TYPE 
LOT* VAR  LOT* SET (TYPE  VAR)  TBO*TYPE  TBO*VAR 
TBO* SET (TYPE  VAR)  TYPE* VAR  LOT*TBO*TYPE 
LOT* TBO* VAR  LOT*TYPE*VAR  TBO* TYPE* VAR 
LOT*TBO*TYPE*VAR ; 

MEANS  LOT  TBO  TYPE  VAR  LOT*TBO  LOT*TYPE  LOT* VAR 
TYPE*VAR  TBO*TYPE  TBO*VAR  LOT*TBO*TYPE 
LOT*TBO*VAR  LOT*TYPE*VAR  LOT*TBO*TYPE*VAR 
/  TUKEY ; 

OUTPUT  OUT=PLOTDATA  P=YPRED  R=YRESID; 

PROC  UNIVARIATE  NORMAL  PLOT; 

VAR  YRESID; 

PROC  PLOT; 

PLOT  YRESID* YPRED; 

PLOT  YRESID*LOT; 

PLOT  YRESID*TBO; 
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SAS 

HALTS !S  Of  VARIANCE  PROCEDURE 
CUSS  LEVEL  IVFQ1AU7I0I 


it: 3.  icokesday  a?r:l  i 


CUSS 

LEVELS 

MitlES 

SET 

M 

123(51 

TBO 

5 

3  4  5  8  10 

LOT 

5 

1  2  3  4  5 

VAR 

3 

1  2  3 

TYPE 

2 

0  i 

MOMaER  OF  OBSERVATION  IK  DATA  SET  *  900 


iv.  im 


35  36 
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source 

DF 

SUM  OF  SQUARES 

HEAR  SQUARE 

t  VALUE 

n  ;  F 

i- SQUARE  C.V 

MODEL 

403 

3086092.21616365 

9146.63080934 

19.91 

0.0001 

0.941794  45.4268 

mot 

490 

227810.79103130 

459.29594966 

100?  KSE 

cos:  EAK 

CORRECTED  TOTAL 

899 

39 13903. 0C7 19496 

21,42.19105 

47. 17739990 

som 

DF 

AIOVA  SS 

F  VALUE 

n  )  F 

LOT 

4 

946113.07232383 

514.98 

0.0001 

TBO 

4 

489507.17354140 

266.44 

0.0001 

TT?E 

1 

0.03925525 

0.00 

0.9926 

m 

2 

35043.05922935 

38.15 

0.0001 

SETIIU'TTPEI 

31 

50760.91525577 

3.57 

O.OOO) 

TBOtlOT 

16 

1779305.47696152 

242.12 

0.0001 

LOT<TYK 

4 

1395.86336107 

0.76 

0.5516 

LOT*  VAfi 

B 

10879.86607234 

2.96 

0.0030 

SETtLOTHU«TT?EI 

124 

99400.84236123 

1.75 

0.0001 

tbo.ttte 

4 

1593.27970625 

0.87 

0.4834 

TBO-VIB 

e 

15280.23991237 

4.16 

0.000) 

SEWBO(Vifi*TTPEl 

124 

12e3f 5. 34069452 

2.11 

C  0001 

BAB'TTPE 

1 

2878.77064521 

6.27 

0.0126 

TBO.LOTiTTfE 

16 

11624.15516630 

: -58 

0.0693 

TBO-LOT.Bil 

32 

55271.14750878 

3.76 

0.0001 

L0T»?4fi«TTPE 

4 

13509  09291719 

7.41 

0.0001 

TBO-?il*TTK 

4 

14363.25119396 

7.83 

0.0001 

no’LOTitU’Tm 

16 

36680.58979438 

5.26 

0.0001 
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SIS  It : 31  KMESDtf.  iHIL  16.  .666  3 

U1LTS1S  OF  fillUCE  PSOCEDUBE 


DEPOOET  IU1UU:  SBOST 


SOOBCE 

or 

SUM  Of  SCDAI £S 

Kit  MMSE 

F  VALUE 

PS  >  F 

i-S0l'AH 

c.r. 

mode 

403 

171562.19431764 

125.  T12640M 

3.69 

0.0001 

0.110815 

60.6670 

SB 801 

4M 

57225.6442458? 

115.311282?! 

BOOT  MS£ 

SBOST  KEAN 

COBBECTED  TOTAL 

8M 

228787.83856351 

10.74124214 

17.64706862 

SOOBCE 

DF 

moil  ss 

F  BALDS 

PB  )  F 

LOT 

4 

30186. 16066813 

86.04 

0.0001 

WO 

4 

16013.01083266 

34.76 

0.0001 

TTPE 

1 

123.1060110? 

1.07 

0.3008 

VU 

2 

16516. 10611106 

80.50 

0.0001 

SETIHMTPE) 

3] 

15010.11310513 

4.22 

0.0001 

WO«LOT 

16 

20311.03800102 

11.00 

0.0001 

LOT'TTPE 

4 

152.21153181 

0.98 

0.4180 

urnni 

8 

1613.06013211 

1.97 

0.0469 

SET*LOT(»U*mi 

124 

31586.11021663 

2.21 

0.0001 

TBO*TTPE 

4 

106.61155651 

1.08 

0.3676 

wo*m 

6 

6062.21100111 

6.5? 

0.0001 

SET'TBOIVU.TIPEl 

124 

20508.03601112 

1.44 

0.0036 

Ell'TTPE 

1 

1266.88733000 

10.00 

0.00 10 

TBO<lOT<m 

16 

1116.31153801 

0.62 

0.6679 

TBO'WWU 

32 

1201.18255351 

1.16 

0.2513 

L0T«1M>TTPE 

4 

612.33112821 

1.46 

0.2)42 

TBOHmUPE 

4 

002  00221611 

2.15 

0.0734 

WO*LOT<»U*nPE 

16 

1560.38107011 

0.85 

0.6336 
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SAS 


It:-,  ftMESIiAT.  AHI-  10  Mi  4 


halts!  $  of  vajiaics  procedure 


DEFE1DER  VARIABLE:  PEIS  HOB? 


SOURCE 

DF 

SUV  Of  SQUARES 

l£it  SQUARE 

MODEL 

403 

261450.71276121 

606.36886541 

ERROR 

406 

114574.00837702 

230.00508463 

COBBECTED  TOTAL 

600 

306024.72113013 

SOURCE 

OF 

ASOVA  SS 

?  VALUE 

PR  >  F 

LOT 

4 

20555. 52434BOO 

22.25 

0.0001 

T60 

4 

4602.80351021 

4.06 

0.0006 

TYPE 

1 

7677.56636450 

34.10 

0.0001 

VAI 

2 

52003.03702266 

112.56 

0.0001 

SETIVAMYPE) 

31 

3B044. 011 34356 

5.44 

0.0001 

TBO*LOT 

16 

20306.04038346 

5.49 

O.COOl 

LOT* TYPE 

4 

633.07670372 

0.60 

0.6016 

LOT* VAR 

8 

3066.40358621 

1.66 

0.1056 

Sn*LOTlVAR«TYPEi 

124 

43655.56215600 

1.52 

0.0009 

TBO»TYPE 

4 

1120.23150214 

1.22 

0.3004 

TBOtVAR 

8 

5343.62661746 

2.89 

0.0037 

SET*TBO(YAI»TY?E) 

124 

70060.30072504 

2.45 

0.0001 

VAR* TYPE 

1 

0.00000000 

0.00 

1.0000 

TBO*LOT«TYPE 

16 

3305.08393570 

0.60 

0.5761 

TBO«LOT«7AR 

32 

6903.08742218 

0.03 

0.5739 

.07* VAR* TYPE 

4 

309.68028756 

0.43 

0.7651 

iBO*VAR»TYPE 

4 

4700.05004523 

5.00 

0.0005 

T60«L0TtVAR»TYPE 

16 

3713.50679050 

1.00 

0.4409 

F  VALUE 
3.02 


PI  >  F 
0.0001 
BOOT  USE 
15.10655206 


i-SQUARE  C.V. 

0.710600  57.6525 

PERS SOB?  KEAV 
26.27121767 
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S  AS 


16:3!  ME&ESI'IY.  iPEIL  10.  .969  5 


ANALYSIS  OF  VAB1AICE  PBOCEDOBE 


DEPEIDE1?  VARIABLE:  DELEBB 


SOUBCE 

DF 

SUM  OF  SQUARES 

£A»  SQUARE 

MODE 

403 

430070.22330010 

1069.52000671 

EUOB 

400 

60054.58264467 

179.54552003 

COBKECTED  total 

600 

526133.60015307 

SOUBCE 

DF 

UOM  SS 

f  pile 

PE  >  F 

LOT 

4 

1 306 . 584 13127 

1.82 

0.1238 

TBO 

4 

&  1306  18270409 

127.14 

0.0001 

m 

1 

19(1.08912491 

10.64 

0.0011 

m 

2 

JOS! . 1569685! 

5.71 

0.0035 

SETIMBiTTPE) 

31 

1(1211.50648509 

25.37 

0.0001 

TBO»LOT 

16 

2300. 41726(41 

1.15 

0.3065 

LOTiTYPE 

4 

1552.03621285 

2.16 

0.0123 

LOT*  MI 

6 

1166.05613204 

1.24 

0.2115 

SET  •  LOTI  MI*  TYPE  1 

124 

39328. 33911134 

1.77 

0.0001 

TBOiTTPE 

4 

1444.48697099 

2.01 

0.0911 

TB0*V11 

6 

8039.00529322 

6.20 

0.000! 

SET*T50(MI»TTPEI 

124 

65304.65735081 

3.13 

o.ooo; 

MB-TTPE 

1 

28815.72348101 

161.05 

0.000) 

TBOiLOTiTWE 

16 

4459.82212102 

1.55 

0.0716 

TBOiLOTiMB 

32 

6196.29688461 

1.06 

0.3553 

LOTUiBiTTPE 

4 

216.24212807 

0.30 

0.8772 

TBO>SUi?T?E 

4 

20702.15150809 

26.63 

0  0001 

TBOiLOTiMI'TTPE 

16 

0.00000000 

0.00 

1.0000 

PB  >  F 
O.OOOl 
BOOT  ICE 
13.30040000 


R-SODABE  C.V. 

0. 03 1370  60.6640 

Dam  KEiH 
16.57016404 
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JE'EKEr  TA2UELE:  PSSESI 


SOURCE 

5F 

s'jk  cf  33': ui: 

ME  AS  SCLARE 

F  VALUE 

•t  F 

••e:':aee  :  v. 

MODEL 

403 

424949.80901476 

;C54.466Cr36 

26.76 

0.000; 

0.956033  12. 6495 

£i;OH 

496 

19542. 05598C10 

39.40:32456 

ROC?  M£E 

FEEERE  MEAL 

CORRECTED  TOTAL 

899 

444492.66600466 

£.27704744 

45.8:95552: 

SOURCE 

OF 

ASOVA  SS 

F  YAi.DE 

fl  :  F 

LOT 

4 

234.71011152 

1.49 

0.2042 

TBO 

4 

12250. 02582361 

71.73 

(  000! 

TYPE 

1 

12817.17504600 

325.30 

0  0001 

VAR 

2 

219300.40496487 

2742.91 

(.0001 

SE7<VAB*TYPE) 

31 

156217.40130964 

129.53 

0.0001 

TBO* LOT 

16 

841.88154835 

1.34 

0.1105 

LOT 'TYPE 

4 

226.87175761 

1.44 

0.2191 

LOT* VAR 

8 

589  60354735 

1.81 

0.0525 

SET* LOT (VAB* TYPE) 

124 

8082.94662557 

1.65 

0.0001 

TBO»TYPE 

4 

32.96368533 

0.21 

0.9333 

TBOiVll 

8 

3465.19045606 

10.99 

0.0001 

SET*TBO(VAR*TTP£l 

124 

14380.82704713 

2.94 

e.mi 

VAS'TYPE 

1 

0.00000000 

o.oo 

1.0000 

TBO«LOT»TTPE 

16 

875.81329835 

1.39 

0,1416 

TBO* LOT* VAB 

32 

3083.52190212 

2.45 

0.0001 

LOT*VAR*TY?E 

4 

14.20115316 

0.09 

0.9855 

?BO*VAB*TYPE 

4 

2455.22493057 

15.58 

0.0001 

TBO*LOT*VAB«TYPE 

16 

0.00000000 

0.00 

:  .oooo 
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sas  it-ii  «DK£.r'iy.  i?s:;  i 

HALTS  IS  OF  VAR1A1CE  PROCEDURE 

TUKEY'S  STUDEPIZED  IAIOE  iHSDl  TEST  FOP  VARIABLE:  COST 
BOTE*  THIS  TEST  COPROLS  THE  TYPE  1  EXPEBUEKTVISE  EHBOR  BATE 
EOT  SEVERALLY  BAS  A  HIGHER  TT?E  11  EIBOR  BATE  THAV  RE3*» 

ALPBAr0.05  DF>4M  USE1 4 50. 296 
CRITICAL  VALUE  OF  STCDEHTJZED  BAI0E*3.872 
Id VI HUB  S1GVIFICAT  D1FFRSEVCE*6.185I 

KEHS  KITH  TEE  SAME  LETTER  ABE  MOT  SIGIIFICAPLY  DIFFEBEP. 


GROUPING 

NEAV 

V 

LOT 

A 

110.333 

160 

5 

B 

45.125 

160 

1 

C 

20.362 

180 

3 

C 

C 

27.652 

180 

2 

c 

c 

23.415 

180 

4 

0.  1065  ? 
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$AS 


-v  tOkEil'AY.  APE 


MALTS!  S  OF  VARiAICE  PROCEDURE 

TUKEY’S  STUDEITIZED  RAISE  (BSDi  TEST  FOR  VARIABLE:  SHORT 
IOTE:  THIS  TEST  COITBOtS  THE  TYPE  I  UPEBIHEITWISE  ERROR  RATE. 
BCT  QEfc'ERALLT  HAS  A  HIGHER  TYPE  II  ERROR  RATE  THAV  REMfc 

ALPHA- 0 . 05  DFM9C  K£E=  U5.374 
CRITICAL  VALUE  OF  STOIiETIZE!'  RAISED. 672 
HI II IKH  SIGHFICAT  DIFfEB£IC£:3.0999 


MEAIS  WITH  THE  SAME  LETTER  ARE  HOT  SI3KIFICAHTLY  5IFFEREIT. 


GROUP! 10 

WEAR 

I 

A 

24.557 

160 

A 

B  A 

2i.549 

160 

B 

B 

19.666 

160 

C 

14.010 

180 

0 

6.253 

160 
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SAS  It:31  «1'I££64Y  APEi: 

AIALYS1S  OF  VARIAICE  PEOCEl'CEE 

TDXEY'S  ST0BEJTTJ2SD  HAMGE  (HSI'i  TEST  FOE  VARIABLE:  PEHS80ST 
IOTE:  THIS  TEST  CQITBOLS  THE  TTPE  1  EXPEBUCHTWISE  ERBOE  HATE. 

BUT  SEME  HALLY  HAS  A  H1GHEB  TYPE  11  EHBOB  HATE  THAJi  HEW 


ALPHA* 0.05  DF*4tt  KE*230.W6 
CRITICAL  VALUE  OF  STTOEHTIZSD  RAMGE=3.872 
Kill  MUM  SlOliFlCAVT  D1FFEHEICE*4.3863 

IE  AMS  WITH  THE  SAME  LETTER  AHE  10?  SIGIIFICAITLY 


GROUP I 1G 

ME  AM 

X 

A 

32.614 

160 

A 

B  A 

29. 191 

160 

B 

B  C 

26.956 

180 

C 

C 

23.535 

160 

0 

16.661 

180 

D1FFEHEIT. 

LOT 

4 
2 
3 
1 

5 


i.  .96S1  k 
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$15  li:i.  HEMI-DA T.  APRIL  IP.  >89  if 

HILTS  IS  OF  TAB]  MCE  PROCEDURE 

TWIT'S  STUDENT! ZED  BAKE  <&SU=  TEST  FOE  TAB! ABLE:  DELERR 
WTE:  THIS  TEST  COITBOLS  THE  TYPE  I  EIPERIHEHTHISE  EBBOE  RITE 
BUT  QENEBAU.Y  HAS  A  HIGHER  TYPE  I!  EBBOE  BATE  THAI'  REOVQ 

AL?HA*0.05  DF*496  BE*  W.  548 
CRITICAL  VALUE  OF  STUDEMTIZEL  BAME=3.872 
HI  II I  HUH  SIGIIFICAIT  D1FFEREICE=2.6671 

KEA1S  HITS  THE  SAC  LETTER  ABE  HOT  SIOKIFICATLY  DIFFESEIT. 

TUEEY  GBCUPJ  KG  CAH  H  LOT 

A  It. 223  160  4 

A 

A  17.490  160  2 

A 

A  16.421  160  1 

A 

A  15.960  160  5 

A 

A  14.75?  160  3 
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SIS  18:3:  BMESMT  IPEII.  16.  IMS 

UliTSlS  OF  ViBlilCE  PB9CESUBE 

TIBET'S  STIJBEIT’Zn  B4I3E  lESOi  TEST  FOB  IIB11BLE:  FEBEBB 
MOTE:  THIS  TEST  CQIT80LS  TIE  TYPE  I  EXPEBIKEMTV1SE  EBBOP.  BITE. 

BET  GEIEBAU.T  HIS  t  IISKEB  TYPE  II  EBBOB  BITE  THIS  BE9M0 

ILPU>O.OS  OF'4M  MS£=34.4013 
CBITICAL  VIUJE  OF  ST0DEMTI2E0  BUOEO.BY2 


minimum  sigiificat  diffem»ce=  i .ei u 

tens  WITH  THE  SAME  LETTEB  ABE  MO?  SI MIFICATTLY 

OIFFEBEI? 

TITKIY  GBOUPIMO  ME  AS 

K 

l« 

A  46.2964 

i 

16f 

2 

A 

A  46.276? 

• 

iec 

4 

a 

A  46.9421 

» 

18C 

5 

A 

1  45.0T13 

180 

1 

1 

1  44.9093 

18C 

3 
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US  lt:3!  KMESiAT.  APRIL  IS  Mi 

AIALTS1S  or  U1IAICE  PROCEDURE 

TIBET'S  STroOTlZH)  RAKE  IHSD!  TEST  FOR  VARIABLE:  COST 
I0TE  THIS  TEST  COITIOLS  THE  TTPE  I  EBERIICITWISE  ERROR  RITE. 

BET  GEBERALL!  HIS  1  HIGHER  TYPE  II  ERROR  RITE  THAR  REGNO 


ALPHAl0. 05  OPMOf  KSEM5S.295 
CRITICAL  VALUE  OF  STOPESTIZET  RAIGE=3.615 
H1IINUN  S1GS1FICAIT  DIFFER£ICEI6. 1651 

CASS  KITH  THE  SAJE  LETTER  ARE  I0T  SIGIIFICAITLT  OIFFEREHT, 


GBO'JMIG 

KEAI 

V 

TBO 

k 

91.032 

140 

10 

B 

45.206 

160 

8 

B 

B 

41.132 

160 

2 

C 

20.501 

160 

6 

C 

C 

2B.025 

180 

4 
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S>£  :<:2.  REOEEEIAY 

UiLYSIS  Of  MB1UCE  FIOCEL'UIE 

TUKEY'S  STKiEriZEr  UISE  <ES!)i  TEST  fOE  VlilliLE  SHORT 
KITE:  THE?  TEST  CorSOlS  THE  TYPE  1  EXPES1EEITVISE  EBIOl  RITE. 

BUT  5EHEEA1.L*  HIS  I  HIUEEi  TYPE  I!  EEiOt  BATE  THUI  tit* 


ALPHA-0 .05  DF'IM  HSEMI5.3H 
CRITICAL  VALUE  OF  STUtERTIZED  RAHOE-3.BT2 
BIIIIMI  SISH1F1CUT  DIFFEEEICE=2  09M 

IEAJS  (ITS  THE  SUE  LETTS!  ARE  HOT  SKUIFICUTLY  MFFEBEHT. 


GiODPIIG 

MEM 

II 

TBO 

i 

23.075 

1 80 

10 

B 

19.616 

180 

6 

B 

B 

IB. 610 

1B0 

6 

B 

B 

16.622 

ieo 

4 

C 

10.312 

180 

2 
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sas  it:i.  eaiEi'ty  mi.  is.  :ssv  n 

AHALT'IS  OF  VAF1AKE  MOCEWHE 

TOKET •  S  STBDE8TIZE0  SAKE  (HSP  TEST  FOE  IAH83LE;  PEBSHOBT 
ion  THIS  TEST  comoii  THE  TYPE  1  EBEBIIttTVlSE  E88QE  8ATE. 

BUT  QE8E2ALLY  HAS  A  EI5HES  TTEE  I!  E8801  HATE  THAI  ujW 

AEPHJ--0  05  bf-iK  HSE>230  098 
CS1TICAL  VALUE  OF  STOOEriZEO  8AI3£=3.S72 
IIIIHDH  S1GI1F1CAIT  DIFFEiElCE-4.3863 

SAKS  MTB  THE  SUE  LETTES  AH  HOT  SIMIFICAHTLT  t’FFEHEIT 

TCKET  GHOTHIf  HE  AH  H  TBO 

A  30.003  180  10 

A 

8  A  07.801  180  « 

8 

8  28.440  180  2 

8 

8  24.211  180  8 

I 

8  24.082  180  8 
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sis  !«-3!  KMES51T.  tftl.  It.  ;48S  1 

UILTSIS  Of  liEUICE  FSOCEOOEE 

TUXfT'S  STUOEHTtaO  SHOE  'BO!  TEST  FOE  UIIEBLE:  OELESi 
DTI:  THIS  TEST  CMTSOIS  TEE  TTEE  I  EE?£SI*MSE  EEEOE  SITE. 

EOT  3EIEtU.1T  US  1  HIQHEI  TTEE  11  EEEOE  EETE  TUI  EEOHQ 

ILfU-O.OS  DF>4M  HSE'lTg  .546 

CE1T1CEL  JEi.CE  OF  stoeitize:  SEIOEO.STJ 

mufti  SIGIJFICUT  DIFFEEEKE'2  StTl 

mm  an  the  sue  iettee  mi  kt  sioeificuttit  oiffeeeit. 


ooufuo 

til 

1 

TBO 

i 

25.58$ 

180 

10 

B 

25.408 

180 

8 

C 

16.538 

180 

6 

D 

8.938 

180 

4 

E 

2.376 

180 

2 
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SIS  16:3.  KOSISMT  Afi:i  :i  :M 

AIALTSIS  OF  IAUAKE  PK-CEDM 

TUKET'S  STTOEITIlir  UK!  iBSD'  TEST  FOE  VA11ABLE:  FEEEM 
■OTE:  THIS  TEST  COITEOLS  THE  TYPE  !  nPQllETOS!  EliOH  UTE. 

HOT  GEIEULLY  HAS  A  EIQHEB  TTTE  I!  EiiOB  UTE  THAI  SEOK 

ALPEA*C.05  OF-4K  BE=3«  40)3 
CEITICAL  VALUE  OF  ST 05 DT! lit  UM£=3.!T2 
KII1HW  SiaiFICAT  OIFTnnCEM.Slie 

ICAIS  HITE  TEE  SA£  LETTE1  AEE  EOT  SIOIFICAITLT  OIFFEEEIT. 


GBCtJPIlG 

mu 

» 

TB0 

1 

46.0441 

160 

10 

A 

A 

49.4533 

160 

8 

B 

46.2735 

160 

6 

C 

43.1060 

160 

4 

d 

40.3166 

160 

2 
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SiS  £6:3!  KMSSOIT  tfEIL  IS.  !»Si 

UILTS1S  OF  IttIUCE  'IOCE50HE 

TDxiT's  srmrizK  iuge  ibsd>  test  foe  ububle:  cost 

VOTE ■  THIS  TEST  COETBOLS  TEE  TTPE  ]  EEPS81IMT»!SE  EESOt  HITE 
BUT  GEKULLT  US  i  BIGHEE  TTPE  1!  ESiOB  BATE  THU  EE6U0 

BLPHU0.05  DF>40G  IGE»454.2M 
CIITICiL  WOE  OF  STUDEITI2EI  IUG£'2  HO 
MIIIHDH  SIGI1F1CUT  DIFEEmCE'2.0221 

HUHIIG:  OEU  SIZES  UE  HOT  EOS1L. 

HUV9I1C  KEU  OF  CELL  SIZES’415.278 

CUE  HITS  THE  SUE  LFTTEJ.  BEE  HOT  SIGHIFICUTLT  IIFFEBir. 


TCKEY  QBOUPJIG 

HI 

1  TTPE 

i 

IT.  EM 

575  0 

k 

1 

IT. 169 

325  1 
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SAE  KEI'SE*  1AT  APE 

AULTS  IS  OF  UBIAKE  PBOCED’JK 

WEST'S  STCDETIZEI  BAIGE  '.SSC'  TEST  FOE  MKULE:  SHOK 
IOTE  THIS  TEST  COITHOLS  THE  TTPE  I  EXmiKEPKISE  EB5CE  KATE. 

BUT  OEIEBAILT  HAS  A  HIOHEE  TTPE  II  EUO&  BATE  THAI  UM 

ALPHA =0.05  DF*4K  ISE* 1 15.374 
CRITICAL  VALUE  OF  STU9EITI2EC  IAIW*2.T7« 

Kill  KM  SIGIIFICAT  OIFFEBEICE  4046 

•ABi:  10-  CELL  SIZES  ABE  KT  EQUAL. 

HABKOKIC  £AS  OF  CELL  SIZES** 15.276 

CAMS  KITE  THE  SAKE  LET7E2  ABE  HOT  SlflUFICAIW  DIFFEBEIT. 

TDIET  ffltOOPIM  £AI  I  TTPE 

17.9250  575  0 


A 

A 

A 


17.1536  325  1 
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SIS  !t:3!  KtSESlAT.  AfECL  .y  iySy  H 

UbTSIS  OF  VAEIAKE  PSOCEttJEE 

toit's  stoeetizei  iajke  <bsd>  test  foe  vasiaele  pebsbobt 

I0TE:  THIS  TEST  C01TS0L5  TEE  TYPE  I  EEPEEU6Htr.SE  ESEOS  1ATE. 

BUT  GEIEEALL*  SEE  i  BIOSES  TTPE  I.’  EISOE  SITE  THAI  EESK 

ALPHA*0.05  OFMM  «SE=!J0.M6 
CEITICAL  VALUE  OF  ST!I0EIT!2EP  EAME=2.778 
UI1HM  SIOIIFICAIT  DIFF£BEICE=2 . 0T2S 

VAUIH:  CELL  SIZES  AEE  I0T  EQUAL. 

HABIOIIC  HEAI  OF  CELL  SIZESM1E.27S 


KUI£  «TH  TEE  SUE  LETTEE  AEE  MT  SIMIFICAITLY  OIFFEEEIT. 


TUXTT  CB0r?IIG 

mi 

i  rm 

A 

28.495 

5*»5  0 

B 

22.336 

325  1 
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SAS  it  3;  KLKil'iY.  LlE 

UALTSIS  OF  YAB1AICE  PBOCEDUBE 

TUKIY'S  STTOErilED  BAKE  <HSD'  TEST  FOE  VA51ABLE  DELEHE 
IOTE:  THIS  TEST  COFTEOIS  THE  TT7E  i  EIPEBIIttWSE  EEBOE  BATE. 

BUT  OEXEEALLT  HAS  A  HI3HEB  TYPE  II  EBBOS  BATE  THAI  BEWG 

ALPHA-0 . 05  DFHM  KSE-179.546 
CRITICAL  VALE  0?  STM1T12EI  kAXQ!=2.77g 
MIKIMUM  S1QHIFICAXT  DIFFEBEICE5! .827 

■AM  I  IQ-  CELL  SIZES  ABE  1C!  EODAL . 

HAfiHOi::  MEM  Of  CELL  SIZESMJ5.276 


MEAIS  MIT5  THE  SAME  LETTER  AEE  HOT  SIQKIFICAKTLY  DIFFEFIK! . 


TUKEY  QBOUPIIQ 

A 

B 


MEAI  X  TYPE 

18.5256  325  1 

15  4644  575  0 
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S AS  IC:’l  WEMESIlAT  AFBIL  IS  !46‘  2: 

A1IALTS1S  Of  VAB1AICE  PEOCEOCRE 

TWEE'S  STODE1TIZEO  RAM!  (ESDI  TEST  FOE  IABIABLE:  PERIEJ 
IOTE;  IBIS  TEST  COmOlS  THE  TTPE  I  EEPERIHERTV1SE  ERROR  RATE. 

BUT  GEBERALLY  HAS  *  HIGHER  THE  I!  ERROR  BATE  THAI  REGWQ 

ALPHA-0.05  0FMB6  I6E*M.401J 
CBIT1CAL  YAL0E  OF  STDDEHTI2ED  RAIO£=2.77S 
H2HIMUH  S1MIFICAIT  OIFFEIUCE- .  BS58T 

WARHIIG:  CELL  SIZES  ABE  HOT  EQUAL 

BAUDIIC  MEAI  OF  CELL  S!ZIS=415.276 


KEAKS  WITH  THE  SAME  LETTEK  AH  ACT  SIQKIFICATTLT  fcIFFESSK 

TOKST  GROUP 1 IQ 

£111 

I  TTPE 

1 

46.£5«7 

575  0 

G 

40.6000 

325  1 
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s*s  !{;j:  ig*sm*.  in::  u.  mu 

AIALTSIS  OF  VAS.'AKE  PSOCEDCET 

TWET'S  Sn®E«?:ZIt  IAIGI  I  ESDI  TEST  FOE  V1S1ABLE:  COST 
I0TE:  THIS  TEST  COITtOLS  TEE  TTPE  I  EEFEE1KEVTV1SE  EiKiE  EATE 

ILF8A=0.05  CO»F1DEICE=0.«5  DF--M6  KE‘455.296 
CIITICAL  VALUE  OF  STUOEF1ZEL  illOE'l.JH 

COIMS1SOIS  SIM1F1CAT  IT  TE  0.0!  LEVEL  Ai£  IltICITEi  BY 


SIMULTANEOUS 

SIMULTANEOUS 

10«B 

DIFFERENCE 

l'??Ei 

m 

CONFIDENCE 

SETMEEK 

CONFIDENCE 

COMPABISOM 

LIMIT 

MEANS 

LIMIT 

3 

-  2 

6.932 

14.367 

19.602 

3 

-  1 

13.  664 

16.702 

23.740 

2 

-  3 

-19.602 

-14.367 

-6.932 

2 

-  1 

0.553 

4.335 

8.117 

l 

*  3 

-23. MO 

-16.702 

-13.664 

1 

-  2 

-8.11? 

-4.335 

-0.553 
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SAS  16:3.  KDIESt'AY  APRIL  lfr  1989  23 

HALTS  IS  OF  VAEIAICE  PROCEDURE 

TDKEY'S  STUDE1TIZED  HAKE  (BSD)  TEST  FOB  ¥  Ail  ABLE:  SHORT 
BOTE :  THIS  TEST  COITIOLS  THE  TYPE  I  EXPERIMEITKISE  ERKOR  KATE 

ALPHA* 0.05  COIFID£1CE*0.95  DF*4tt  KSE--115.374 
CRITICAL  VALUE  OF  STMKTIZED  KAKEO.325 

COMPARISONS  S10IIFICA1T  AT  THE  0.05  LEVEL  AKE  IIMCATED  BY 


VAE 

COMPARISON 

S1K0LTAIE0DS 
LOVER 
COA7I DEICE 
LIMIT 

tIFFESEICE 

(ITKE* 

KEUS 

SIMCLTAIEOCS 
UPPER 
COAT.’ DEICE 
LIMIT 

3 

-  2 

6.2010 

4.0148 

11.7381 

3 

-  1 

10.8952 

13.4202 

15.9453 

2 

♦  3 

-11.1387 

-4.0148 

-6.2910 

2 

•  I 

2.5097 

4.4054 

6.3011 

l 

-  3 

-15.9453 

-13.4202 

•10.8952 

l 

-  2 

-6.3011 

-4.4054 

-2.5097 
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us  16:3*  CMBSi'ST.  IftEL  !*  :w* 

UilTSIS  OF  H1IUCE  PBOCEBUHE 

TUXET'S  STU5EITI2E1  HUGE  iHSEI  TUT  FOI  UEI9BLE;  peeshoit 
KTE;  THIS  TEST  COITJOLS  THE  TTPE  I  EXPEIIWTV1SE  ESiOB  HiTE 

ilpm=;.o‘  cor : mice-- o. !i  OF'«e  CE^sc  soe 
C1ITICAL  VALUE  OF  STUtEITIZU  1AME-3.325 

corisisor  SIOIIFICUT  iT  THE  0.05  LEVEL  UE  IIDIC1TK  ST  '•••' 


SIKILTUEOIIS 

SIMJlTilEOOS 

LOW 

DIFFEBEICI 

0PPE1 

m 

COHFIKICE 

BFTWSEI 

COIF 1  DEUCE 

Col£r  AE1  SOB 

LIMIT 

MUMS 

LIMIT 

3 

-  2 

11.4351 

15.2893 

19.1434 

3 

-  i 

IS. 9225 

22.4954 

26.0683 

2 

-  3 

•19.1434 

>15.2693 

-11.435! 

2 

-  1 

4.5239 

7.2062 

9.8685 

1 

•  3 

-26. 0585 

•22.4954 

•16.9225 

1 

-  2 

•9.9985 

-7.2062 

-4.5238 
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S AS  16:31  CMTESliA?  W..  : 

ANALYSIS  OF  fiRWKE  Pi OCEDt/K 

TUKET’S  STODEIT1ZEI)  BAKE  (BSI»)  TEST  FOB  Fill  ABLE 1  DElEiB 
MOTE:  THIS  TEST  COfTBOLS  TEE  TT?E  1  EXPEilICBTWISE  EB10B  BATE 

ALPHA =0.05  COIF  JIiEICE*  0.95  DF*406  MSI* ITS. 546 
CB1T1CAL  TILDE  OF  STUOET1ZEO  BAIQE^.325 

COM? Ail S01S  SIQIIFICAIT  A?  THE  0.05  LEVEL  ABE  IWDICATET  BT 


SIMDLTAIEOUS 

S1MCLTAIE0DS 

LOVES 

D1FFEBEMCE 

UPPEi 

m 

COM?  IDEIICE 

BETtf-EII 

COEFIDEMCE 

COMP Ail SOI 

LIMIT 

MUMS 

LIMIT 

3 

-  2 

0.7936 

4.1915 

7.5895 

lit 

3 

-  1 

1.2925 

4.4424 

7.5924 

Hi 

2 

-  3 

-7.5805 

-4.1915 

-0.7036 

•  •• 

2 

-  1 

•2.1140 

0.2509 

2.6157 

i 

•  3 

-T.5924 

-4  4424 

-1.2025 

1 

-  2 

-2.6157 

-0.2509 

2.1140 

0  1080  25 
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SIS  16:3;  (EMESiAT.  AFE1L  19.  1989 

U  ALTS  IS  OF  VAJIIAICE  PBOCEDBEE 

TOKEY ■  S  STTOEBTIZED  SAIGE  (BSD)  TEST  FOE  FAS! ABLE  PEREIB 
KTt:  TS1S  TEST  C0IT10LS  TSE  TYPE  1  EEFElUeiflCSE  EUOB  SATE 

ALPBA'C  0!  COEIDEECE--O.S5  DF'496  ICE' Si  4013 
C11T1CAL  I ALOE  OF  STDDEHTIZED  RAISE'S  325 

COMFAS1SOIS  SiaSlFlCAT  AT  THE  0.05  LEVEL  AIE  IIDICATEP  BV 


VAB 

COMP  AI! SON 

sucimnous 

LOMEB 

COMFIDEICE 

LIMIT 

OIFFEEUCE 

BETIEEI 

£AIS 

SIMULTANEOUS 

DPPEB 

COBTIDEBCE 

LIMIT 

3 

-  2 

32.3536 

33.9455 

35.5373 

3 

•  1 

46.1264 

45.6050 

46.0777 

2 

-  3 

•35.5373 

-33,9455 

-32.3538 

2 

-  1 

11.5487 

13.6565 

13.7643 

1 

*  3 

-48.0777 

-46.6030 

-45.1264 

1 

-  2 

-12.7643 

-12.6565 

-i 1.5407 

147 


SIS  je-3:  KDIESDI!  I till  1*  1999 

IIALTSIS  OF  VUIUCE  PROCEDURE 

CANS 


T80 

LOT 

1 

COST 

SH05T 

pershort 

DELEU 

PEEEU 

2 

1 

36 

3b. 35 1202 

11.5868706 

27.6259877 

3.1506861 

40.7559350 

2 

2 

36 

45.068724 

8 . 2236613 

21.6120518 

1.4616322 

4C  0517133 

2 

3 

36 

41.637606 

11.6052362 

25.6601053 

3  4516604 

40  6616026 

2 

4 

36 

30.587900 

10.9064004 

25.2490050 

1.3991*51 

39.1407972 

2 

5 

36 

41.015926 

0.1463710 

24.0202803 

2.4092640 

40.3542586 

4 

1 

36 

29.613956 

12 .7320075 

22.128'T416 

1C. 4 183065 

43  1647575 

4 

2 

36 

29.166954 

16.6790863 

28.-220869 

6.2017279 

42  5987644 

4 

3 

36 

26.660043 

!  6. 3*74902 

22.4658052 

fc  4512304 

42.10*7975 

4 

4 

36 

22.110066 

20.0525064 

20.2517525 

10.2022562 

42.3655681 

4 

K 

36 

32.546886 

14.3568956 

24.6356208 

6.4141020 

42.0640456 

6 

! 

36 

36.870735 

14.0820004 

17.7349140 

15  8985332 

4!  5631914 

6 

2 

36 

22.513733 

22.6783861 

26  5706106 

1*V58,542I 

4"  4626722 

6 

3 

36 

24.535388 

16.6623167 

23.6119099 

15.4453766 

45. 0573 i  00 

6 

4 

36 

19.666236 

26.1516047 

32.3546905 

20.9582918 

47.8711109 

6 

5 

36 

43.268506 

1) .4764432 

20.7745380 

12.6003435 

44.6034000 

6 

1 

36 

46.045366 

15.2200011 

18.262624! 

22  5282606 

47.5732508 

B 

2 

30 

22.264563 

25.0883553 

31.0475554 

26.1608404 

50.2401910 

fi 

3 

36 

27.932226 

21.2023257 

51.2731321 

23.0071438 

48.0390100 

6 

4 

30 

10  449422 

29  8533942 

21.7745999 

25  7930292 

40.8674560 

B 

i 

36 

111.31944! 

6  7497322 

17.3331541 

27.5418358 

50.9307650 

10 

I 

36 

24. 74  .’920 

10.4237752 

3ck:»x 

2C.U1M  5 

51  6091047 

10 

2 

30 

18  104296 

34.1702801 

37.60C2202 

3!  0091 103 

5C.  109*1353 

10 

3 

30 

27.032730 

29.3212882 

31  6428205 

23.4262855 

4*.  65004  IS 

10 

4 

30 

16. 263381 

34.9200342 

42.4370074 

32 . 759*7800 

5C  516*253 

10 

5 

36 

323.010550 

0.5343915 

7.5319063 

30.6301549 

5:  5319711 

LOT 

TT?E 

1 

COST 

SHORT 

PERSHOIT 

DELEU 

PEIEU 

1 

0 

115 

43.270024 

14.8156566 

35.5573693 

15.2574145 

48.5217337 

1 

65 

48.30D002 

12.565362] 

10.0563306 

16.4606570 

40.6261306 

0 

115 

28.076302 

21.7104675 

31.2755307 

16.5695387 

40.2044003 

1 

65 

26.806655 

21.2478440 

25.5021644 

10.0630507 

40.0076766 

0 

115 

29.568842 

20.3116004 

30.2502503 

13.2643606 

47.5000052 

1 

65 

26.060157 

10.0767320 

21.1276352 

17.3616836 

40.3240700 

0 

115 

24.063500 

23.8606455 

34.1635841 

15.6083463 

46.5910562 

1 

65 

22.268770 

25.7890233 

30.4254802 

22.6683172 

42.1805666 

0 

US 

110.002056 

6.0220331 

21.2305810 

16.4021020 

40.3744530 

1 

65 

100.327463 

7.0607877 

14.6664231 

15.0103606 

39  6694423 

LOT 

TU 

1 

COST 

SKIT 

PQSB01T 

DELEU 

mm 

1 

100 

36.040171 

10.4301340 

16.2016460 

15.5406420 

35.2260771 

2 

55 

46  678455 

17  0631200 

26.2607320 

17.5644117 

46.6701460 

3 

25 

74  004226 

21.5766230 

46.4003644 

17.3467400 

60.3046826 

1 

100 

25.861700 

17.9250780 

24.3463727 

15.9092632 

35.4119560 

2 

55 

26.014043 

21.7665711 

32.1724424 

17.7106345 

49.0620656 

3 

25 

36.431874 

35.5501764 

42.0082060 

23.3074812 

63  7642036 

1 

100 

26.524868 

16.5054616 

21.0655850 

14  8556350 

35.0566523 

2 

55 

28.353247 

10.7504256 

27  7668521 

14.2001070 

47.2847678 

3 

25 

42  026468 

33.5604694 

45.120:660 

16.1677480 

70.0036476 
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SIS 


16:2.  SMV  API!:  iv  108*  26 


UALTSIS  OF  f  All  MCE  PkOCEDUIE 

£ilS 


LOT 

Til 

1 

COS? 

SHOt? 

PEISSOS? 

DELEU 

PEHESB 

1 

100 

21.735546 

21.1604820 

20.5038136 

H.OB5M53 

36.5013256 

2 

55 

23.186107 

26.2467084 

32.7371812 

I5.IM1T78 

47.6476520 

3 

25 

30.634740 

34.4257400 

46.2216624 

21.5158456 

82.3623:72 

1 

100 

105  214565 

6.1542216 

12.6174220 

14. 1621535 

35.2019325 

2 

55 

112.616401 

0.3834076 

21.8074668 

15  865886! 

47.6077382 

3 

25 

125  786387 

14.1624164 

37.3524583 

23.5351481 

84.7042600 

m 

TYPE 

1 

COS? 

SHOT? 

PEISB01? 

DElEfil 

FEIEtt 

i 

0 

300 

41.5524285 

14.3540646 

21.1855728 

10.4273310 

34.3050300 

i 

1 

200 

45.8008260 

14.5615458 

20.5640610 

23.0603071 

37.1055773 

2 

0 

150 

46.1200251 

16.7034131 

30.6611387 

21.3720770 

40.3230077 

2 

1 

125 

40.3426766 

21.3012867 

25.1301303 

0.6327150 

46.7)10101 

3 

0 

125 

61.0571411 

27.8572846 

43.4403082 

20.3160520 

82.0616262 

?BO 

?TPE 

■ 

COST 

SKIT 

rusniT 

DELEU 

fEilli 

0 

m 

42.35 13970 

10.6265609 

26  9136661 

2.79C1226 

13.0526959 

I 

55 

J8.9T5UT6 

9.3973062 

19.2192957 

1.8591992 

35  1716)91 

0 

:n 

2!  7258621 

15.6363511 

29.1172102 

8.3021069 

45.8324511 

1 

55 

25.7850057 

16.5966313 

21.3870793 

10  0599559 

38.2877854 

0 

115 

26.2918011 

19.828275) 

25  ’796293 

15.3837928 

19.2387689 

1 

55 

31. 8691951 

18.1552185 

21  1311709 

II  5801083 

11.0326991 

0 

115 

11.1166057 

19.6975162 

27.0009367 

23  5829931 

52  5065196 

1 

65 

17.0811581 

19  1720381 

16  9181219 

28.6369235 

11  0176098 

10 

0 

115 

92  3951885 

22  6367581 

31.3356307 

27.2616572 

52  6199997 

10 

1 

65 

91.1119591 

23.8175737 

27.6997130 

33  7059021 

15.1567785 

?B0 

Til 

1 

cor 

SB01? 

PE1SHOI? 

DELEU 

rum 

1 

100 

11.338955 

0.1807235 

10.4537500 

I. 5167119 

32.1304086 

a 

35 

11.705217 

11.1143088 

35.8924474 

2.3605612 

42.7337873 

3 

as 

39  015318 

13.0804833 

48.3873034 

5.6618727 

71.7832168 

i 

100 

2!  078919 

13.2378304 

10.5461367 

7.9651888 

33.1103380 

a 

33 

27.768673 

10.8305303 

33.4704148 

6.6527798 

445 134778 

3 

as 

10.39)128 

13.0800814 

40.0083316 

17.6522118 

80.0087260 

1 

100 

25.6008)7 

13.8145103 

18.0083042 

16.3633716 

33  6238636 

a 

33 

28.857509 

10.3385737 

33.3730704 

12.1653655 

48  6004300 

3 

as 

15.101632 

36  5508883 

43  4083672 

W  1983616 

17.1408202 

1 

100 

37  588998 

14. 33751  la 

10.3530813 

22.7961016 

37  3645034 

a 

33 

1710605! 

31.7)48176 

25.0480082 

26  6921221 

52  2300863 

3 

13 

71.576865 

35.3138400 

38.0040860 

32.5626985 

00.6351032 

10 

1 

100 

86.692220 

11.6141017 

37  1703681 

36.7989257 

30.7697400 

10 

a 

53 

91.713679 

33.0640033 

31.0718451 

32.6661766 

34  4846360 

!0 

3 

13 

113.388885 

31.3716311 

30  1230830 

19.2965970 

00  6324360 
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Si s  i6:2i  «:k-:.ay  Mil  a 

Al  ALTS IS  OF  VARIAICE  PROCEDURE 

ICilS 


TBO 

LOT 

TYPE 

1 

COST 

SH09T 

PttSBOBT 

OtLEU 

PERERR 

1 

0 

23 

41.371047 

11.5563140 

31.7391566 

3. 5406620 

43.5479783 

1 

1 

13 

33.007586 

11.6409322 

20.3488423 

2.4607285 

35.8161662 

2 

0 

23 

43.744716 

9.4864454 

25.026942? 

2.5149629 

43.0366178 

2 

1 

13 

47.411202 

5.9895602 

15.5695525 

-0.3465S82 

34  7707285 

3 

0 

23 

43.785632 

12.4505865 

29.5076300 

4.2145021 

4:. 5141885 

3 

1 

23 

37.872700 

1C. 3566473 

16.9065600 

2.1020174 

35.6703577 

4 

0 

23 

30.524310 

10.6157763 

3.  ^050638 

0.8062994 

42  0C97394 

4 

1 

13 

39  700390 

11.0676942 

22.1366966 

2  4445151 

35.52C3769 

2 

5 

0 

23 

43.350622 

0.6342623 

26.737618? 

2.0721883 

43  0999685 

ft 

3 

13 

y  HSJI1 

7.9292972 

16.2376048 

1 .5903031 

35  4964677 

1 

0 

S3 

30.031055 

13.0669956 

2*  6954937 

8.7516646 

46.4335355 

1 

1 

13 

28.874418 

12.1004129 

16.5875848 

11.5977498 

39.0430715 

2 

0 

23 

30.398780 

17.6419660 

31.7427596 

9.6020431 

46.9429139 

2 

1 

13 

27.042955 

14  9754911 

22.8239736 

8.4934780 

27.6822477 

3 

0 

23 

26.257756 

17.2402766 

35.3206377 

7.3350028 

45.2473374 

3 

1 

13 

24.62747! 

20.4172006 

27.4145939 

JO  4260946 

30.2119962 

4 

0 

23 

24.279529 

20. 8913361 

27.8439374 

fi. 8573622 

4!  7!’94B3 

4 

1 

13 

21.041015 

21. 0609856 

31.7423638 

12.5639817 

39.1862800 

5 

0 

S3 

32  6*4283 

14.3192257 

26  4631727 

5.9595618 

44.760536! 

5 

1 

13 

31.339166 

14  4290809 

21.3669000 

7.2184655 

37.3133315 

1 

0 

S3 

34  469644 

ir.tfWSJO? 

20  0851703 

14. 1850922 

40.18*2730 

1 

1 

13 

41.118818 

'  6961928 

13.576768! 

16.9300956 

40.9577392 

2 

c 

22 

23.616962 

22.9899191 

24.6012409 

18.1326599 

51.8502196 

2 

1 

13 

23.250635 

22.1289815 

30.0550482 

14.8521044 

39  7377606 

3 

0 

23 

24.625066 

19  4927530 

28.3521514 

13  9059178 

45.2993648 

3 

1 

13 

24.376688 

17. 1930831 

15.2419442 

10.1690(05 

40.9046746 

4 

0 

23 

20.361613 

26.1091013 

32  9874873 

18.6926481 

50.3356887 

4 

1 

13 

18  435953 

25  0682754 

31.2351270 

24.9667382 

43.5103492 

5 

0 

23 

36.184163 

12  7690915 

22.8720964 

11.0016469 

47.5330683 

3 

1 

13 

52.263881 

0  1694499 

17.0634691 

15.9826520 

39  6738515 

1 

0 

S3 

42.78795? 

15  4632644 

21.1513441 

10.1724037 

49.4641078 

1 

1 

13 

51.806476 

14.6000553 

13.2077504 

30.2347766 

44  227886! 

2 

0 

S3 

23.550608 

24.2077000 

34.9526161 

24.7710207 

52.6125083 

2 

1 

13 

20.044284 

26  9793531 

25  6001403 

34.1821616 

46.0679400 

3 

0 

23 

28  140454 

22.4776647 

25  9634655 

23.1329235 

52.2972661 

3 

1 

13 

27.566501 

19  111721! 

J 4. 6085962 

22. 7646106 

42.1663369 

4 

0 

S3 

16.958640 

29.2236013 

32.4408209 

22.3015001 

52  7158009 

4 

1 

13 

17.548620 

30.9676431 

30.5859012 

31.9703498 

44.6260792 

3 

0 

23 

107.295101 

7.0250400 

21.3764466 

29.5371192 

55  4533948 

3 

1 

13 

116  439321 

3.4934194 

10.1796363 

24.0117191 

42.9458046 

JO 

1 

0 

S3 

07.733916 

16.2763010 

28. 1)56917 

30.6372496 

54  9987726 

10 

1 

1 

13 

•7.140712 

16  0813172 

34.0607282 

29.1600243 

43  0957077 

10 

2 

0 

23 

18.876696 

34.1621043 

40.0521393 

26.9280077 

52.0201922 

10 

2 

1 

13 

16  734202 

34.165633! 

33.2622095 

36.2330777 

46.7296962 

10 

3 

0 

23 

2!  155230 

29.6960513 

32.96’:*99 

17.8335566 

48  1473513 

10 

3 

1 

13 

30.357230 

26  3026073 

29  2643617 

33.3266547 

44  5649046 

10 

4 

0 

S3 

1’. 193436 

32.1739126 

45.8405912 

27.8219217 

52.0604026 

10 

4 

1 

13 

14.617915 

39.7102165 

36  4171100 

41.4959911 

47  7*575*? 

!0 

3 

0 

S3 

322.9:50*2 

5  *825259 

i  8925M7 

3?  O925504 

5?  005279? 

10 

3 

1 

13 

20?  ?99-3i 

:. 3076923 

5  4945054 

2?  29?‘,?2fc 

43  61765*2 
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AULTS!?  OF  V61JUCE  PROCELUII 

CHS 


T80 

LOT 

m 

1 

COST 

SMB 

PUSBOBT 

Mlt»t 

mm 

2 

1 

i 

20 

35 .052615 

10. 2024183 

20 . 7079948 

1.1661634 

30.9969655 

2 

1 

2 

n 

41.600691 

14.0012802 

27.3707455 

4.4610651 

43.5607145 

2 

1 

3 

5 

41.557320 

11.6808524 

55.8594920 

6.1401430 

73.6213380 

2 

2 

1 

20 

47.600213 

6.8613858 

It. 7695575 

0.0131525 

30.6992040 

2 

2 

2 

1! 

44.792256 

1.8682455 

22.215V76 

2.1314600 

42.596660C 

2 

2 

3 

5 

35 . 550996 

13.1125118 

39.6584320 

5.865664} 

71.8584680 

2 

3 

t 

20 

41.004741 

1C  6119542 

3.2646666 

31.6641095 

2 

3 

2 

11 

41.565353 

1:  8368606 

24.7918253 

2.5342203 

42.8284255 

2 

3 

3 

5 

44.326674 

15.4161560 

46.3435940 

6.2111646 

72.0292140 

2 

4 

) 

20 

41.729544 

10.2501381 

21.4077739 

1.3334266 

31.0405510 

4 

2 

11 

40.166849 

11.5014361 

29.5219966 

0.5150652 

41.9735809 

4 

3 

5 

29.699238 

12.2088160 

52.6196340 

3.4160514 

69.6292560 

5 

1 

20 

40.707662 

1.1816863 

17.3804241 

1.1650656 

31.2511130 

5 

2 

11 

40.176935 

10.1860665 

25.5608918 

2  0146805 

42.6674564 

5 

3 

5 

44.090362 

12.2612614 

47.2551600 

5.6015044 

71  6776060 

1 

1 

20 

26.336100 

6.0118118 

12.1702713 

6  4642635 

33.6553610 

1 

2 

11 

27.966749 

18.8340863 

25.5925355 

8.6225463 

45.5071055 

1 

3 

5 

46.217230 

11  1035550 

61.5422760 

18.1651506 

80.3691780 

2 

1 

20 

26.156653 

13.5316664 

20.2139254 

6.6514014 

32.7959850 

2 

2 

11 

29.261732 

18.8101135 

35.9806564 

6.1661523 

44.5734627 

2 

3 

5 

4MJ5640 

24.5550540 

45.J4544M 

24.8141004 

84.5656900 

3 

1 

20 

23.236006 

16.2868216 

24.3073828 

1.6612165 

33.0090885 

3 

2 

11 

25.979620 

16.1165663 

40.8742818 

6.6660612 

44.1567055 

3 

3 

5 

34.717476 

21.3318380 

46.6008460 

14  6343534 

78.3034360 

4 

1 

20 

20.660239 

11.2168133 

25.6360427 

10.030110! 

33.6315635 

4 

2 

11 

22  605149 

26.1005810 

28.5626696 

1.5211353 

44.6306436 

4 

3 

3 

33.940186 

24.5616660 

45.2302140 

16  1866040 

7e.  8624200 

5 

1 

20 

26.946743 

11.0346836 

15.4130612 

5.1636160 

32.2596960 

5 

2 

11 

33.000916 

11.2356815 

31.3417306 

5.2646331 

43.6924718 

5 

3 

3 

45.946590 

21.2461511 

46.7724620 

13. 1766636 

77.8429060 

1 

1 

20 

31.965145 

10.0038611 

12.7921809 

11  4766021 

36. 1300145 

1 

2 

11 

37.846075 

I1.3IT1226 

20.5436114 

12.0525151 

46  2744445 

1 

3 

3 

54.287344 

36.4161800 

31.3260520 

16.0366076 

81.7711420 

2 

1 

20 

20.555307 

15.5066313 

22.5390190 

15.3736412 

35  5405940 

2 

2 

11 

21.105462 

26.4164208 

24.6995645 

10.6154646 

46  2545855 

2 

3 

3 

39  625566 

42.8865040 

46.6137)00 

41.3686402 

97.9527780 

3 

1 

20 

21.721764 

13.6811285 

19.2675513 

16.1133133 

36.0313805 

3 

2 

11 

23.282272 

16.1661116 

19.6243564 

10.3066660 

45.7837164 

3 

3 

3 

36.390740 

36.6253800 

46  3651626 

21.6760430 

83.1629340 

4 

1 

20 

17.067665 

21.5)36860 

27  7856258 

20.6123110 

37.3547355 

4 

2 

11 

16  663296 

26.0303016 

36  3856848 

17.3050606 

46.7086627 

4 

3 

3 

32.222102 

60.3114260 

11.8427620 

20.1375220 

66.0940420 

3 

1 

20 

36  673963 

1.1461111 

12.6265942 

11.4373305 

34.0725945 

3 

2 

11 

46.690417 

13.1062515 

25  4218347 

11.4001544 

46.4736409 

5 

3 

3 

62.116390 

23  6630264 

43.1342500 

21.1427754 

84.0472520 

1 

1 

20 

43.264639 

7.8271643 

6  8315352 

30.0047000 

36  5985290 

1 

2 

11 

36.396169 

21.8303420 

21.7263349 

26.3062070 

53.0412682 

1 

3 

3 

73.912066 

26.6164100 

44.5122660 

11  3543004 

79.4425000 

2 

1 

20 

18.603469 

21.0664645 

28.6597792 

26.0531210 

36.3990800 

2 

2 

n 

21  265656 

26.1804755 

35  4075536 

30  4010403 

54.1605227 

2 

3 

3 

36  406094 

45.1031740 

35.3266640 

31.3315136 

69.0447120 

3 

i 

20 

23  895520 

14  6323610 

14.2323318 

30.1253005 

36.7727260 
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SIS  16:3;  ttMESDAY.  mil  19.  1969 

HALTS  IS  OF  VtilAKCE  PBOCEDORE 


MEilS 


TBO 

LOT 

TU 

1 

COST 

sbok 

PEBSHOBT 

DELES* 

PEBEfifi 

6 

3 

2 

11 

27.546029 

20.2852078 

26.2901736 

23.0644005 

50.4077764 

e 

3 

3 

5 

44.935882 

48.6303040 

39.1168420 

34.3645163 

92.2171600 

8 

4 

1 

20 

15.895039 

24.2880340 

32.8275536 

24.5120144 

36.0240930 

8 

4 

2 

11 

17.873941 

34.6861464 

27.1692995 

22.1681112 

50.3316236 

8 

4 

3 

5 

29.933370 

41.4523800 

37.6944460 

36.8919077 

96.2197460 

e 

5 

1 

20 

65.966324 

4.4138723 

12.2186063 

22.6327446 

36.0285490 

6 

5 

*> 

« 

11 

130.426257 

5.4738664 

19.1371193 

30.4527605 

52.25*1516 

6 

5 

3 

5 

170. 695912 

11. 6885564 

33.6226220 

41.7741663 

97.2513780 

10 

1 

1 

20 

47.420356 

15.8641267 

25.9562476 

28.9722724 

36.7595255 

10 

1 

2 

11 

86.360573 

18.0721481 

36.1151327 

34.4916344 

56.0:22016 

10 

1 

3 

5 

154.067180 

12.3558466 

39.2117460 

25.0267266 

86. 7692560 

10 

2 

1 

20 

16.191307 

32.4835060 

33.5494623 

31.2422929 

39.6249215 

10 

2 

2 

11 

17.  429388 

28.3838382 

42.5574600 

36.6349952 

57.7231982 

10 

2 

3 

5 

27.241050 

51.4265720 

42.6972440 

12.6594764 

75.2993720 

10 

3 

» 

20 

22.766305 

26.7063487 

31.0177101 

25.2673494 

37.6059570 

10 

3 

2 

11 

23.412759 

30.3463181 

27.0636235 

26.6569306 

53.2452155 

10 

3 

3 

5 

52.069566 

37.5258740 

44.2174960 

4. 044C 172 

66.0354740 

10 

4 

1 

20 

13.305022 

32.5238085 

39.6820716 

36.6410049 

42.2555660 

10 

4 

2 

11 

16.589749 

30.8080800 

42.0462555 

26.0107093 

52.5937991 

10 

4 

3 

5 

27.378846 

53.3333300 

53.5211560 

19  6628360 

79.0061200 

10 

5 

1 

20 

333  778110 

0.3660666 

5.4464285 

29.0517086 

40.4027100 

10 

5 

2 

11 

312.785927 

0.3030303 

7.5757573 

30.1346145 

52.6467700 

10 

5 

3 

5 

306.077680 

1.7142858 

15.7777774 

34.6773286 

93.1520580 

LOT 

m 

TYPE 

1 

COST 

SHOW 

PEBSHOBT 

DELEBB 

PEBE1B 

l 

0 

60 

30.653670 

1 l . 3792636 

16.2729969 

10.8061813 

34.4807863 

i 

1 

40 

46.069923 

8.7288388 

16.3196196 

22.8653331 

36.3440132 

2 

0 

30 

42.525395 

15.8528480 

26.6819516 

22.4171087 

50.0428373 

2 

1 

25 

52.102126 

18.7556800 

25.7750685 

11.7851753 

47.4627184 

3 

0 

25 

74.004228 

21.8786230 

46.4903644 

17.3487408 

80.3946828 

1 

0 

60 

23.700616 

17.8461584 

25.2176700 

9.9437368 

34.0533748 

1 

1 

40 

29.103552 

18  4857060 

23.0394267 

24.8576000 

37.4498300 

2 

0 

30 

29.672711 

18.8328832 

34.4473138 

24.2628534 

51.0519603 

2 

1 

25 

23.365622 

25.8512846 

29.4425968 

9.0437718 

46.6742384 

3 

0 

25 

36.431674 

35.5581764 

42.0082980 

23.3074612 

63.7642036 

1 

0 

60 

23.266741 

16.2201810 

22.8273606 

9.5955491 

34.1374246 

1 

1 

40 

31.412057 

16.833400! 

20.6728918 

22.2457861 

36.4354935 

2 

0 

30 

31.116666 

17.4540810 

32.6968750 

18.2425698 

47.9005607 

2 

1 

25 

25.037116 

22.5060631 

21.8552247 

9.5471515 

46.5457924 

3 

0 

25 

42.928468 

33.5604684 

45.1291880 

18.1877488 

78.0836476 

1 

0 

60 

21.357070 

18.1881218 

28.5435471 

11.6339671 

34.7488730 

1 

1 

40 

22.303260 

24.1040244 

30.9442133 

30.2632577 

39.1263550 

2 

0 

30 

24.000355 

24.3814473 

35.3552429 

18.9793152 

46.1339550 

2 

1 

25 

22.213608 

26.4850217 

19.5955072 

10.5664124 

47.0641104 

3 

0 

25 

30.634749 

34. 8257400 

46.2216824 

21.5156456 

62.3623172 

1 

0 

60 

106.584046 

7.2205870 

13.0662691 

10.4572035 

34.5580908 

1 

1 

40 

100.160343 

4.5450508 

11.9441535 

19.7695765 

36.1701950 

2 

0 

30 

103.134476 

7.8458181 

24.1243104 

22.9430425 

49.4906550 

2 

1 

25 

123.994909 

11.1883838 

19.0272545 

7.4190640 

45.7682360 

3 

0 

25 

125.786387 

14.1824164 

37.3524583 

23.2351481 

84.7942800 
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SIS  16:3!  KDIESDIY  UIIS  It,  1988  32 

UOITSIS  OF  V4BIAICE  PSOCEDOEE 

KUS 


TBO 

LOT 

?ii 

TTPE 

t 

COST 

SBOIT 

PESEKIST 

D£L£B£ 

nuit 

2 

1 

1 

0 

12 

JO. 173050 

11.5235074 

23.1087583 

0.4345740 

31.1560492 

2 

1 

1 

1 

8 

30.171040 

0.2222072 

18.0718525 

2.200)04 

30.7410650 

3 

) 

2 

0 

0 

44.013005 

11.5300030 

20.7108703 

5. 0)9403 

43.2489033 

2 

1 

2 

1 

5 

30.140040 

17.1107040 

25.7520280 

2. ’002104 

43.0:46680 

2 

1 

3 

0 

5 

41.557330 

11.0000524 

55.8504020 

6. 140^430 

73.62:3360 

2 

2 

1 

0 

12 

45.01744! 

7.0473557 

18.3284207 

0. Ofc 14300 

3!. 1506575 

2 

2 

1 

1 

8 

50.574305 

5.5574550 

14.4345127 

*0.0642700 

30.0220237 

2 

2 

2 

0 

8 

40.0:7357 

0.0420470 

26.2427457 

4. 5720246 

42.7903300 

2 

2 

2 

1 

5 

43.350140 

0.0000272 

17.3856160 

-0.7082316 

42.3666560 

2 

2 

3 

0 

5 

35.550000 

13.7125770 

30.6504320 

•  8850042 

71.6564660 

2 

3 

1 

0 

12 

44.330731 

11.5204758 

23.5177703 

3.7707050 

32.2187256 

2 

3 

1 

1 

8 

30.100750 

0.3001900 

17.2307052 

2.407*130 

30.8321850 

2 

3 

2 

0 

6 

43.300433 

11.7770050 

27.4573633 

3.4211647 

42.3425650 

2 

3 

2 

1 

5 

40.500000 

11.0085000 

21.5031706 

1.4006630 

43.4114340 

2 

3 

3 

0 

5 

44.330074 

15  4707500 

46.3435040 

6.2171040 

72.0292140 

2 

4 

1 

0 

12 

43.750407 

0.7107544 

21.0538650 

-0.4104007 

30.9479275 

2 

4 

1 

1 

8 

40.180114 

11.0457177 

20.5886374 

3.0550733 

31.1797362 

2 

4 

2 

0 

6 

41  343530 

11.0441537 

33.6110002 

1.0320680 

41.3470633 

2 

4 

2 

1 

5 

30.033333 

11.1033704 

24.6141172 

0.0266616 

42.7254020 

2 

4 

3 

0 

3 

30.000330 

13.2008700 

52.8108340 

3.4700574 

69.6262580 

5 

1 

0 

12 

44.541005 

8.0050250 

18.5775025 

J. 0525240 

31.7013000 

5 

1 

1 

8 

34 .050350 

0  1207070 

15.5846715 

1.7009524 

30. 5756325 

5 

2 

0 

6 

40.351005 

0.8053127 

25.0507200 

2.0370775 

42.0624333 

3 

2 

1 

5 

30.071030 

10.0134414 

25.0822080 

1.4004042 

43.3694640 

3 

3 

0 

3 

44.000303 

12.2072014 

47,2551600 

5. 60 J 5044 

71.6776060 

1 

1 

0 

12 

33.075303 

10.1302041 

13  0344022 

7.1707271 

33.3093617 

1 

1 

1 

8 

31.503300 

4.0500035 

10.0730550 

12.0045082 

34.1743600 

1 

2 

0 

6 

30.857533 

15.1020532 

21 . 1451002 

7.8856345 

44.4021650 

1 

2 

1 

5 

34.541010 

23.3314000 

30.0203404 

6.5066404 

46.6330100 

1 

3 

0 

5 

40.317230 

17.7033550 

61.5422700 

18.1051506 

60.3661760 

2 

1 

0 

12 

23.040200 

10.5120403 

23.0706167 

4.106)227 

32.3324983 

2 

1 

1 

8 

20.024220 

0.0747335 

14.7108006 

11.1263165 

33.4912150 

2 

2 

0 

3 

34.332412 

14.1301040 

36.1304763 

7.7420005 

44.7260663 

2 

2 

1 

3 

22.012010 

24.4100000 

35.7005100 

4.2600316 

44.3670000 

2 

3 

0 

3 

41.133040 

24.3530540 

45.3454400 

24.0747004 

14.6656600 

3 

1 

0 

12 

31.010310 

14.2033440 

24  5830455 

5.3400074 

32.7683523 

3 

1 

1 

8 

25.302255 

10.3072800 

23.0025300 

11.6682850 

33.3701012 

3 

2 

0 

0 

20.000402 

14.7104007 

47.3046150 

5.2200015 

42.6505567 

3 

2 

1 

3 

23.431020 

22.1030014 

33.0406820 

8.3605600 

45.0566640 

3 

3 

0 

3 

34  717470 

17.3370580 

46.6000460 

14  0343534 

78.3034360 

4 

1 

0 

12 

20.404230 

10.8727007 

21.0066004 

6.7415143 

33.2337017 

4 

1 

] 

1 

21.004245 

17.7302034 

32.4002002 

14.0045053 

34.7282662 

4 

2 

0 

6 

23.070500 

15  0000045 

26.0100477 

6  5104306 

43.2162617 

4 

2 

1 

5 

20.055840 

20.3772000 

30.5450100 

8.7231000 

46.3243020 

4 

3 

0 

3 

33.040100 

24.5010000 

45.2302140 

16.7609040 

76.6624200 

5 

1 

0 

12 

20.500023 

11.0243820 

13.0477720 

3.2500157 

32.0667506 

3 

1 

] 

1 

32  407422 

11.1001002 

17.0100050 

6.0044730 

32.5461137 

3 

2 

0 

0 

33.723303 

15  1330720 

34.0402105 

4.1416673 

42.6521317 

3 

2 

1 

5 

32.133050 

10.7533540 

27.3703400 

5.6166521 

44.0406100 

5 

3 

0 

3 

45  040500 

11. 34875)0 

46.7724020 

13.7700658 

77.0420060 

) 

I 

0 

12 

23.317744 

13  0402050 

12  7107373 

13.576540) 

35.5016602 

1 

1 

l 

0 

44  000240 

4.2301100 

12.0023462 

23.3316616 

37.0475325 
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SIS  16:21  HEDfESDiY  IPEIL  10.  108i  22 

A1AITSIS  OF  VUWCE  PBDCFDUSI 


HE  IIS 


TYPE 

1 

COST 

shost 

PEISBOBT 

0EL£tLB 

PEtf&B 

0 

6 

38.875202 

0.7222217 

25.4506350 

12.190)740 

45.4014233 

1 

5 

35.410932 

13.2323238 

14.6558430 

11.8873244 

47.2140700 

0 

5 

54.267344 

36.476)900 

31.3260520 

18.0346976 

61.7711420 

0 

12 

17. 177731 

14.6900733 

19.28044:7 

10. 5042260 

34 . 4303792 

1 

e 

25.621642 

16.7234837 

27.4268850 

22.4262640 

37.2050162 

0 

6 

23.762530 

22.6956233 

16.7324483 

17.8465332 

48.3004367 

1 

5 

19.457024 

30.7777780 

34.2601040 

2  73C2490 

42 . 786**640 

0 

5 

39.825586 

42.9965040 

46.8137100 

4 : . 3668402 

07.9527760 

0 

12 

17.917209 

13.7037033 

23.1301968 

12.0629000 

35.0275702 

J 

8 

27.428461 

14.1687612 

13.4735831 

23.6590832 

37.5370625 

0 

6 

26.402622 

14.9603167 

21.2852163 

11.0751658 

45.1861283 

1 

5 

19.493852 

21.9999980 

18.07:3220 

9.3849722 

46.5008220 

0 

5 

38.590740 

36.8253960 

49.3651620 

21 .6760630 

83.8620340 

0 

12 

16.128342 

21.9675917 

31.3270012 

15.6526455 

36.1663683 

1 

6 

18.477109 

20.8333325 

22.4222)26 

26.60)8259 

39.1372862 

0 

e 

18.944413 

22.6851833 

26.9272636 

16.4019243 

47.2006400 

1 

5 

18.369960 

34.4444440 

45.3357000 

16.8305980 

50.5072500 

0 

5 

32.222102 

40.5714260 

41.8427620 

28.7375228 

86.0040420 

0 

12 

25.163864 

9.6754082 

17.1746405 

6. 3680503 

33.1848642 

1 

8 

53.939162 

2.9071070 

5.6095247 

10.0417732 

35.4041600 

0 

6 

44.279573 

0.0939157 

17.3816603 

11.6176097 

45.6010267 

1 

5 

49.583430 

10.2410546 

35.0697800 

i;  0880600 

4\  2856180 

0 

5 

62.118390 

23.8839264 

43.1342500 

21.1427754 

64.0472520 

0 

12 

32.805423 

7.0675424 

10.1267763 

11.2624813 

34.6125233 

1 

8 

59.003462 

0.1866471 

0.3866736 

34  7660495 

30.2775375 

0 

6 

36.816265 

2C. 3703700 

23.7330530 

32.6739293 

53  7852833 

1 

5 

40.296498 

23.8023084 

10.3182732 

22.0615400 

52.1484500 

0 

5 

73.912066 

29.6764706 

44.5122560 

17  3543864 

70.4425000 

0 

12 

16.826624 

18.6507933 

30.9505942 

15.6702480 

35.6612433 

1 

8 

21.768736 

24.7700212 

25.2235567 

41.6339475 

42.4758350 

0 

6 

24.610770 

18.2530667 

42  6449533 

37.3413203 

56.1148683 

1 

5 

17.265160 

35.7142860 

26.7226740 

22.2610042 

51.6153080 

0 

5 

36.406004 

45. J031740 

35.3266640 

31.5265156 

69.0447120 

0 

12 

17.878703 

14.6329367 

14.57)9613 

14.1325960 

35.4636900 

1 

8 

32.920745 

15.4315475 

13.7228876 

29.1)43764 

36.7362600 

0 

0 

34.667767 

16.3690478 

34.4103267 

31.7739158 

52.6979000 

1 

5 

18.999944 

25.0000000 

16.5450900 

12.6560821 

47.6596280 

0 

5 

44.935682 

48.6363640 

39.1166420 

34.3645163 

92.2171600 

0 

12 

14.974919 

20.4661)90 

28  9520525 

11.5336307 

34.8376275 

1 

8 

17.275219 

30.0099200 

38.6408052 

43.0705900 

42.6037912 

0 

6 

17.700507 

36.5079350 

35.0403367 

30.01)8992 

52.2186600 

1 

5 

17.986062 

32.5000000 

17.7240548 

12.7555656 

48.0669400 

0 

5 

20.033370 

41.4523800 

37.8044460 

38.6919077 

96.2197460 

0 

12 

78  323194 

5.9603004 

12.5063405 

16.0564024 

36.3837067 

i 

0 

07.431020 

2.0942301 

11.7870050 

32.40725T0 

40.4056125 

0 

0 

112.404303 

5  2586387 

28  7448467 

47.1343466 

56.7611183 

1 

5 

152.052602 

5.7321196 

7.6078464 

10.4348570 

46.6657920 

0 

5 

170.696912 

11.6998564 

33.8226220 

40.7741663 

07.2513760 

0 

12 

35.096220 

15.2777783 

22.2862325 

21.5845751 

37.6134283 

l 

8 

65.906560 

17.0436492 

31  4612700 

40.0538185 

40.4788712 

0 

0 

61.064922 

21.5476183 

34.36)1933 

53.4)63665 

63.2663017 

1 

5 

121.115354 

16.1015860 

38.2196600 

11  7819535 

47.2831740 
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SAS  ii:l-  KL'KESWY  APRIL  it.  :988  34 

AM ALTS  IS  OF  VARIARCE  PROCEDURE 

lEAVS 


TBO 

LOT 

VAR 

TYPE 

If 

COST 

SHORT 

PERSHORT 

DELEiR 

llslii 

10 

1 

3 

0 

5 

154.067100 

12.3559466 

39.2117460 

25.0287266 

6( .7692560 

10 

2 

1 

0 

12 

15.032989 

29.9206342 

33.6542769 

19.2946555 

36.6720958 

10 

2 

1 

1 

6 

17.928786 

36.3526132 

33.3922904 

49.1637475 

44.0541600 

10 

2 

2 

0 

6 

19.601488 

26.3333333 

50. 4709433 

53.9106295 

63.3:70683 

10 

2 

2 

1 

5 

14.822866 

30.6666640 

33.0540800 

20.7440059 

51.0:05540 

10 

2 

3 

0 

5 

27.241050 

5i.4285720 

42.8972440 

12.8594764 

75.2993720 

10 

3 

1 

0 

12 

14.496465 

26  9444450 

28.3330208 

12.63S9357 

35.2087758 

10 

3 

l 

1 

b 

35.166067 

26.3492059 

35. 044*440 

44.2599700 

41.7317287 

10 

3 

2 

0 

6 

24  040148 

29.4444450 

32.9368515 

39.7135813 

56.6177233 

10 

3 

2 

1 

5 

22.659692 

31.4285680 

20.0157500 

15.6333502 

48.1961940 

10 

3 

3 

0 

5 

52.069566 

37.5259740 

44.2174960 

4.0446172 

68.0254T4C 

10 

4 

1 

0 

12 

12.519354 

26.9444450 

39.4173164 

24.6585457 

36.5641800 

10 

4 

1 

1 

8 

M. 463524 

40.8928562 

40.5792050 

59.6146937 

47.7926950 

10 

4 

2 

0 

6 

18.053757 

25.0000000 

52.2866700 

40.9312446 

56.6747500 

10 

4 

2 

1 

5 

14.832646 

37.9999980 

29.7577580 

12.5060667 

47.6966580 

10 

4 

3 

0 

5 

27.378846 

53.3333300 

53.5211560 

19.6828350 

79.0061200 

10 

5 

1 

0 

12 

368.291683 

0.2777777 

3.1250000 

24.7492242 

39.4518125 

10 

5 

1 

1 

8 

282. 00*750 

0.5000000 

6.9265712 

37.5054352 

41.8290562 

10 

! 

2 

0 

6 

264.813433 

C. 5555555 

13.8888883 

46.2642210 

58.1565650 

10 

5 

2 

1 

5 

346.232920 

0.0000000 

0.0000000 

6.3550668 

46  4794160 

10 

j 

3 

0 

5 

306.077680 

i  '142656 

25.7777774 

34.8773286 

93.1520560 

APPENDIX  E 


ANOVA  RESULTS  —  GROUP  II  DATA 


The  following  text  provides  the  code  used  in  the  SAS 
routine.  Output  consists  of  all  subsequent  pages. 


DATA; 

INPUT  SET  TYPE  VAR  COFVAR  TBO  LOT  COST  SHORT 
PERSHORT  DELERR  PERERR  BIAS; 

DROP  TYPE  VAR  COFVAR  BIAS; 

CARDS; 

PROC  ANOVA; 

CLASS  SET  TBO  LOT; 

MODEL  COST  SHORT  PERSHORT  DELERR  PERERR  =  LOT  TBO 
SET  LOT* TBO  LOT* SET  TBO* SET; 

MEANS  LOT  /  TUKEY  E=LOT*SET; 

MEANS  TBO  /  TUKEY  E=TBO*SET; 

MEANS  SET  LOT* TBO  LOT* SET  TBO*SET  /  TUKEY; 

TEST  H=LOT  E=LOT*SET; 

TEST  H=TBO  E=TBO*SET; 

OUTPUT  OUT=PLOTDATA  P=YPRED  R=YRESID; 

PROC  UNIVARIATE  NORMAL  PLOT; 

VAR  YRESID; 

PROC  PLOT; 

PLOT  YRESID*YPRED ; 

PLOT  YRESID*LOT ; 

PLOT  YRESID*TBO; 
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SAS 

6EMEBAL  LIME  Alt  MODELS  PHOCEDUfcE 
CLASS  LEVEL  IFGHMATIOI 
CLASS  LEVELS  VALUES 
SET  5  1  2  3  4  5 

TBO  5  1  4  6  6  10 

LOT  5  1  2  3  4  5 

vmn  OF  uLSEKVATlOKS  II!  DATA  SET  *  !25 
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12:;*  ?SESL‘AT  MAY  i.  1 : 


GEIEBAL  LI  IE  At  MODELS  PBOCEDUBE 


DEPQDEiT  VABIABLI:  COST 


SOUBCE 

DF 

SOM  Or  SOUABES 

£11  SOUlil 

F  VALUE 

PB  -•  F 

B-Stl’ABE 

c  v. 

IDDEL 

80 

565030.40037688 

6411.17333961 

20.27 

0.0001 

0.950001 

41.4280 

EUOB 

64 

29737.76416453 

464. 65256561 

BOOT  MS* 

COST  MEAD 

COBBECTED  TOTH 

124 

584766.16454142 

55580119 

50.80425250 

SOUBCE 

DF 

TYPE  1  SS 

F  VALUE 

PE  >  F 

DF 

type  :::  ss 

F  VALUE 

VC  .  f 

LOT 

4 

170805.61714428 

91.90 

0.0001 

4 

170805. ei714426 

91.90 

0.0001 

TBO 

4 

91515.99144414 

49.24 

0.0001 

4 

9:515.95144414 

49.24 

0.0001 

SET 

4 

6459  4927226! 

3.48 

0.0124 

4 

6459.49272261 

2.46 

0.0124 

T8O*L0T 

16 

246895.12529861 

33.21 

0.0001 

16 

246895.12529861 

33.21 

0.0001 

SET*LOT 

16 

27499.52674268 

3.70 

0.0001 

16 

27499.52674266 

3.10 

0.0001 

SET* TBO 

16 

21853.44702451 

2.94 

0.0011 

16 

21653. 447C2457 

2.94 

0.0011 
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SIS  13:11  TUESDAY.  DAT  5.  198il  3 

QEIEEAL  LIlEil  HDEli  PBOCEt'JSE 

DE?EIDEI7  MI1ASLE:  SHOE? 

SOCBCE  Of  SUE  Of  SOCAEES  «il  S03ASS 

»sa  00  181DT. 45899103  189.95784985 

UllOl  84  5413.31004373  14.58313588 

CO EJECTED  TOTAL  124  21810.77903478 


Of  TYPE  II!  SS  f  YAL3E  fE  )  f 

4  4589.73863595  13.57  0.0001 
4  3084.87142157  9.12  0.0001 
4  380.04155787  1.12  0.3534 
16  4442.28461492  3.28  0.0004 
16  2040.80953508  1.51  0.1247 
18  >659.71222483  1.23  0.2738 


soubce 

OF 

TTPE  I  SS 

f  VALUE 

fE  >  f 

UJT 

4 

4589. 73863595 

13.57 

0.0001 

TBO 

4 

3084.87242157 

9.12 

0.0001 

SET 

4 

380.04155767 

1.12 

0.3534 

TBQ*IQT 

16 

4442.28461492 

3.28 

0.0004 

smor 

16 

2040.80953608 

1.51 

0.1247 

srr»T8G 

16 

1659.71222483 

1.23 

0.2738 

F  VALUE 

pe  >  r 

1-SQOABE 

C.V. 

3.19 

0.0001 

0.749508 

56.2294 

IOC?  USE 

SBOBT  CAM 

9.19690849 

16.35604909 

13:27  TUESDAY .  IUT  2.  1909  4 


S  AS 

geieba:  liras  kkls  fbocedtjbi 


DEPEIDEIT  VAB1ABLE:  PEBSBQBT 


SOOBCE 

DF 

sm  of  sqoues 

Ell  S90UI 

F  VALUE 

PB  >  F 

B-SOUABE 

C.V. 

DDEL 

eo 

36667.66279924 

611.46104665 

2.30 

0.0006 

0.663531 

59.3199 

EBBOB 

04 

16966.06941667 

265 . 40760026 

BOOT  16E 

FEES HOST  ME U 

CORRECTED  TOTAL 

124 

53673.74921561 

16.29133513 

27.46351057 

SOUBCE 

DF 

TYPE  !  SS 

F  VALUE 

H  '  ? 

It 

TTFE  Ji:  SS 

?  VALUE 

rE  F 

LOT 

4 

4415.65662667 

4.16 

0.0047 

4 

4415.65662667 

4.16 

0.0047 

TBO 

4 

34 17. 49657626 

3.22 

0.0160 

4 

3417  49651626 

2.22 

O.OibU 

SET 

4 

4955.99207772 

4.67 

0  0023 

4 

4955.99237772 

4.6? 

0.0023 

TBO*LOT 

16 

4391.22395075 

1.03 

0.4349 

16 

439:. 22395675 

1.03 

0.4549 

SET*LOT 

16 

9261.27023573 

2.10 

0.0146 

16 

9261.27623573 

2.16 

0.0146 

SET*TBO 

16 

10245.81032408 

2.41 

0.0007 

10 

10245.61932496 

2.41 

0.9067 
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Si $  13:27  TUESDAY  MAY  2.  1989  5 

QEIERA1  UKA*  MODELS  PtOCEDURE 


DEPETOEP  VAfilABLE:  Dtim 


SOUBCE 

DF 

SOM  OF  SQUARES 

MEA1I  SQUARE 

F  VALUE 

PR  >  F 

B-SOUAEE 

C.V 

«ML 

60 

45730.57998833 

762.17633314 

3.04 

0.0001 

0.740138 

61.7653 

asoi 

84 

16056.00994200 

250.07515534 

SOOT  NSE 

DELE  EE  HEAR 

COSEECTE0  TOTiL 

124 

61786.58993033 

15.83903897 

25.64392810 

S0D3CE 

DF 

TYPE  :  SS 

F  VALUE 

PE  /  F 

DF 

TYPE  HI  SS 

F  VALUE 

PI  F 

LOT 

4 

1238.06194197 

1.23 

0.3055 

4 

1236.06194197 

1.23 

0.3055 

TBO 

4 

17706.37837516 

17.64 

0.0001 

4 

17706. 3762^5 19 

17.64 

0.0001 

SET 

4 

12737.24578276 

12.69 

0.0001 

4 

12737.24576276 

12.66 

0.0001 

TBO«LOT 

16 

1345.14003183 

0.34 

0.9911 

16 

1345.14003183 

0.34 

C.6911 

SET«LOT 

16 

7350.04017087 

1.63 

0.0459 

16 

7350.04017087 

1.83 

0.0450 

SETtTBO 

16 

5353.71368571 

1.33 

0.2053 

16 

5353  71368571 

1.33 

0.2053 
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SIS  12:2”  TUESl'AT.  MAT 

sieral  likas  models  procedure 


DEPEIDEIT  VIII ABLE:  PERE1R 


SOURCE 

DF 

SUM  OF  SQUARES 

IEAI  SQUARE 

F  VALUE 

PR  >  F 

R- SQUARE 

MODEL 

60 

13415. 76612311 

222.59646872 

7.63 

0.0001 

c.er*363 

ERROR 

64 

1875.25080257 

25.30079520 

ROOT  MSE 

CORRECTED  TOTAL 

124 

1 529 1.0390 1568 

5.41302085 

SOJJKt 

DF 

TTPE  l  SS 

F  VALUE 

n )  ? 

DF 

tt?i  :::  ss 

F  VALUE 

LOT 

4 

142.06486092 

1.21 

0.3144 

4 

142.06486092 

1.21 

TBO 

4 

1701.34034705 

14.52 

0.0001 

4 

1701.34034705 

14.52 

SET 

4 

9915.78741060 

84.64 

0.0001 

4 

9919.76741060 

84  64 

TBO* LOT 

16 

148.12926649 

0.32 

0.9935 

16 

148.12926649 

0.32 

SET *LOT 

16 

853.60554973 

1.82 

0.0474 

16 

653.60654973 

1.62 

SET’TBO 

16 

650.85968830 

1.39 

0.1763 

16 

650.65968830 

1.39 

158$  t 


C.V. 
14.465ft 
PERTH  ME  AH 
37.415C£5b4 


n .  ? 

0.3144 

0.000] 

0.000] 

0.9935 

0.0474 

0.1763 
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SiS  53:2*  TrE£I)AY  MAY  .  Iwi  7 

GEIERAL  LIIEAE  MODELS  PBOCEDQEE 

TDXEY'S  STULEITIZEL  HAIQE  (HSU:  TEST  FOE  VABIABLE-  COS? 

I0TE:  THIS  TEST  COITBOLS  THE  TT?E  I  EXPEEIKITtflSE  E&BOE  HATE. 

BUT  GEJfElAUY  HAS  A  HI  GHEE  TT?E  II  EHBOE  HATE  THAI  EEGtQ 

ALPHA* 0.05  D?*i C  MS£*1718.72 
CRITICAL  VALUE  OF  STUUETIZEI-  kAKH*4.333 
Kill MIV  SI3IIFICAT  D!FFEEEKE*35.824 

KEAIS  KITH  THE  SAME  LETTES  ARE  KJT  SIQIIFICATL?  DIFFERED. 

TUKEY  GROUP I IQ  MEAJi  I  LOT 


A 

121.43 

25 

5 

B 

54.29 

25 

1 

B 

B 

ft.  02 

25 

4 

B 

B 

25.77 

25 

3 

B 

£ 

25.52 

25 

2 
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SIS  33:2“  Tl'Ei'.il  IIB! 

geeef.a:  lieu  cdees  peocedobe 
tcket  e  sTn>ar:2£t  uke  <Bst>>  test  for  tajubee  sbobt 

KITE:  THIS  TEST  COITION  TE  TTPE  !  EEEBIcmSE  EBEOS  SITE. 

BUT  QEEH1LLY  BAS  t  BIGEE  TEE  II  EEBOB  UTE  THU  EW 

U.PU-0.05  OF- If  EE-121. 551 
CEITICAL  ULffi  OF  STTOEIH2ED  RUGE-4.333 
MIBIIIGM  SIGHIF1C1IT  DIFFEHEBCE-fi  1863 

CUE  WITH  TE  SIC  LETTEi  1EI  HOT  SlGlIFtCUTl!  DIFFEBEIT 

TOE!  9B02F1IG  CU  k  LOT 

t  12.131  15  1 

1 

1  11.110  25  3 

8 

i  10.551  25  4 

1 

B  i  12.523  25  1 

B 

E  6 .300  25  5 
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SAS  !2:2*  TESPAT.  HAE  l.  \ii 5  $ 

GESEEAl  LIIEAE  ICDLLS  PROCEDURE 

TUKEY’S  STUI’EIJT  I  ZED  RAKE  <HSD!  TEST  FOR  VARIABLE:  PERSBOBT 
IOTE;  THIS  TEST  COITBOLS  THE  TYPE  I  EZPEE1MDTTVISE  ERROR  BATE. 

BUT  GEERALLT  HAS  A  HIGHER  TYPE  II  ERROR  RATE  THAM  REG*) 

ALPHAS. 05  DF*16  MSE=576.83 
CRITICAL  VALUE  OF  STUDEriZED  RAMEM.333 
HIIIKM  SIGIIF1CAI?  DIFFER£KE=20.847 

CARS  KTH  TEE  SAC  LETTER  ARE  MT  SIQIIFICATLT  DIFFEREIT. 

TUXET  G30GPIK  ICAR  V  LOT 

A  33.874  25  4 

A 

A  31.02!  25  2 

A 

A  28. BOO  25  3 

A 

A  24.257  25  ! 

A 

A  17.467  25  5 
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SAS  13.2*  riESiAI  HAT  2 

GEIEHAL  LI  tut  DUELS  PBOCEDUBE 

TEXET'S  OTMJTIZEI  BAIGE  (BSD)  TEST  FOB  VAB1ABLE:  OELEBB 
IOTE:  THIS  TEST  COITIOLS  THE  TT?£  !  EXBEBIEmi’SE  EBBOB  BATE. 

BUT  OEIEBALLY  H1S  A  HI  GHEE  TTTE  11  EEBOB  UTE  TEAS  BE3W 

ALPU>0.05  OF*  16  EE-459.318 
CBITICAL  VALUE  OF  STEOEKTIEHS  BAlGE-4.333 

mimn  sioiificait  diffebike-is.sib 

EUS  KITH  THE  SAE  LETTES  ABE  ACT  SI0BIF3CABTLY  DIFFEBEHT. 

tueey  ato'jFiBG  ffiit  #  lot 

A  29.016  25  1 

A 

A  26.264  25  3 

A 

A  26.142  25  2 

A 

A  23.611  2!  4 

A 

A  22.526  23  5 


Mi  10 
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SAS  12:2?  TUESDAY  HAT  1.  :9ii 

SEVERAL  1 1  NEAR  ICDELS  PROCEDURE 

TDKEY'S  STCDEITIZED  RANGE  tHSDs  TEST  FOR  VARIABLE:  PERES* 

NOTE.  THIS  TEST  CONTROLS  TEE  TYPE  I  EXPERIKENTV1SE  ERROR  RATE. 

BUT  GENERALLY  US  A  HIGHER  TYPE  II  ERROR  RATE  THAN  REGNO 

ALPHA*0.05  OF* 16  CE*52.35D4 
CRITICAL  VALUE  OF  STUDENTI2EI)  RAIGE*4.333 
mim  SIGIIFICAV?  DIFFEHICE*£.3282 

HEARS  HITE  TEE  SAC  LETTER  ARE  HOT  SIGNIFICANTLY  DIFFERENT. 

TUXEY  GROUPING  CAf  K  LOT 

A  36.036  25  I 

A 

A  36.166  25  3 

A 

A  3?. 452  25  2 

A 

A  36.474  25  4 

A 

A  36.044  25  5 
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GE1EUL  LllEli  MODELS  FHOCEDHHE 

T’JXET'S  STHDEHTIZir  HAKE  !HSBi  TEST  FOB  VABIA81.E:  COST 
IOT£  THIS  TEST  CorBOLS  THE  TYPE  !  EEFEBIUEHTIISE  EBBOB  RATE. 
BUT  GE1EBELL?  HIS  i  HIGHEB  TYPE  II  EBBO&  BITE  THU  HEWS 

ALPHA-0. 05  OF=it  BEE' 1365  84 
CRITICAL  VALUE  OF  STC0EHTI2EL'  BAJGE-4.JJJ 
HIlIHlm  SIGHIFICAIT  MFFEBE1CE-3!  024 

£US  WITH  THE  SUE  LETTER  ABE  HOT  SIMIF’CUTIT  LIFFEBEIT. 


GROUPING 

IEU 

N 

TBO 

i 

102. 15 

25 

)0 

B 

52.24 

25 

6 

h 

B 

36.05 

25 

2 

B 

B 

32.66 

25 

6 

B 

B 

27.70 

25 

4 
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SAS  12:2*  ?TtS: AY  HA:  1 

QEMEUAL  L1IEAB  MODELS  PROCEDURE 

TOKEY’S  STBDBfTIZED  HUGE  lESDi  TEST  FOE  TAE1ABLE:  SHORT 
MOTE:  THIS  TEST  COITBOLS  THE  TYPE  I  EIPERIHEKTilSE  ERROR  RATE. 

BIT?  GE  HE  RALLY  HAS  A  HIGHER  TYPE  I!  ERROR  RATE  THAR  REGMQ 

AL?BA*0.05  DF*16  USE =103  732 
CRITICAL  VALUE  OF  STUl'EHTIZEi  RAME*4  333 
MIKIHUM  SIGIIFICAJT  RIFF£R£ICE= 6 . B254 

CAMS  HITE  TEE  SAJC  LETTER  ABE  MOT  SIGIIFICAMTLT  LIFFEBEJT 

TEXET  5&0UMHG  HE  AM  M  TBO 

A  22.355  25  10 

A 

A  16.672  25  6 

A 

A  16.927  25  6 

A 

B  A  14.766  25  4 

B 

B  7.659  25  2 
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SiS  13:27  TUESDAY .  MAT  2.  1010  14 

GEIERAL  LIRA*  MODELS  PROCEDURE 

TtJKEY'S  STUDEWTI2XD  RAISE  (BSD)  TEST  FOE  VARIABLE:  PEBSHOBT 
IOTE:  THIS  TEST  COITROLS  THE  TYPE  !  EIPEB1IGITWSE  EEBOS  BATE. 

BUT  SEVERALLY  HAS  A  HIGHER  TYPE  II  ERROR  RATE  THAI’  &£<M 

AL?RA*0.05  DF*lfc  MSE=640.363 
CRITICAL  VALUE  OF  STUDEITIZED  RAISE *4. 3 33 
MIIIMUM  SIGIIFICAI?  DIFFER£ICE*2I.02fi 

MEADS  r  TH  THE  SAME  LETTER  ARE  10?  S1QIIF2CAITLY  LIFFEREIY. 

Tl'KEY  GROUP  IIS  MEAI  I  TBO 

A  35.661  25  10 

A 

A  30.453  25  6 

A 

A  26.205  25  6 

A 

4  23.677  25  2 

A 

A  20.921  25  4 
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S AS  1 3:2*  T3ES5AY.  MAT  2 

GEIEBAL  LIIEAB  MODELS  PfiOCEDUBE 

•WET'S  STUDEITIZED  BAMGE  (BSD'  TEST  FOE  TAB! ABLE:  DEIEBE 
MOTE:  THIS  TEST  COITBOIS  THE  TYPE  I  EPEBIMEMTMISE  EBBOE  BATE. 

BP  GEME BALLY  HAS  A  HIGHER  TYPE  II  EBBOR  BATE  THAM  KEQMQ 

ALPHAS  05  DF*)6  M$E*334.60*r 
CEITICAL  VALUE  OF  STUDEKTIZEI-  BAM3EM.333 
M1HIMU*  S3G1IFICAJT  D1FFEBEICES15.851 

CAMS  WITH  THE  SAME  IETTEB  ARE  ICT  SIGS1FICATLY  DIFFERED 


GBCUFIIG 

MEAM 

V 

TBO 

A 

38  476 

25 

10 

A 

A 

37.063 

25 

8 

A 

B 

A 

26.033 

25 

8 

B 

B 

C 

20.875 

5 

4 

c 

c 

5  950 

25 

2 

:96S  15 
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SAS  12  y  TrSEI-AV  BAY!,  19B?  it 

GENERAL  LIKAS  urns  PROCEDURE 

mrr*  stupentizep  raws  'esd>  test  for  variable  pesers 

IOTI  THIS  TEST  CONTROLS  THE  TYPE  !  EXPERIMENT*! SE  ERROR  RATE. 

BITT  GENERALLY  HAS  A  fi!GHEE  TYPE  I!  ERROR  RATE  THAN  HGW 

ALPHAS. 05  DF*lf  MSEM0.8W 
CRITICAL  VALUE  Of  STUDEATIZED  RANGED.  333 
MINIMUM  SIGNIFICANT  DIFFERENCE1 5. 5286 

CANS  WITH  THE  SAME  LETTER  AH  NOT  SIGNIFICANTLY  DIFFERENT. 

TOKEN  GROUPING  MEAN  K  TB: 

A 
A 
A 
A 
A 
A 

B  A 

B 


41.435  25  10 

41.05B  2f  f 

37.21$  25  6 

!5  94?  25  4 

v;  3?7  25  2 
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SAS 

SEEHAL  IVRISi  MOOELS  PROCEDURE 

TREY'S  STUOEPIZE!1  RAMGE  (BSD'  TEST  FOR  VARIABLE-  COST 
IOTE-  THIS  TEST  WOOLS  THE  TT?E  I  EPERIKEOISE  ERROr  RATE. 
BUT  SEVERALLY  HAS  A  HIGHS?  TYPE  I!  ERROR  RAT?  THAS  REG* 

ALPHA* 0.05  DF=64  KEM64  553 
CRITICAL  VALUE  OF  STUl-tlTIZEI  RDIGE=2.S70 
l!V!M»  SIGIIFICAJT  J>!FFEl£JICE*n.  JH 

cars  with  tfi  same  letter  are  mot  sigkficaply  l-iffesep. 

7'JKEY  GROUP  IMG  CAM  K  SET 


60.1*7  2?  4 

*6  95f  25  3 

45.910  25  2 

45.162  25  i 

42.614  25  5 
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S  AS  '.3:2'  Tv  ESDI!  MAT  2.  ;y89 

GEIERAL  L1IEAR  MODELS  PEOC£DOBI 

TDKEY’S  STUDEKTIZED  EAliQE  (ESDI  TEST  FOR  VARIABLE-  SHORT 
IOTE:  THIS  TEST  COITROLS  THE  TYPE  I  EXPEBIMEVTWISE  ERROR  RATE. 

BIT  GENERALLY  EAS  A  HI3HEE  TYPE  II  ERROE  RATE  TEAR  REGNO 

ALPBA=0.O5  DF-64  MSE*84.5831 
CRITICAL  VALUE  OF  STUDEHTIZEl  RAKGEO.970 
MlVltfim  SIGIIFICUT  DIFFER£RCE^.30ifi 


ICAIiS  WITH  TEE  SAME  LETTER  ARE  ROT  SIGIilFICAVTLT  DIFFEREMT 


GbO’JPI  MG 

ME  All 

K 

SET 

A 

A 

;5.380 

25 

3 

A 

A 

It  W 

25 

4 

A 

A 

1C. 026 

25 

2 

A 

A 

15.005 

25 

5 

A 

14.424 

25 

1 
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sas  12. r  tteelay  may  :  .«*  r? 

GEKfiA-  LIKAS  KOL'ELS  ?&0 CC38E 

TUXEY'S  S?'JT*r:2ir  RASGE  :HST  TEST  FOE  ?A2IA£LEr  FEES83B? 

SOT::  THIS  TETT  C0T5K.S  TEE  TT?E  1  EKEEIOHJISE  KS01  SATE. 

BYT  GSliEBALLY  HAS  A  SlOEEE  TYPE  II  **£0E  EAT:  TEAL  S£G« 

A1?HA= S. 05  $F«6«  £E*2fS.fOS 
valve  or  5T::ei.t::e:  eaj 
k:s:o  sissmcAsr  ::?Fr-E!:::=:..;!4 

CAS:  Sir-  THE  £AC  LETTEE  AEE  i'CT  £  I Z'.:'.  VAST:  Y  1 .771111: 

tyke?  se::*:sj  kas  e  -e: 

A  St.CS  if  ? 

A 

P  A  2..SM  2f  4 

B  A 

B  A  27.435  2*  2 

B 

B  22. AM  2f  3 

B 

B  Ife.Wfc  If  ! 
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s as  \i-.r  Tins:  ay  hay  : 

GEIEiAL  LI  mi  MODELS  PBOCETDBE 

TOT'S  STU&EITIZH  EAlGE  IESM  TEST  FOS  VARIABLE:  DELEER 
I0TE:  THIS  TEST  COITROLS  THE  TTPE  1  EXPERIMEXTMISE  EBSOR  BATE. 

6CT  GEIEiALLY  HAS  A  HIGKL*  TTPE  II  ERROR  SATE  THAA’  RE3MC 


ALPHA1 0.05  DF*64  MSE=250.675 
CRITICAL  VALUE  OF  ST’JDEITIZEI  BAUGH* 3. 970 
Kill  MUM  SIGI1FICAIT  DIFFE£EMC£= 12 . 575 


MEANS  MITE  THE  SAME  LETT: 

TUKEY  GROUP  I IG 

A 

A 

B  A 

B  A 

B  A 

B 

B 

C 


ARE  HOT  SIGKIFICAATLY 
MEAL  K 

36.232  25 

22.114  25 

2? . 656  25 

21.555  25 

6.862  25 


Mr  FERE  IT. 
SET 


4 

2 

1 

5 
2 


:&S9  20 
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GEIERAl  LIRAS  MODELS  PEOCEDCBE 

TUKEY'S  STUDEJTTIZED  EAKjE  (ESD)  TEST  FOE  VAEIASLE:  PEERS 
IOTE:  THIS  TEST  C0ITB01S  TEE  TYPE  1  EXPEEIMESTWISE  E&EOE  SATE. 
BUT  GEIEBALLY  BAS  A  El GHEE  TYPE  II  EEEOE  SATE  THAI  HEWQ 

ALPHA* 0. 05  DF*64  MSE=29.3006 
CEITICAL  VALUE  OF  STUDEMT1ZED  EAIGE*3.970 
MIKIMUM  SIGSIFICAT  BIFFEBEMCEM.2976 

MEAKS  MITE  THE  SAIE  LETTEE  ABE  HO?  SIGAIFICAHTLY  DIFFESEIT. 


GSCUFIIG 

ME  AS 

I 

SET 

A 

52.717 

25 

2 

B 

39.1)4 

25 

4 

B 

B 

36.155 

25 

5 

C 

29.606 

25 

1 

C 

C 

27.502 

25 

3 

Mi 
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SA£ 


12:2*  TCESi'IY  Itt?  j.  1966  22 


uiivasute 


VARIABLE1  TEES  ID 


13  BEITS 

00IITIL££!DEF=4I 

HTEEIES 

I  125  SUB  BQTS 

125 

10Q1  MAX 

46434  Oil 

44.3286 

LOBES? 

HIGHEST 

CAN  6.472E-1J  SUB 

8.091E-11 

752  03 

i. 45969  951 

25.2995 

-45.9912 

26.3648 

STD  DEV  15.488!  VASIAHCE 

239.821 

50:  BED 

0.729511  90X 

17.6527 

-45  6433 

3.. 0922 

SKEWESS  -0.25 0999  MITOSIS 

1.43372 

252  01 

-8.99172  101 

-16.2701 

-43.9996 

34.0022 

USS  29737.8  CSS 

29737.6 

01  Mil 

-45.1912  51 

-26.1132 

-39.847 

38.3373 

CV  99999  STD  EA1 

1.38512 

11 

-45.97fi‘r 

•21.444 

46.434 

T:KEU=0  4.8735-13  PSOB)' 

1 

EAIGE 

92.4252 

SGI  iAIX  84.5  PE03>‘S‘ 

0.8747 

03-01 

17.4514 

UK  “«  0  225 

BODE 

*45.9912 

0:808*81.  0.0755447  PBOBID 

0.08 

STEM  LEAF 

• 

B0IPL0T 

KORBAL  PROBABI 

-ITT  PLOT 

4  0 

1 

0 

47.5* 

• 

3  0 

1 

0 

3  14 

2 

2  68 

2 

**♦ 

2  0134 

4 

1  566578880 

0 

1  011123444 

9 

•»** 

0  55556677777788899 

17 

♦ . 

•  Ml 

0  1111112222223344444 
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